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ABSTRACT 


This  report  presents  the  results  of  a  study  of  the  combined  effects  of 
nose  bluntness  and  real-gas  phenomena  on  blunted  slender  body  flows,  with 
particular  emphasis  on  the  flows  past  spherically  blunted  cones.  One  of  the 
primary  objectives  of  the  study  has  been  to  assess  the  effects  of  not  matching 
all  the  similitude  parameters  in  hypervelocity  test  flows  about  the  slender 
bodies  that  are  under  consideration.  In  the  first  phase  of  the  study,  direct 
comparisons  between  the  flight  and  wind-tunnel  flow  fields  are  made  for  a 
blunted  slender  cone  for  several  flight  cases  by  using  a  stream-tube  computer 
program.  In  this  part  of  the  study  the  pressure  distribution  on  the  surface  of 
the  cone  is  assumed  to  be  insensitive  to  nonequilibrium  effects.  In  the  second 
phase  of  the  study,  the  influence  of  real-gas  effects  on  blunted  slender  body 
flow  fields  is  studied  analytically  by  using  a  thin- shock-layer  approach.  It  is 
found  that  real-gas  effects  on  the  pressure  distribution  and  the  shock  shape 
for  a  blunted  slender  cone  are  smail  for  the  range  of  free-stream  conditions 
of  interest.  The  implications  of  various  existing  similitude  laws  to  hyper¬ 
velocity  testing  are  discussed  briefly,  and  a  restricted  similitude  applicable 
to  real-gas  flows  over  blunted  slender  bodies  is  derived. 
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SECTION  I 
INTRODUCTION 

In  a  recent  study  of  hypervelocity  flows  and  similitudes^  it  was  pointed 
out  that,  m  the  combined  presence  of  nonequilibrium  flow  and  nose  bluntness 
effects,  similitude  between  two  slender  bou/  flows  is  theoretically  impossible. 
That  is,  the  similitude  requirements  dictate  complete  duplication  of  ambient 
thermo-chemical  state,  free- stream  velocity  and  vehicle  size.  However,  at 
high  flight  velocities  and  low  altitudes  it  is  impossible  to  duplicate  the  ambient 
thermo-chemical  state  in  existing  wind  tunnels.  In  view  of  this  state  of  affairs, 
it  is  necessary  that  the  similitude  requirements  for  hypervelocity  flows  be 
carefully  examined  from  the  standpoint  of  delineating  the  more  important 
similarity  parameters  and  determining  the  degree  to  which  the  less  important 
similitude  requirements  can  be  relaxed.  For  the  hypervelocity  flight  regime 
considered  in  the  study  reported  in  Ref.  1,  it  was  found  that  there  do  not  exist 
sufficient  experimental  data  or  theoretical  solutions  tc  evaluate  the  sensitivity 
of  various  flow  quantities  to  variations  in  the  similitude  parameters  or  to 
fully  determine  whether  some  parameters  are  less  important  than  others. 

Thus,  the  purpose  of  the  present  investigation  has  been  to  obtain  new  theo¬ 
retical  solutions  for  blunted  slender  bodies  with  the  specific  objective  of 
studying  the  effects  of  relaxing  certain  similitude  requirements. 

There  are  two  phases  to  the  present  work.  The  first  phase  is  con¬ 
cerned  with  a  direct  comparison  of  flight  and  wind-tunnel  flow  fields  about  a 
blunted  slender  cone.  For  this  purpose,  the  hypothetical  wund-tunnel  per¬ 
formance  defined  in  Ref.  1  has  been  used.  The  comparisons  are  made  at 
s*  rai  flight  conditions  for  which  complete  duplication  can  not  be  achieved 
in  the  wind  tunnel,  i.  e.  at  low  altitudes.  Partial  results  of  the  study  have 
been  reported  in  Ref.  2.  The  second  phase  of  the  work  is  concerned  with 
flight  conditions  at  higher  altitudes  where  it  is  assumed  that  a  wind  tunnel 
could  fully  duplicate  the  ambient  altitude  thermodynamic  state  if  there  were 
no  nozzle  noneauiliorium  effects.  Here  the  emphasis  is  on  an  investigation 

A  »  © 

of  r  onequilibrium  flows  about  blunted  slender  bodies  within  the  framework  of 
thin  shock  layer  theorv  ar.d  the  blast-wave  analogy.  Both  studies  are  restricted 
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to  bodies  at  zero  angle  of  attack. 

In  the  following  Section  (Sec.  II)  of  this  report,  the  flight  regime  of 
interest  will  be  defined  and  related  to  wind-tunnel  performance  capabilities. 
Because  of  the  limitations  of  existing  wind  tunnels,  a  hypothetical  wind-tunnel 
performance  is  defined  which  minimizes  nozzle  nonequilibrium  effects  over 
most  of  the  flight  region.  Test  conditions  in  the  hypothetical  wind  tunnel  are 
compared  with  flight  conditions  for  various  velocities  and  altitudes.  These 
comparisons  include  estimates  of  the  effects  of  nonequilibrium  phenomena  on 
the  nozzle  expansion. 

In  Sec.  Ill,  some  of  the  general  features  of  the  flow  about  blunt-nosed 
slender  bodies  are  discussed  with  particular  attention  to  rarefaction  effects, 
viscous  interactions  and  flow  chemistry.  This  is  followed  by  a  brief  review 
(Sec.  IV)  of  the  similitudes  for  such  flews.  The  next  two  sect-~ns  (Sec.  V,  VI) 
are  concerned  with  the  two  studies  carried  out  as  part  of  this  investigation. 
The  concluding  section  summarizes  the  general  findings  of  these  studies. 
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SECTION  II 

COMPARISON  OF  WIND  TUNNEL  FLOWS  WITH  FLIGHT  CONDITIONS 

In  this  section  the  flight  region  of  interest  will  be  defined  and  related 
to  wind  tunnel  performance  requirements.  For  the  purposes  of  later  compar¬ 
isons  of  flight  and  wind-tunnel  flow  fie'ds,  a  hypothetical  wind  tunnel  perfor¬ 
mance  is  postulated.  Test  section  conditions  in  this  hypothetical  wind  tunnel 
are  compared  to  flight  conditions.  The  effects  of  nonequilibrium  nozzle  flow 
on  wind-tunnel  test  conditions  are  estimated.  Finally,  the  test  conditions 
available  in  present-day  wind  tunnels  are  compared  with  flight  conditions  for 
several  velocities  and  altitudes. 

The  hypervelocity  flight  regime  of  interest  in  this  study  is  the  altitude 
range  from  50,  000  ft.  to  250,  000  ft.  for  velocities  from  IQ,  000  fps  to  36,  000 
fps  (see  Fig.  1).  This  regime  encompasses  all  of  the  continuous  flight  corridor 
lying  above  10,  000  fps  velocity  and  below  250,  000  ft.  altitude  and  all  of  the 
reentry  corridor  for  manned  satellites  and  lunar  vehicles  (except  portions  of 
skip  trajectories).  These  flight  regions  are  also  shown  in  Fig.  i.  This 
altitude-velocity  map  further  serves  to  define  the  range  of  values  for  various 
similarity  parameters  and  will  be  referred  to  frequently  in  subsequent  sections. 

The  severe  requirements  that  duplication  of  ambient  altitude  conditions 
at  hypersonic  speeds  imposes  on  a  wind  tunnel  are  illustrated  in  Fig.  2.  This 
altitude -velocity  map  shows  the  nozzle  reservoir  pressures  and  temperatures 
necessary  for  duplicating  flight  conditions  assuming  an  isentropic  expansion 
of  real  air  in  thermodynamic  equilibrium.  It  should  be  noted  that  there  are 
few  facilities  presently  operating  at  reservoir  pressures  greater  than  2000 
atm;  therefore,  reservoir  pressures  above  4000  atm  have  not  been  shown  in 
Fig.  2.  In  existing  facilities  velocities  greater  than  20,  000  fps  can  be 
obtained,  but  pressure,  temperature  and  density  for  altitudes  below  250,  000  ft. 
cannot  be  duplicated  simultaneously. 

In  order  to  assess  the  magnitudes  of  the  mismatches  that  can  be  expected 
in  a  wind  tunnel,  it  was  necessary  to  specify  a  wind-tunnel  performance  cap¬ 
ability.  This  was  done  with  due  regard  to  nozzle-flow  nonequilibrium  effects. 

In  descriptions  of  recently  proposed  new  facilities,  several  authors^’ ^ 
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adopted  a  maximum  entropy  criterion  for  minimizing  nonequilibrium  effects 

in  the  nozzle  flow.  This  criterion  is  based  on  the  fact  that  the  nonequilibrium 

effects  on  the  expansion  of  air  from  a  stagnant  reservoir  to  high  velocity  can 

5  6 

be  correlated  in  terms  of  entropy  *  .  While  specifying  a  maximum  entropy 
of,  say,  S/R  =  31.  5  implies  that  the  nonequilibrium  effects  will  be  small  in  the 
nozzle  flow,  it  also  limits  the  maximum  altitude  at  which  full  duplication  of 
ambient  conditions  can  be  achieved  to  190,  000  ft.  The  conclusion  is  that,  if 
air  is  expanded  in  a  nozzle  from  a  high  temperature,  high  pressure  reservoir, 
nonequilibrium  effects  cannot  be  eliminated  when  testing  at  high -altitude, 
hypervelocity  conditions. 

The  hypothetical  wind-tunnel  capability  postulated  in  the  present  study*’ ^ 
is  defined  by  the  full-duplication  boundary  line  shown  in  Fig.  2.  Above  this 
full -duplication  boundary,  it  is  assumed  that  both  flight  velocity  and  ambient 
altitude  conditions  could  be  duplicated  in  the  test  section  if  the  nozzle  expansion 
remained  in  thermochemical  equilibrium.  Below  this  boundary  line  it  is  assumed 
that  the  flight  velocity  and  only  one  of  the  ambient  altitude  conditions  (e.  g.  , 
density-altitude  or  free-stream  temperature)  can  be  duplicated.  In  Ref.  1  it 
was  shown  that  nozzle  nonequilibrium  effects  on  test  section  velocity  and  density 
are  negligible  below  the  full-duplication  boundary  for  this  hypothetical  wind 
tunnel.  Estimates  of  tne  nonequilibrium  effects  on  pressure,  temperature  and 
Mach  number  in  this  partial  simulation  region  are  presented  in  this  section. 

Above  the  full -duplication  boundary  of  Fig.  2,  duplication  of  velocity 
and  ambient  altitude  thermochemical  state  is  compromised  by  nozzle  flow  non¬ 
equilibrium.  The  effect  on  velocity  was  indicated  in  Ref.  1.  Nonequilibrium 
effects  on  the  test  conditions  will  be  discussed  in  greater  detail  later  in  this 
section. 


For  the  flow  phenomena  that  are  of  interest  in  the  present  study,  the 
more  important  parameters  to  be  duplicated  in  the  partial  simulation  region 
are  the  density-altitude  {i.e.  ,  the  free-stream  density)  and  the  flight  velocity. 
These  two  parameters  insure  the  duplication  of  the  fluxes  of  free-stream  mass, 
momentum  and  energy  provided  that  A^sinO-  (where  <T  is  the  shock  angle) 
is  sufficiently  large  so  the  bow  shock  wave  is  everywhere  strong. 
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In  Section  V  a  comparison  is  made  of  flight  and  wind-tunnel  flow  fields  about 


a  blunted  slender  cone.  The  three  flight  conditions  considered  in  that  study 
are  velocities  of  15,  000  fps  and  25,  000  fps  at  a.n  altitude  of  80,  000  ft.  and 
35:  000  fps  at  180,  000  ft.  altitude.  It  can  be  seen  in  Fig.  2  that  these  con¬ 
ditions  are  substantially  beyond  the  duplication  capabilities  of  the  hypothesized 
wind  tunnel.  From  the  performance  capability  hypothesized  in  Fig.  2,  it  is 
possible  to  determine  the  wind-tunnel  test  conditions  when  the  flight  velocity 
and  ambient  density  altitude  are  duplicated  under  the  assumption  that  the  nozzle 
expansion  is  isentropic  and  in  thermochemical  equilibrium.  The  validity  of  the 
isentropic-equilibrium  assumption  is  then  evaluated.  The  flight  and  wind-tunnel 
conditions  and  related  parameters  for  the  three  cases  are  compared  in  the 

following  tables  for  the  case  of  an  equilibrium  nozzle  flow. 

Table  II- I 

COMPARISON  OF  FLIGHT  AND  HYPOTHETICAL  WIND  TUNNEL  CONDITIONS 
Case  I:  Un  =  15,000  fps,  80,000  ft  Alt. 


Table  H-2 


COMPARISON  OF  FLIGHT  AND  HYPOTHETICAL  WIND  TUNNEL  CONDITIONS 


Case  U:  Um=  25,000  fps,  80,000  ft  Alt. 


PARAMETER 

FLIGHT 

WIND  TUHNEL 

VELOCITY,  64,  (fos) 

25,000 

25,000 

DEHSITY,  /4,  (ama) 

3.420xl0"2 

3. 420x1 0'2 

PRESSURE,  ^  (atm) 

2.765xl0"2 

1. 513x10*’ 

TEMPERATURE,  T  (eK) 

220.9 

1209 

HACH  HUMBER,  M 

25.57 

11.25 

H 

TOTAL  EHTHALPY,  (°K) 

1.019xl05 

!.056xl05 

ENTROPY,  S/R 

26.44 

30.99 
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Table  U-3 

COMPARISON  OF  FLIGHT  AND  HYPOTHETICAL  WIND  TUNNEL  CONDITIONS 
Case  HI:  Ua=  35,000  fps,  180,000  ft  ALT. 


PARAMETER 

FLI6HT 

HIND  TUNNEL 

VELOCITY,  Ua  (fps) 
DENSITY,  (tat) 

PRESSURE,  Pa  (ato) 
TEMPERATURE.  Ta  (*K) 

MACH  HUMBER,  M 

CO  Jy 

TOTAL  ENTHAl  PY ,  ( *  K) 

EHTR0PY,  S/R 

35,000 

4.tH3xlO‘‘4 

9.292xl0”4 

265.9 

32.63 

1 .990x10® 

31.25 

35,000 

N.1l3xl0'4 

6.l35xlO*4 

380.0 

27. 3« 

l .9Sdxl05 

32.15 

It  is  seen  that,  in  the  first  two  cases,  appreciable  differences  in  free- 
stream  pressure,  temperature  and  Mach  number  exist  between  the  flight  and 
wind-tunnel  conditions.  Along  the  lower-boundary  of  the  altitude-velocity 
map  in  Fig.  2,  the  mismatch  in  static  temperature  can  be  as  great  a  factor 
of  7  for  the  wind-tunnel  performance  hypothesized  here.  This  emphasizes 
the  importance  of  understanding  the  effects  of  partial  simulation. 

Because  the  free- stream  conditions  in  the  wind  tunnel  tabulated  above 
have  been  computed  under  the  assumption  of  an  isentropic,  equilibrium  nozzle 
flow,  it  is  pertinent  to  investigate  the  effects  the  nonequilibrium  phenomena 
would  have  on  these  conditions.  Lordi  and  Mates^  have  shown  that  for  a  given 
nozzle  scale  and  geometry  the  nonequilibrium  effects  on  composition  and 
thermodynamic  state  of  the  gas  could  be  correlated  with  reservoir  entropy. 
Thus,  their  results  can  be  used  to  estimate  the  nonequilibrium  effects  for  the 
present  wind-tunnel  conditions.  These  estimates  are  summarized  in  the 
following  table. 

Table  H-4 

NOZZLE  NONEQUILIBRIUM  EFFECTS  IN  HYPOTHETICAL  WIND  TUNNEL 


CASE  I 

CASE  n 

CASE  HI 

FREE-STREAM  PARAMETER 

U  =  25,000  fps 

OS 

U =  35,000  fps 

CD 

(  UJ  NOHEQ/fCfJ  EQ 

1.00 

0.995 

0.793 

(*„>  H0NE */(*£,)  ^ 

1.00 

0.815 

0.700 

(  Tm)  H0MEQ/ (  )  EQ 

1.00 

0.790 

0.670 

(  MJ  N0HEQ/{  MJ  EQ 

1.00 

1.10 

1.21 

FREE-STREAM  composition 

{HOLES / gm  OF  HIXTURE) 

NITR0GEH,  H2 

2.6xl0'2 

2.6xl0'2 

2.6xl0”2 

oxygen,  o2 

5.3xl0*3 

5.3xl0'3 

5.0xl0*3 

NITRIC  OXIDE,  NO 

3.3xl0'3 

2.9xl0*3 

2.6xl0'3 

ATOMIC  OXYGEN,  0 

1 .7xl0~4 

7.7xl0'4 

1.7xl0'3 
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Even  when  nonequilibrium  effects  are  energetically  unimportant  as  in 
the  above  three  cases,  the  gas  composition  can  be  altered  from  that  of  undis¬ 
sociated  air,  particularly  through  the  formation  of  nitric  oxide,  as  is  shown 
in  Table  II  -  4.  The  fact  that  the  wind-tunnel  flow  composition  is  altered  by 
the  amounts  shown  in  Table  II  -  4  will  not  significantly  affect  the  flow  about 
a  body  if  the  flow  achieves  thermochemical  equilibrium  after  passing  through 
the  bow  shock  wave.  This  is  illustrated  in  Sec.  V  where  the  equilibrium 
stagnation- point  conditions  are  presented  for  the  three  foregoing  cases 
(Tables  V  -  1,  V  -  2  and  V  -  3).  A  similar  result  is  shown  in  Ref.  7.  Further 
more,  the  effect  of  nozzle-flow  nonequilibrium  on  free-stream  pressure, 
temperature  and  Mach  number  is  to  decrease  the  differences  between  wind 
tunnel  and  flight  conditions  shown  in  Tables  II  -  1,  -  2,  and  -  3. 

Next,  we  consider  the  effects  of  nonequilibrium  nozzle  flow  in  the  region 
above  of  the  full-duplication  boundary,  Fig.  2.  For  convenience,  we  will  con¬ 
sider  just  the  five  velocity-altitude  conditions  listed  in  the  following  table. 


Table  II-5 

VELOCITY- ALTITUDE  CONDITIONS  IN  FULL-DUPLICATION  REGION 


CASE 

VELOCITY  (f ps ) 

ALTITUDE  (ft) 

IT 

15,000 

180,000 

T 

25,000 

180,000 

ZE 

15,000 

250,000 

HI 

25,000 

250,000 

znr 

35,000 

250,000 

Using  the  frozen-enthalpy  correlation  with  entropy  giver,  in  Ref.  5  and 
in  addition  other  correlations  derivable  from  the  results  presented  there, 
estimates  have  been  made  of  the  nonequilibrium  effects  at  the  above  velocity- 
altitude  conditions.  The  results  are  given  in  the  following  table. 
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Table  H-6 


NOZZLE  NONEQUILIBRIUM  EFFECTS  IN  HYPOTHETICAL  WIND  TUNNEL 


FREE-STREAH  PARAMETER 

CASE  n 

CASE  I 

CASE  n 

CASE  HI 

CASE  HIT 

(4o)  N0NEQ/  (  Ueo)  EQ 

0.977 

0.992 

0.945 

0.980 

0.990 

(  P»)  NONEQ/^®)  EQ 

0.770 

0.770 

0.550 

0.550 

0.550 

(  Ta> )  NONEQ/f7® )  EQ 

0.758 

0.758 

0.506 

0.506 

0.506 

(Ma>)  KOHEQ/^®')  EQ 

I.U 

1.13 

1.29 

1.34 

1.36 

FREE-STREAH  COMPOSITION  (M0lES/gm  OF  ORIGINAL  AIR) 

NITROGEN,  N2 

2.57xl0“2 

2.57xl0“2 

2.58xl0“2 

2.58xl0"2 

2.58xl0"2 

OXYGEN,  02 

5. 45x1 O'3 

5.45xl0“3 

4.32xl0"3 

4.32xl0”3 

4.32x10~2 

NITRIC  OXIDE,  NO 

2.54x10"3 

2.54x10‘3 

2,33x10"3 

2.33x10’3 

2.33x!0“3 

ATOMIC  OXYGEN,  0 

1.03xl0"3 

1.03xl0‘3 

3.53xl0“3 

3.53xlO“3 

3. 53x1  O'3 

The  effects  of  nozzle -flow  nonequilibrium  on  free -stream  velocity  (and 
5 

hence  on  density  also  )  are  not  large  for  any  of  the  above  cases  and  could  be 

suitably  compensated  by  adjustment  in  the  total  enthalpy.  On  the  other  hand, 

appreciable  effects  on  pressure  and  temperature  occur,  particularly  at  the 

higher  altitude  (Cases  VI  -  VIII).  However,  the  resulting  mismatches  are 

not  as  severe  as  those  occurring  at  80,  000  ft.  due  to  tunnel  performance 

limitations.  The  most  serious  nonequilibrium  effect  for  the  above  conditions  is 

the  alteration  of  the  test-gas  composition,  particularly  at  the  higher  altitude 

where  it  is  unlikely  that  thermochemical  equilibrium  will  exist  behind  a  bow 

shock  wave.  Thus,  the  existence  of  free-strearr.  dissociation  will  influence 

8 

the  flow-field  chemistry  about  a  model.  Inger  has  investigated  the  influence 

of  a  dissociated  free -stream  on  the  flow  in  the  stagnation  region.  However, 

his  method  applies  only  for  highly  nonequilibrium  free -stream  conditions  and 

does  not  permit  isolating  the  effects  of  free- stream  dissociation  from  mismatches 

9 

in  free-stream  Mach  number.  On  the  other  hand,  Gibson  has  suggested  a 
‘subtraction  rule"  for  accounting  for  the  effects  of  free-stre^m  dissociation 
separate  from  other  effects. 

Finally,  it  is  of  interest  to  compare  test  conditions  that  might  be 
attainable  in  present-day  wind  tunnels  with  flight  conditions  for  the  velocity- 
altitude  cases  listed  in  Table  II  -  5.  For  this  purpose  we  select  a  maximum 
reservoir  pressure  capability  of  2000  atm  and  assume  that  reservoir  temperatures 
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up  to  17,  000  °K  are  available.  This  would  provide  flow  velocities  up  to 
36,  000  fps  for  isentropic,  equilibrium  nozzle  flow. 

Inspection  of  Fig.  2  indicates  that,  in  the  absence  of  nonequilibrium 
effects,  such  a  wind  tunnel  would  provide  full  duplication  at  both  altitudes 
for  XJ^  =  15,  000  fps  (Cases  IV  and  VI).  The  nonequilibrium  effects  on  test 
conditions  for  these  two  cases  would  be  the  same  as  listed  in  Table  II  -  6 
for  our  hypothetical  wind  tunnel.  For  the  remaining  cases,  the  flight  and 
wind-tunnel  conditions  are  compared  in  the  following  tables.  Since  full  dupli¬ 
cation  is  not  possible  for  these  cases,  the  wind-tunnel  is  assumed  to  operate 
at  pc  =  2000  atm  and  to  duplicate  free-stream  velocity  and  density.  Here, 
the  wind-tunnel  conditions  are  based  on  an  isentropic,  equilibrium  nozzle 
flow.  The  nonequilibrium  effects  will  be  discussed  subsequently. 


Table  11-7 

COMPARISON  OF  FLIGHT  AND  2000-AIM  WIND  TUNNEL  CONDITIONS 
Case  V:  Um  =  25,000  fps,  180,000  ft  Alt. 


PARAMETER 

FLIGHT 

WIND  TUNNEL 

VELOCITY,  64,  (fps) 
DENSITY,  ^(ama) 
PRESSURE,  (atm) 

TEMPERATURE,  ^(“X) 

MACH  NUMBER,  M 

CO 

TOTAL  ENTHALPY,  HjR°* 
ENTROPY,  S/R 

25,000 

4.41 3xl0-4 
4.292xl0-4 

265.9 

23.30 

1.020xl05 

31.25 

25,000 

4. 413X10"1* 
4.161x1 0— 3 

2475 

8.325 

1 . 1 37  x  1 05 

38.77 
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Table  i_L— 6 

CCXPASiSOH  OF  FLIGHT  AMD  2000- A1H  W1HD  TUKKEL  COHDmGXS 
Case  ill:  U^.  =  25,000  fps,  250,000  ft  Alt. 

PmjgTiS  |  FLiST  tfiKD  TUB  EL  | 

TElCCIlr,  £4(f?»)  J  25,009  25.  KQ « 

D3SITT.  ^{ts*}  I  2.307x50"=  2.S07xl0~:>  I 

?2tSSU2£.  C  (itiij  *  2.CC7xlG'5  1.233x50"*  ' 

»  j  , 

TEX?£2JttU2r.  r!*lj  i  535.5  !  i205 

O  •  #  |  [ 

XiC3  JTO»e£2.  M  27. 15  t  10.35  ■ 

*  j  -  i 

TOTAL  OrtSALPr .*,:£{*?)  1.013x10=  j  1.055x10^  j 

D!T£0?T.  S/5  33.2?  1  23.05 

1  I 


Table  U-3 

CCHPAR.SOH  OF  FLIGHT  AMD  2GOO-ATH  WIND  TUNNEL  C0HDIT10HS 


Case  vill:  L/^  -  35.0G0  fps,  250.000  ft  Alt. 


Clearly  the  mismatches  in  pa,,  and  Ma>  are  large  for  these 

three  cases.  The  question  is  whether  they  are  too  large  to  permit  simulation 
of  a  Mach  number  independent  flow.  The  simplest  test  of  this  is  to  compare 
icuilibrium  and  frozen  flow  conditions  behind  a  normal  shock  wave.  The 
equilibrium  conditions  for  the  flight  case  are  given  in  Ref.  10.  For  the  2000-atm 
wind  tunnel  they  have  been  calculated.  The  equilibrium  flow  -results  are  compared 
in  the  following  table. 
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Table  j_l-!G 


CCHPAR! SOX  OF  E0U1 LI BRfUM  COXDITIOXS  3EHIH9 
A  XORKAL  SHOCK  WAVE 
IX  FLIGHT  AXD  2000- ATH  W1XD  TUHHEL 

PARAMETER  j  FLIS3T  I  WIKD  TUKXEl 


CASE  7:  -  25.000  fps  ISO. COO  F7  AITI7U3E 


PRESSURE,  p  (ala) 
DESS1TY.  pr  (a«a) 
T3<?E2ATU2£,  7;(*A) 

0.3072 

7.273xl0*? 

£801 

0.2053 

6.751x10“ 

6952 

_ 

j  CASE  ZH:  =  25.000  fps  25C. COO  FT  ALTITUDE 

PRESSURE.  ^  (aln) 
DEKSITT.  /Ds  (aaa) 
TEMPERA IURE.  7^(*K) 

0.01957 

5. 196x1 0“3 

5959 

.01S35 

5.  lOlxlO"5’ 

6022 

CASE  HU:  0^  =  35.0C0  fps  250,000  FT  ALTITUDE 

PRESSURE,  ps  (ats) 
DENSITY,  y^(aaa) 
TEMPERATURE,  Ts  (’K) 

0.03837 

5.815xl0“3 

10,120 

0.03856 

5.591x10“” 

10.320 

For  frozen  flow  the  ideal-gas  shock  equations  are  used  with  the  proper  values 
of  the  specific-heat  ratio  for  the  wind  tunnel  conditions.  These  are  Y  =  1.416, 
1.401  and  1.465  for  Cases  V,  VII  and  VIII,  respectively.  The  frozen  flow 
results  are  tabulated  below. 
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Table  II-ll 

COHPAR I SON  OF  FROZEN- FLOS  CONDITIONS  BEHIND 
A  NORHAL  SHOCK  HAVE 

IH  FLIGHT  AND  2000- ATX  WIND  TUNNEL 


PARAMETER 

FLISiT  WIND  TUNNEL 

CASE  J:  U =  25,000  fps  ISO. 000  FT  ALTITUDE 

PRESSURE,  f>s  (ata) 
DENSITY,  /?s(*3i) 
TEMPERATURE,  7*  (*K) 

0.271S 

2.629X10-3 

28.310 

0.2655 

2.354xl0'3 

29.510 

j  CASE  213:  U^  =  25,000  fps  25C. 00  FT  ALTITUDE 

PRESSURE,  -e  (ata) 
DENSITY,  />s  (aaa) 
TEMPERATURE,  7^  (*K) 

0.01730 

1.673x10"’ 

28,260 

0.01935 

1 .61 2x1 O-’ 

29,420 

CASE  203:  U*,  =  35,000  fps  250.000  FT  ALTITUDE 

PRESSURE,  -ps  (ata) 
DENSITY,  ps  (aaa) 
TEMPERATURE,  7^’K) 

0.03390 

1.678x10’’ 

55,220 

0.03300 

1 .405x10”’ 

53,200 

I 
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It  is  oDservea  tnai  xor  uoui  equiiiunu 


oth  eauilibrium  and  frozen  flow  the  conditions 


behind  a  normal  shock  wave  in  the  2000-atm  wind  tunnel  correspond  closely 
to  the  flight  conditions  for  these  three  cases  when  nonequilibrium  nozzle  flow 


effects  are  neglected.  Thus  Mach  number  independent  flow  fields  would  be 
quite  well  duplicated  in  the  wind  tunnel  if  nozzle  nonequilibrium  effects  were 


not  present. 


When  the  flow  field  is  not  Mach  number  independent,  one  must  expect 
differences  as  a  result  of  mismatches  in  the  free-stream  conditions.  For 
Case  VII  (Table  II- 3)  the  relative  free-stream  differences  between  the  2000-a 
wind  tunnel  and  flight  are  nearly  the  same  as  those  for  the  high  performance 
hypothetical  wind  tunnel  in  Case  II  (Table  II- 2 ) .  Thus,  the  results  reported 
for. this  latter  case  in  Sec.  V  would  be  indicative  of  the  relative  differences 
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oetween  such  flows  in  the  2000-atm  wind  tunnel  for  Case  VII-  In  Cases  V  and 
Vffi  the  free  -  stream  pressures  and  temperatures  in  the  2000-atm  tunnel  are 
about  10  and  15  times  greater  than  the  corresponding  flight  values,  however 
nozzle  nonequilibrium  effects  will  reduce  these  differences  appreciably. 

Estimates  of  the  nonequilibrium  effects  on  the  free-stream  conditions 
in  the  2000-atm  wind  tunnel  have  been  made  using  the  entropy  correlations 
described  previously.  The  results  are  summarized  in  the  following  table. 


Table  11-12 

NOZZLE  NONEOU I  LIBRIUM  EFFECTS  IN  2000-ATM  WINS  TUNNEL 


FREE-STREAM  PARAMETER 

CASE  I 

CASE  m 

CASE  TiTT 

f  ^Jhoheq  •  (Uahq 

0.99 

0.9'4 

0.93 

(/LWeq  1  (/Ueq 

1,01 

1.06 

I.C7 

(ftj'HOMEQ  1  (wbj  EQ 

0.55 

0.18 

0.10 

i  'kJi'OHEQ  1  (^)e<5 

0.50 

0.14 

0.075 

(HJkGNEQ  / 

1.3 

2.2 

2.6 

FREE-STRE/M  COMPOSITION  (HOlES/go  OF  ORIGINAL  AIR)  J 

NITROGEN,  H2 

2.6x1 0”2 

2 

2.6x10  ^ 

1.6xl0"2 

OXYGEN,  02 

1.9xlO'5 

1x10'** 

NITRIC  OXIDE,  NO 

1x10'** 

6. 3x10"  ** 

3.5x10"® 

ATOMIC  HITROGEN,  N 

3.7x10"** 

lxlO"4 

1 .9x1  O'2 

A  TAtl  1  r  nvycpu  A 

Mivni  v  ua  i  Olii  j  v 

! .5x10” 2 

i  ii „ in- 2 

•  •  1  /\  »  V 

i  c-in-2 

•  •  JA  •  V 

When  nozzle  nonequilibrium  effects  are  accounted  for,  the  large 
differences  in  free-stream  pressure,  temperature  and  Mach  number  between 
flight  and  wind  tunnel  are  drastically  reduced.  For  example,  in  Case  VIII, 
the  wind-tunnel  Mach  number  is  increased  from  9.  11  to  23.  7  compared  to  a 
flight  value  of  38.05.  On  the  other  hand,  the  nozzle  nonequilibrium  alters 
the  test  gas  composition  in  an  undesirable  manner.  In  Cases  V  and  VII  the 
test  gas  is  primarily  a  mixture  of  molecular  nitrogen  and  atomic  oxygen; 
whereas  in  Case  VIII  it  is  composed  of  nearly  equal  parts  of  molecular  nitrogen, 
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atomic  nitrogen  and  atomic  oxygen.  Fortunately,  large  differences  in 
gas  composition  do  not  necessarily  result  in  proportionately  large  effects 
on  pressure,  density  and  temperature.  To  illustrate  this  fact,  the  limiting 
cases  of  equilibrium  and  frozen  flow  behind  a  normal  shock  wave  have  been 
computed  for  the  nonequilibrium  free- stream  conditions  of  the  2000-atm 
wind  tunnel  (Table  11-12).  The  results  are  shown  in  Table  II—  1 3. 


Table  11-13 

EOUILlBRiUH  AND  FROZEN  FLOW  CONDITIONS  BEHIND  A 
NORMAL  'SHOCK  WAVE  IN  A  N0NE0UI LIBRIUM  2000- ATM  WIND  TUNNEL 


PARAMETER  J  EQUILIBRIUM  j  FROZE*  FLOW  | 

CASE  Y:  'J  =  25.000  fps  180,000  FT  ALTITUDE 


1  PRESSURE,  (ato)  i  C.3C33  |  0.2446  | 

1 DEKSITY,  [sat)  '  6.tH2xlO"3  !  2-23*lxl0"3 

1  TEMPERATURE,  7j  ['K)  [  7135 _ j  26,400 

j  CASE  UL:  Um  =  25,000  fps  250,000  FT  ALTITUDE 


!  PRESSURE.  (ato)  0.01828  0.01606 

(  DEKSITY,  J0  (aaa)  4.8l7xl0"4  1.565x10"’ 

TEHPERATURE,  Ts  (’K)  5952  23,120 

CASE  HU:  U^  =  35,000  fps  250,000  FT  ALTITUDE 

r RESSURE,  p  (atra)  0.03495  0.03048 

DEKSITY,  p  (ana)  4.095x10"^  1.384xl0_t* 

TEHPERATURE^,  T$  (*K)  10,200  40,020 

Comparison  of  the  equilibrium-flow  limit  with  Table  II- 10  shows 
relatively  little  effect  of  the  dissociated  free-strearn  on  conditions  behind 
the  shock  wave.  The  differences  from  the  flight  case  are,  however, 
somewhat  larger  than  when  nozzle  nonequilibrium  effects  were  not  taken 
into  account. 


Comparing  Tables  II- 1 1  and  11-13  indicates  that  in  the  frozen-flow 
limit,  the  effects  of  nozzle  flow  nonequilibrium  are  more  pronounced  than 
in  the  above  equilibrium  limit  -  particularly  on  the  temperature  Ts .  Although 
a  direct  comparison  with  Ref.  8  is  difficult,  these  results  are  generally  in 
qualitative  though  not  quantitative  agreement  with  Inger's  findings. 
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SECTION  III 

GENERAL  FEATURES  OF  BLUNT-NOSED  SLENDER  BODY  FLOWS 


The  slender  cone  having  a  blunted  nose  is  of  particular  interest  because 
it  typifies  more  general  blunt-nosed  slender  bodies.  In  this  section  the  various 
features  of  the  flow  field  about  a  slightly  blunted  cone  are  discussed  for  the 
flight  region  shown  in  Fig.  i.  The  important  items  to  be  considered  are  rare¬ 
faction  effects,  viscous  effects  and  flow  chemistry.  Since  these  phenomena 
depend  on  geometry  and  scale,  the  discussion  will  pertain  to  the  specific  body- 
treated  in  Sec.  V.  This  body  has  a  cone  half-angle  of  9s  ana  a  nose  radius  of 
either  0.2  or  1.5  inches.  Where  a  length  is  needed,  a  value  of  88  inches  has 
been  assumed.  In  discussing  the  various  phenomena  listed  above,  it  is  con¬ 
venient  to  consider  the  flow  in  the  blunt-nose  region  separately  from  the  after 
body  flow  along  the  cone.  This  distinction  is  reasonable  since  the  flows  in 
these  two  regions  may  be  quite  dissimilar. 

Because  the  flight  region  of  interest  (Fig.  1)  extends  to  an  altitude  of 
250  Kft,  the  first  concern  is  about  the  possible  presence  of  rarefaction  effects. 
As  long  as  the  bow  shock  wave  is  thin  relative  to  the  shock  stand-off  distance 
at  the  nose  and  the  shock-layer  thickness  along  the  cone,  the  flow  may  be 
treated  as  a  continuum  fluid  and  rarefaction  effects  may  be  considered  unim¬ 
portant.  Hayes  and  Pro'ostein^  have  considered  the  various  rarefaction  and 
viscous  interaction  flow  regimes  and  have  established  simple  criteria  which 
will  be  adopted  here. 

In  the  region  of  the  blunt-nose,  the  mean-free -path  behind  the  shock 
serves  as  a  measure  of  the  shock  thickness  and  the  following  quantities  are 
all  of  the  same  order  of  magnitude 

A  S  A- 03  ^  Meo 

where  X3  is  the  mean-free-path  behind  the  bow  shock  wave,  A  is  the  shock 
stand-off  distance,  is  the  mean-free-path  in  the  free-stream,  R  is  the 
nose  radius  and  Rep  is  Reynolds  number  based  on  free-stream  conditions 

and  nose  radius.  A  suitable  criterion  for  continuum  flow  in  the  nose  region  is 
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that  these  ratios  be  small  compared  to  unity.  Fig.  3  shows  results  for  the 
two  nose  radii  of  interest  and  indicates  the  various  flow  regimes  defined  by 
Hayes  and  Probstein**.  For  each  Rw  ,  the  upper  bound  of  the  shaded  region 

is  defined  by  RN  and  the  lower  bound  by  /Re^  .  This  indicates 

9  • 

the  degree  to  which  these  ratios  differ.  The  boundaries  of  the  flow  regimes 
shown  in  Fig.  3  are  based  on  the  criteria  of  Ref.  1 1  with  the  interpretation 
that  1/10  is  much  less  than  one.  Clearly,  these  boundaries  are  not  distinct 
as  shown  in  Fig.  3,  but  should  have  some  width  analogous  to  the  shaded  areas 
shown  for  the  two  nose  r?uii. 


From  Fig.  3  it  is  seen  that  continuum  flow  exists  at  all  altitudes  of 
interest  for  the  larger  nose  radius  and  below  about  225  Kft  for  the  smaller 
nose  radius.  It  Should  be  pointed  out  that  in  the  incipient  merged  layer  regime 
the  shock  wave  is  not  a  thin  discontinuity  governed  by  the  Rankine-Hugoniot 
relations  even  though  the  flow  may  be  treated  as  a  continuum.  Above  225  Kft 
for  the  smaller  nose  radius,  the  continuum  flow  concept  begins  to  breakdown. 


Having  established  that  rarefaction  effects  are  generally  unimportant 
for  the  bodies  and  flight  conditions  of  interest,  we  next  consider  the  importance 
of  viscous  effects.  For  the  flow  about  the  blunt-nose,  the  regions  of  important 
viscous  effects  are  defined  in  Fig.  3.  Considering  just  the  smaller  nose  radius, 
it  is  seen  that  the  boundary  layer  is  thin  below  100  Kft.  From  100  Kft  to 
about  160  Kft  vorticity  in  the  inviscid  flow  influences  the  boundary  layer. 


r  rom  lo0  Kft  to  about  190  Kft,  the  shock  layer  becomes  fully  viscous  although 


the  bow  shock  may  still  be  considered  as  a  discontinuity.  Above  190  Kft,  the 


rarefaction  effects  on  shock  structure  begin  to  enter  as  just  described. 


Similar  viscous  interaction  effects  occur  on  slender  bodies  as  a 
result  of  the  relative  thicknesses  of  the  boundary  layer  and  the  shock  layer. 
These  viscous  phenomena  take  the  form  of  weak  and  strong  pressure  inter¬ 
actions  between  the  boundary  layer  and  the  outer  inviscid  flow  vorticity 
interaction,  transverse  curvature  effects  and  ultimately  a  fully  viscous 
shock  layer.  For  slender  bodies  having  a  blunt  nose,  the  inviscid  flow 
consists  of  an  inner  region  of  high  entropy  emanating  from  the  nearly  normal 
portion  of  the  detached  bow  shock  and  an  outer  region  of  low  entropy  flow 


16 


which  has  passed  through  the  oblique  portions  of  the  bow  shock.  Thus,  when 
the  boundary  layer  is  thin,  it  grows  within  the  high  entropy  layer.  Eventually 
all  of  the  entropy  layer  is  "swallowed'’  by  the  boundary  layer.  Then  the 
outer  edge  conditions  of  the  boundary  layer  are  governed  by  the  flow  passing 
through  a  much  weaker  bow  shock.  Clearly,  the  downstream  influence  of 
the  blunt-nose  complicates  the  problem  of  delineating  these  various  viscous 
regimes. 

Theoretical  studies  of  these  various  effects  in  the  flow  of  an  ideal 
gas  about  sharp  and  blunted  cones  and  "power-law"  bodies  have  been  reported 
in  Refs.  11  -  19.  The  more  recent  of  these  studies,  which  are  also  the  most 
complete  in  their  inclusion  of  the  effects  described  above,  require  numerical 
solutions  in  order  to  achieve  quantitative  results.  In  view  of  the  limited 
results  reported  and  the  lack  of  consideration  of  real-gas  effects,  only  qual¬ 
itative  estimates  can  be  made  of  the  flight  conditions  for  which  viscous  inter¬ 
actions  are  important.  Based  on  consideration  of  the  resuits  presented  an 
Refs.  11,  18,  and  19,  it  would  appear  that  for  the  smaller  value  of  R^.  (0.20  in) 
viscous  interaction  effects  are  negligible  below  an.  altitude  of  100  Kft,  become 
important  between  150  Kft  and  200  Kft  and  are  dominant  at  250  Kft. 

Finally,  we  must  consider  the  nature  of  the  air  chemistry  for 

the  flight  conditions  of  interest.  When  the  free-stream  static  temperature 

2  0 

is  low  (cr300*K),  the  equilibrium-air  solutions  for  cones  indicate 
that  at  values  of  M„  sin<T>  5  (where  (7*  is  the  shock  wave  angle),  the 
temperature  and  density  in  the  flow  behind  the  bow  shock  wave  differ 
by  10%  or  more  from  the  ideal-gas  values.  For  a  9*  half-angle  cone 
having  blunted  nose,  real-gas  effects  would  be  present  throughout  the 
flow  field  for  lvfm>  29.  At  M,*  =  25,  they  would  be  present  wherever 
(P  £  11,5*  or  over  most  of  the  flow  field  since  the  asymptotic  shock  angle 
is  just  10.  I*.  At  M &  ~  15,  real-gas  effects  would  be  present  for  <f  >  19.  5*. 

Since  the  detached  bow  shock  decays  to  this  value  within  three  nose  radii 
and  the  asymptotic  shock  angle  is  10.4*,  the  extent  of  real-gas  effects 
would  be  much  reduced  over  the  two  higher  Mach  number  cases  but  would 
extend  beyond  what  is  typically  termed  the  entropy  layer.  Thus  for  the 
flight  cases  of  interest  real-gas  effects  will  not  necessarily  be  restricted 
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to  the  region  of  the  blunt  nose  and  the  entropy  layer  emanating  from  this 
region.  On  the  other  hand,  whenever  the  wind  tunnel  must  resort  to 
partial  simulation,  the  high  free- stream  temperatures  may  introduce 
even  more  extensive  regions  of  real- gas  effects.  The  situation  in  this 
case  is  much  more  complex,  however,  since  the  high  free- stream 
temperatures  are  accompanied  by  lower  free-stream  Mach  numbers 
than  in  the  flight  case.  The  criterion  Afm sind'>  5  for  real-gas  effects  is 
not  appropriate  when  300* K.  Some  indication  of  the  presence  of 

real-gas  effects  in  the  wind  tunnel  case  when  none  existed  in  the  flight 
case  may  be  found  in  conical  shock  calculations  presented  in  Sec.  V.  C. 

The  flow  behind  the  detached  bow  shock  at  the  nose  must  be  treated 
as  a  real-gas  over  the  entire  range  of  flight  conditions  being  considered. 

The  question  then  becomes  one  of  determining  whether  the  flow  may  be 
treated  as  equilibrium  air  or  whether  finite-rate  chemical  processes  (i.e. 
nonequilibrium  phenomena)  are  important.  Criteria  for  delineating  regions 
of  frozen,  nonequilibrium  and  equilibrium  chemistry  are  rather  arbitrary. 

21  > 

Harney  and  Wittliff  et  alf  have  defined  different  criteria  which  are 

compared  in  Ref.  1  for  the  flow  about  a  blunt-nose  of  1  foot  radius.  At  the 

1,21 

boundary  between  frozen  and  nonequilibrium  chemistry  ’  only  bimo- 

22 

lecular  reactions  are  important.  Thus,  binary  scaling  can  be  used  to 
apply  Harney's  results  for  a  1  -ft  nose  radius  to  the  radii  of  present  interest. 

On  the  other  hand,  binary  scaling  does  not  govern  the  approach  to 
equilibrium  chemistry.  However,  the  nonequilibrium  normal  shock  wave 
solutions  reported  in  Ref.  1  have  been  used  to  define  a  distance  required  to 
reach  equilibrium.  If  this  distance  is  small  compared  t.i  the  shock  standoff 
distance  at  a  blunt  nose,  equilibrium  flow  may  be  assumed.  On  this  basis, 
equilibrium  flow  boundaries  have  been  estimated  for  the  nose  radii  being 
considered  in  the  present  study. 

Fig.  4  shows  the  boundaries  between  frozen,  nonequilibrium  and 
equiiil  rium  flow  for  R  0.  2"  and  1.  5"  as  determined  by  the  above  methods. 
For  the  flight  conditions  and  nose  radii  of  interest,  the  assumption  of  frozen 
flow  is  restricted  to  a  small  low-velocity,  high  altitude  region.  At  low 
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altitudes,  the  assumption  of  equilibrium  flow  is  valid  for  a  fairly  extensive 
range  of  velocities.  The  equilibrium  criteria  specified  in  Ref. .  1  (namely, 
that  the  temperature  be  within  1%  of  the  equilibrium  value)  is  quite  restric¬ 
tive,  however;  and  equilibrium  flow  probably  can  be  assumed  to  exist  over 
a  somewhat  broader  range  of  conditions  than  indicated  in  Fig.  4. 
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SECTION  IV 

REVIEW  OF  HYPERSONIC  SIMILITUDES 


A.  INVISCID  FLOW 

The  inviscid  flow  similitudes  that  are  appropriate  to  slender  bodies 

23 

are  based  on  hypersonic  small-disturbance  theory.  Tsien  showed  that  for 
the  hypersonic,  inviscid,  irrotational  flow  of  an  ideal  gas  with  constant 
specific  heats  past  affinely  related  slender  bodies,  the  requirement  of  dupli¬ 
cation  of  Mach  number,  per  se,  can  be  relaxed  provided  that  the  parameter 

_  24 

M^c  (where  t  is  a  characteristic  thickness  ratio)  is  kept  constant.  Hayes 

showed  the  equivalence  of  a  steady  hypersonic  flow  over  a  slender  body  with 

an  unsteady  flow  in  one  less  space  dimension  and  thus  extended  Tsien's 

result  to  rotational  flow. 

25 

Cheng  extended  this  classical  slender-body  hypersonic  similitude 
to  include  the  effects  of  small  nose  bluntness  and  equilibrium  real-gas 
phenomena.  Similitude  of  equilibrium  real-gas  effects  i?auires  duplication 
of  the  ambient  flow  thermodynamic  and  chemical  state  (i.  e.  pressure  p^  , 
density  ,  temperature  T*  ,  and  species  concentrations  ^  )  in  addition 
to  the  hypersonic  similarity  parameter  Mo, 't  .  Through  application  of  the 
blast-wave  analogy,  Cheng  found  that  for  axisymmetric  bodies  similitude  of  nose 
bluntness  effects  requires  duplication  of  a  biuntness  parameter 


K»  =  m!  C*  f 


(iv-l) 


where  Cp^  is  the  nose  drag  coefficient,  djq  is  the  diameter  of  the 
blunt  nose  and  L  is  the  body  length.  Cheng  states  his  similitude 
in  the  form 


|W(£ ^,/L;X/co)  (IV'2) 


with  similar  expressions  for  all  other  flow  variables. 

It  is  possible  to  rewrite  the  above  similitude  in  the  form 


20 


A  E  DC-T  R-69-36 


’‘•fi’tjS: ;  Ucotj  LL 


(rv-3) 


Eqs.  (IV-2)  and  (IV-3)  are  entirely  equivalent  when  the  free-stream  thermo¬ 
dynamic  state  is  completely  duplicated.  However,  as  mentioned  previously, 
mismatches  in  free-stream  temperature  are  unavoidable  at  many  altitudes 
because  of  wind-tunnel  perform?nce  limitations.  In  such  cases,  the  form 
of  the  similitude  given  in  Eq.  (IV-3)  is  to  be  preferred  over  that  given  in 
Eq.  (IV-2)  since  the  former  (by  virtue  of  the  Mach  number  independence 
principle)  represents  the  correct  similitude  when  M^^"  1,  even  when 

the  free-stream  temperature  is  not  duplicated. 

25 

Cheng  briefly  considered  the  question  of  extension  of  the  hyper¬ 
sonic  similitude  to  nonequilibrium  flows.  He  suggested  that  for  slender 
bodies,  a  nonsquilibrium  similitude  should  exist  for  a  two-component  react¬ 
ing  gas  when  »  1.  The  nonequilibrium  similitude  required  that  the 

particle  transit  time  L/U<»  be  duplicated  in  addition  to  ,  KN,  p ^  , p ^ 

and  X.  .  The  formal  development  of  the  nonequilibrium  hypersonic  simil- 

1  °°  26 
itude  for  a  single-species  diatomic  gas  has  been  carried  out  by  Inger  ,  who 

reached  the  same  conclusion  as  Cheng.  As  in  the  equilibrium  case  this 

similitude  contains  the  very  stringent  requirement  of  duplication  of  the 

free-stream  thermodynamic  state;  a  requirement  that  can  not  be  satisfied 

for  low-altitude  hypervelocity  flight  even  in  the  hypothetical  wind  tunnel 

described  in  Sec.  II.  In  Sec.  VI-D  an  extremely  useful  though  restrictive 

similitude  is  derived  in  which  the  requirement  of  the  duplication  of  the 

free-stream  thermodynamic  state  can  be  relaxed. 

It  should  be  pointed  out  that  Cheng's  nose  bluntness  parameter  does 
not  insure  similitude  of  the  flow  in  the  immediate  region  of  the  blunt  nose, 
but  only  provides  similitude  of  the  nose-bluntness  effects  in  the  flow  far 
downstream  of  the  nose. 
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The  Mach  number  independence  principle  for  the  flow  of  an  ideal 
gas  past  a  blunt  body  was  fir3£  stated  by  Oswatitsch^.  For  a  perfect  gas 
having  a  constant  specific -heat  ratio,  the  Mach  number  independence  principle 
specifies  that  the  flow  field  behind  the  bow  shock  wave  of  a  blunt  body  reaches 

a  limiting  configuration  as  the  free-stream  Mach  number  becomes  very  large. 

1 1 

Hayes  and  Probstein  have  generalized  this  principle.  They  show  that  two 
geometrically  similar  flows  of  a  given  gas  with  the  same  p^  and  U*  and 
large,  but  different,  are  not  merely  similar,  but  are  essentially  identical 

behind  the  bow  shock  wave.  Hence,  the  independence  principle  applies  to  all 
real  fluid  effects.  It  should  be  emphasized,  however,  that  whenever  dissipative 
processes  (i.  e. .  chemical  relaxation  or  viscous  phenomena)  are  present, 
the  bodies  must  be  not  merely  geometrically  similar  fc-.t  identical  since 
particle  transit  time  )  or  length  scale  (L)  must  be  matched  in  addition 

to  C  and  U-,  . 

s  w 

The  combined  requirements  of  the  Mach  number  independence  principle 
and  hypersonic  siender-body  similitude  dictate  that,  when  the  combined  effects 
of  nonequiiibrium  phenomena  and  nose  bluntness  are  present,  complete  dupli¬ 
cation  of  free-stream  thermo-chemical  state,  free-stream  velocity  and  body 
size  and  shape  is  necessary.  However,  for  those  cases  where  these  require¬ 
ments  cannot  be  fully  satisfied  because  of  wind-tunnel  performance  limita¬ 
tions,  it  is  perferrable  to  duplicate  p^  ,  U*  and  body  size. 

Clearly,  the  practical  usefulness  of  similarity  laws  depends  upon  their 
range  of  applicability.  A  similitude  cannot,  by  itself,  define  the  limits  for 
which  it  yields  reasonable  results.  This  must  be  accomplished  by  correlations 
of  theoretical  solutions  or  experimental  data  in  terms  of  the  similarity  para¬ 
meters.  Often,  the  similitude  parameters  are  applicable  over  a  wider  range 
of  Mach  numbers  and  slenderness  ratios  than  might  be  expected  from  the 
assumptions  made  in  their  derivations.  This  is  shown,  for  example,  in  the 
cot  relations  presented  in  Refs.  28  and  29  for  the  case  of  an  ideal  gas.  When 
real  gas  effects  are  present,  it  has  been  shown*’  that,  for  slender  bodies, 
the  onset  of  Mach  number  independence  does  not  occur  until  the  Mach  number 
is  significantly  higher  than  in  the  corresponding  ideal-gas  case. 
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B.  VISCOUS  FLOW 

The  similitudes  that  apply  to  hypersonic  flows  for  cases  where 
transport  effects  due  to  viscosity,  diffusion,  thermal  conduction  and 
the  like  play  a  role  have  been  reviewed  in  Ref.  1.  In  addition  Hayes  and 
Probstem^0  have  presented  a  fundamental  description  of  viscous  hypersonic 
similitude  for  slender  bodies  including  the  implications  of  the  Mach  number 
independence  principle  for  hypersonic  boundary  layers.  A  further  discussion 
of  the  combined  effects  of  nose  bluntness,  boundary- layer  displacement  and 
flow  nonequilibrium  on  slender-body  flow  fields  is  given  in  Sec.  VI-F  of  this 
report. 

30 

For  the  case  of  sharp-nosed  slender  bodies,  Hayes  and  Probstein 
note  that  "it  is  impossible  to  find  a  viscous  similitude  for  a  fluid  obeying 
a  general  equation  of  state.  "  They  go  on  to  point  out  that  when  the  inviscid 
flow  behaves  as  a  perfect  gas  and  only  the  boundary-layer  flow  exhibits 
real-gas  effects,  a  similitude  is  possible.  The  similitude  requirements 
for  affinely  related  bodies  are  that  the  following  quantities  be  duplicated: 

M-Z  ;  f  fa  ,  Z*L  ,  h„/H~ 

and  that  the  gas  be  the  same.  Whenever  /4>  ~  z~  Zi.  ,  the  enthalpy 
ratio  h^/H may  be  replaced  by  T  . 

In  Sec.  VI-F  the  viscous  similitudes  for  blunt-nosed  slender  bodies 
are  discussed  for  the  case  where  the  boundary  layer  is  thinner  than  the 
entropy  layer.  For  cases  of  interest  here,  the  Mach  number  independence 
principle  provides  similitude  of  the  flow  in  the  region  of  the  blunt  nose, 
even  when  viscous  effects  are  present,  if  ,  U*  and  the  nose  radius  are 
duplicated.  However,  the  requirement  that  both ^  and  be  duplicated  implies 
that  T"  is  duplicated  also.  Thus,  for  slender  bodies  having  blunt  noses  no 
similitude  exists  and  full  duplication  is  required.  An  approximate  similitude 
can  be  found  though  for  plane  flows  and  is  presented  in  Sec.  VI-F. 
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SECTION  V 

COMPARISON  OF  FLOWS  OVER  BLUNTED  SLENDER  CONES  FOR 
PARTIALLY  SIMULATED  CONDITIONS 

A.  INTRODUCTION 

We  now  turn  our  attention  to  the  flight  region  below  the  full  duplication 
boundary  of  the  hypothetical  wind  tunnel.  Fig.  2.  The  objective  of  the  study 
is  to  compare  flow  fields  in  flight  and  in  a  wind  tunnel  when  complete  dupli¬ 
cation  is  unattainable.  Consideration  is  restricted  to  a  blunted  slender  cone 

as  a  typical  body  shape  of  interest.  The  only  general  solutions  for  blunted 

31  32* 

cones  are  for  an  ideal  gas  and  involve  certain  approximations  '  .  Exact 

3  3  34 

numerical  solutions  have  been  reported  for  both  an  ideal  gas  ’  ’  and  equi- 
35 

librium  real  air  .  In  addition,  a  few  numerical  solutions  for  nonequilibrium 
real  air  have  been  calculated  (e.g.  Ref.  36).  These  existing  solutions  were 
insufficient  for  the  purposes  of  the  present  study  because  none  of  them  corre¬ 
spond  to  the  cases  of  partial  simulation  in  a  wind  tunnel  being  considered  here. 
Thus,  in  order  to  obtain  a  quantitative  comparison  of  flight  and  wind  tunnel 
flow  fields  over  a  slender  cone,  it  was  necessary  to  undertake  specific  solutions. 
As  indicated  in  Sec.  II  this  was  done  for  three  flight  conditions. 

In  the  present  study  a  cone  half-angle  of  9°  was  chosen  since  it  is  a 
shape  that  has  been  studied  extensively  (e.g.  Ref.  37).  In  order  to  inves¬ 
tigate  nose  bluntness  effects,  nose  radii  of  0.2  and  I.  5  inches  were  selected. 
Calculations  were  made  for  air  with  frozen,  finite-rate  and  equilibrium 
chemistry  to  assess  the  importance  of  real-gas  effects.  Only  the  body  stream¬ 
line  and  a  streamline  passing  through  the  outer  oblique  portion  of  the  bow  shock 
wave  were  considered.  These  represent  the  two  extremes  of  the  total  inviscid 
flow  field. 

The  free-stream  conditions  for  the  flight  and  wind  tunnel  cases  have 
been  compared  in  Sec.  II.  In  keeping  with  the  similitude  discussion  of  Sec.  IV, 
Um  and  9C  are  duplicated  in  the  present  calculations  rather  than  Mm9e  , 

In  this  section,  the  method  of  calculation  is  described  and  the  detailed  results 
are  presented  and  discussed. 


*  See  Sec.  VI  B  for  a  further  discussion  of  these  analyses. 
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B.  METHOD  OF  FLOW  FIELD  ANALYSIS 

1.  General  Approach.  The  flow  field  about  a  blunted  slender  cone  is 
characterized  by  a  detached  bow  shock  wave  at  the  nose  which  decays  to  a 
conical  shock  wave  far  downstream.  The  flow  entering  the  bow  shock  near  the 
axis  experiences  much  higher  entropy  and  temperature  levels  than  the  flow 
passing  through  the  more  oblique  outer  portion  of  the  shock.  Along  the  body, 
the  flow  is  initially  dominated  by  the  presence  of  the  blunt  nose.  However, 
with  increasing  distance  from  the  nose,  the  pressure  gradually  approaches 
the  equivalent  sharp  cone  value.  Rather  than  attempt  an  analysis  of  the 
complete  flow  field,  comparisons  of  the  flight  and  wind-tunnel  cases  have 
been  made  for  just  two  streamlines:  (1)  the  streamline  passing  through  the 
stagnation  point;  (2)  a  streamline  passing  through  the  conical  portion  of  the 
shock  wave. 

Because  the  body  streamline  will  not  differ  greatly  from  adjacent 
streamlines  entering  the  nearly  normal  portion  of  the  bow  shock,  the  com¬ 
parison  between  flight  and  wind  tunnel  for  the  body  streamline  will  also 

36 

typify  other  streamlines  in  the  entropy  layer.  In  addition,  Wood  et  al.  have 
shown  that  a  body  streamtube  calculation  based  on  the  pressure  distribution 
for  equilibrium  flow  is  in  good  agreement  with  an  exact  method  of  character¬ 
istics  solution.  The  calculation  for  the  conical-shock  streamline  typifies 
conditions  downstream  of  the  nose  where  the  flow  conditions  outside  the  entropy 
layer  approach  the  sharp  cone  values. 

The  body  streamtube  flow  has  been  calculated  using  the  CAL  computer 

program  for  the  quasi-one-dimensional,  inviscid  flow  of  reacting  gas  mix- 
38 

tures  .  One  of  the  capabilities  of  this  program  is  to  compute  the  flow  of  a 
general  mixture  of  reacting  gases  through  a  streamtube  having  a  specified 
pressure  distribution.  Thus,  a  first  step  in  the  present  study  was  the  deter¬ 
mination  of  the  pressure  distribution  on  a  blunted  9°  half-angle  cone  at  the 
flight  and  wind-tunnel  conditions  of  interest.  In  addition,  it  was  necessary 
to  specify  initial  conditions  at  the  stagnation  point.  For  the  latter,  thermo¬ 
dynamic  equilibrium  has  been  assumed.  This  assumption  is  well-justified 
for  the  two  velocities  at  80,000  ft  altitude  (Fig.  4);  however,  it  is  sowewhat 
questionable  for  the  case  at  180,  000  ft  particularly  for  the  smaller  nose 
radius.  Since  the  primary  objective  of  the  study  is  a  comparison  of  flight 
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and  wind-tunnel  flow  fields,  the  relative  differences  between  these  flows  are 
more  important  than  the  absolute  values.  On  this  basis,  the  error  introduced 
in  assuming  equilibrium  flow  at  180,  000  ft  altitude  is  believed  to  be  accept¬ 
able.  The  actual  equilibrium  stagnation-point  conditions  were  obtained  from 
an  equilibrium  normal  shock  calculation  followed  by  an  equilibrium,  isentropic 
compression  to  the  total  enthalpy  value  {i.  e.  zero  velocity). 

2.  Blunted  Cone  Pressure  Distributions.  The  pressure  distribution  was 
determined  from  a  study  of  existing  flow  field  solutions.  On  the  spherical  nose, 
the  pressure  was  found  to  be  insensitive  to  flow  chemistry.  This  is  illustrated 

in  Fig.  5  which  shows  pressure  distributions  on  a  spherical  nose  for  an  ideal 

.  33  35  39 

gas  (y  =  1.4),  real  air  in  equilibrium,  ’  and  nonequilibrium  real  air. 

The  nonequilibrium  pressure  distribution  was  obtained  using  the  computer 

program  reported  in  Ref.  40.  The  results  shown  in  Fig.  5  encompass  velocities, 

altitudes  and  nose  radii  appropriate  to  the  present  investigation.  Although  the 

ideal-gas  pressures  are  slightly  higher  than  the  equilibrium  air  values,  no 

appreciable  effect  of  chemistry  on  pressure  is  evident.  The  faired  curve  shown 

in  Fig.  5  has  been  used  in  the  present  study. 

Along  the  conical  afterbody  ideal -gas  and  equilibirum  air  pressure 
distributions  do  not  exhibit  the  same  agreement.  This  is  illustrated  in  Fig.  6 
where  theoretical  solutions  for  an  ideal  gas  (3^  =  1.4)  and  equilibrium  air 
are  correlated  as  Cp/z9^  versus  ( 9* /'J*C0u )(*■  /ZRn)  .  Distinct  dif¬ 

ferences  between  the  two  ideai-gas  solutions  and  the  several  equilibrium  air 
solutions  are  readily  apparent  in  this  figure.  Some  comment  on  this  correla¬ 
tion  is  appropriate  here.  The  exact  similitude  is  of  the  form 


(V- 1) 


If  equilibrium  real-gas  effects  are  to  be  included,  then  /i.  has  to  be  replaced 
byy2.  >  7Z  and  X;m.  Cheng  has  shown  that  in  the  limits  and 

7^-* /  ,  the  above  similitude  becomes: 


fw  _  r(  €k  X.  \ 

TljuM:#  t  ZR„) 


(V-2) 
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y/  __  / 

where  6  =  ypj~  ■  •  This  is  the  result  given  by  zeroth-order  thin  shock  laye  • 
theory.  Griffith  and  Lewis^  showed  that  using  Cp/Z9c  rather  than 

Mj 9c  yielded  better  correlation  of  experimental  data  obtained 

under  essentially  ideal-gas  conditions.  When  p* >>/A»  »  we  note  that  C/>/z&c 
reduces  to  Cheng's  parameter.  In  the  limit  /  both  parameters  pro¬ 

perly  express  the  dependence  of  pressure  on  M<»9C  •  For  values  of  /^,0C 
below  the  Mach  number  independent  regime,  Cp/2&p  appears  to  provide 
better  correlation  of  theoretical  solutions  and  experimental  data  even  though 
the  theoretically  proper  similitude  parameter  is  p# / li. p<*  9*.  •  The 

effect  of  a  finite  value  of  can  be  seen  in  Fig.  6  by  comparing  the  two 

ideal-gas  solutions  for  -  10  and  &=>  at  9C  =  10°.  The  major  differ¬ 

ence  is  in  the  region  of  the  minimum  pressure  where  both  the  magnitude  and 
the  location  are  different.  In  contrast,  the  equilibrium  air  solutions  show 
the  greatest  disagreement  close  to  the  spherical  nose  and  relatively  good 
agreement  near  the  pressure  minimum.  When  these  same  solutions  were 
correlated  in  terms  of  pw / Mp rather  than  Cp/2Qc  ,  greater 
scatter  was  found  without  any  discernable  trend  due  to  McoOc  .  The  solid 
curve  in  Fig.  6  drawn  through  the  equilibrium  air  correlation  was  adopted 
for  the  present  calculations. 

After  the  recompression  region  downstream  of  the  minimum  pressure 

point,  the  pressures  approach  the  equivalent  sharp  cone  pressures.  The 

sharp  cone  pressures  have  been  obtained  from  a  correlation  of  ideal-gas 
42  43  20 

solutions  '  and  equilibrium  air  solutions.  An  approximate  fairing  from 
the  distribution  given  in  Fig.  6  to  the  sharp-cone  pressure  for  each  free- 
stream  Mach  number  has  been  used.  The  final  result  is  shown  in  Fig.  7  as 
p / ftP  versus  .  The  curves  are  labeled  by  Case  number  (See  Sec.  II) 

and  whether  they  correspond  to  flight  (Fit.  )  or  wind  tunnel  (W.  T.  )  conditions 
as  well  as  by  free-stream  Mach  number.  The  pressure  distribution  for 
Case  III  W.  T.  -  27.  34  has  been  omitted  to  avoid  cluttering  the  figure 

since  it  lies  between  the  two  lower  curves..  It  should  be  emphasized  that 
these  curves  do  not  exactly  account  for  the  effect  of  different  values  of 
but  have  inherent  in  them  the  scatter  observed  in  Fig.  6. 

The  final  step  in  specifying  the  pressure  distribution  for  the  computer 
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program  was  to  fit  polynomial  functions  to  the  faired  curves  shown  in  Figs.  5 
and  7  for  each  value  of  Each  pressure  distribution  was  divided  into  8 

intervals  so  that  only  third-order  polynomials  were  needed. 

3.  Streamtube  Computer  Program.  The  computer  program  utilized  for 

the  body  streamtube  calculations  provides  numerical  solutions  for  the  quasi- 

one -dimensional,  inviscid  flow  of  chemically  reacting  gas  mixtures  as  well 

38 

as  the  limiting  cases  of  frozen  and  equilibrium  flow  .  Vibrational  excitation 
can  be  either  frozen  at  the  initial  level  or  in  equilibrium  with  the  translational 
temperature.  The  latter  option  was  utilized  in  the  present  work.  Electronic 
excitation  is  assumed  to  be  in  equilibrium  with  the  translational  temperature 
also.  For  the  flight  cases,  calculations  were  made  for  frozen,  finite-rate, 
and  equilibrium  chemistry.  Whereas  for  the  wind-tunnel  cases,  only  finite- 
rate  or  nonequilibrium  chemistry  calculations  were  performed. 

The  calculations  for  Case  I  ((/«*=  15,000  fps)  were  made  with  an  air 

model  consisting  of  8  species  (N->,  07,  A  ,  N,  O,  NO,  NOT  and  e  )  and  11 

3  8  c  c  r 

chemical  reactions  .  This  simple  air  model  is  justified  because  the  stag¬ 
nation-point  temperature  was  less  than  6000  °K,  hence  NO+was  the  only 
significant  ionized  species.  For  Case  II  (  (J^  -  25,  000  fps),  the  stagnation- 
point  temperatures  approach  9000  "K,  and  it  is  necessary  to  include  other 
ionized  species.  Consequently,  a  12-species,  39-reaction  air  model  was 
utilized.  The  additional  ionized  species  wc  .e  »  N  and  O  .  The 

12-39  air  model  was  also  used  for  Case  III  (  Uco  -  35,  000  fps)  where  the 
stagnation-point  temperature  is  somewhat  higher  than  11,000  °K. 

4.  Conical-Shock  Calculations.  Conical-flow  solutions  have  been  correlated 

in  terms  of  the  hypersonic  similarity  parameter  M„sin  0  for  an  ideal 

44  20  C 

gas  and  for  equilibrium  real  air.  In  the  latter  case,  correlation  obtains 

25 

only  at  a  fixed  free- stream  state  (i.  e. ,  at  a  fixed  altitude);  this  is  equiv¬ 
alent  to  saying  that  the  solutions  for  equilibrium  air  can  be  correlated  in 
terms  of  sin  9^  at  each  altitude  (see  Sec.  II).  In  the  present  study, 

LL  and  &  are  duplicated,  while  M  *  sin  0  is  not. 

v->  C 

A  conical-flow  streamline  experiences  a  discontinuous  compression 
at  the  shock  wave  followed  by  an  additional  isentropic  compression  downstream 
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of  the  shock.  For  siender  cones  at  high  Mach  number,  the  isentropic  compres¬ 
sion  is  small  compared  to  the  shock  compression.  For  the  range  of  conditions 
under  consideration  here,  the  pressure  on  the  cone  surface  is  from  6  percent 
to  13  percent  higher  than  the  pressure  just  behind  the  conical  shock,  the 
density  is  from  4  percent  to  9  percen:  greater,  and  the  temperature  from 
1.  5  percent  to  3.  5  percent  greater.  Thus,  for  present  purposes  of  compar¬ 
ison  of  flight  and  wind-tunnel  situations,  it  is  adequate  to  calculate  only  the 
conditions  immediately  behind  the  shock.  In  these  calculations,  account  must 
be  taken  of  the  high  free-stream  static  temperatures  occurring  in  the  wind- 
tunnel  cases  because  they  introduce  real-gas  effects  at  lower  values  of 
Mwsin  0^  than  predicted  by  the  solutions  reported  in  Ref.  20.  Utilizing 
the  fact  that  the  shock  angle  and  pressure  are  relatively  insensitive  to  real- 
gas  effects,  it  is  possible  to  use  Refs.  20,  42,  43  to  obtain  solutions  for 
equilibrium  air  for  the  wind-tunnel  cases.  The  highest  temperature  computed 
for  these  cases  was  less  than  2300  °K.  At  such  temperatures,  any  nonequi¬ 
librium  phenomena  will  not  introduce  variations  in  the  flow  quantities  that 
are  any  larger  than  the  already  neglected  effects  of  the  isentropic  compression. 
Hence,  it  is  sufficient  to  consider  just  the  equilibrium-air  case. 

For  these  conical- shock  calculations,  any  nozzle  flow  nonequilibrium 
effects  will  influence  the  flow  downstream  of  the  conical  shock  wave.  To 
assess  this  influence,  conical-shock  calculations  were  also  made  for  frozen 
flow  utilizing  the  estimated  nozzle  nonequilibrium  conditions  presented 
earlier.  The  results  are  then  compared  with  the  equilibrium  air  calcula¬ 
tions. 

C.  RESULTS  OF  CALCULATIONS 

1.  Equilibrium  Stagnation- Point  Conditions.  Before  presenting  the  results 
of  the  streamtube  calculations,  the  initial  stagnation-point  conditions  will  be 
compared.  For  the  flight  cases,  these  results  were  obtained  from  Ref.  10 
whereas  the  wind-tunnel  results  were  calculated  from  an  equilibrium  normal- 
shock  solution  followed  by  an  isentropic  compression  to  the  total  enthalpy 
state  (i.  e.  ,  zero  velocity).  Since  the  nozzle-flow  nonequilibrium  effects  are 
not  large  for  these  three  cases  {Table  II-4),  the  wind-tunnel  equilibrium  free- 
stream  conditions  have  been  used  {Tables  II- 1,  -2,  -3).  The  following  tables 
present  the  pressure,  temperature,  density  and  species  concentrations  (moles/ 
of  mixture)  at  the  stagnation  point  for  thermochemical  equilibrium. 
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Table  I- 1 

COMPARISON  OF  STAGNATION  -  POINT  CONDITIONS 
Case  I,  15,000  fps 


THERMODYNAMIC  QUANTITY 

FLIGHT 

WIND  TUNNEL 

PRESSURE,  ^'(atm) 

8.701 

8.727 

TEMPERATURE,  7^'{*K) 

5172 

5650 

DENSITY,  fl'  (mis) 

0.3632 

0.3500 

SPECIES  (MOLES/gm  OF  MIXTURE) 

"2 

2.573xl0‘2 

2.56x10-2 

°2 

2.291xl0-4 

1.663xl0-4 

H 

1. 002xl0-3 

1.109x10-3 

0 

1.256xl0-2 

1.283x10*2 

NO 

1.167xl0-3 

1.323.(10-3 

N0+,  e 

2.123x10“® 

3. 269xlO-6 

Table  1-2 

COMPARISON  OF  STAGNATION-POINT  CONDITIONS 
Case  33,  U„=  25,000  fps 


THERMODYNAMIC  QUANTITY 

FLIGHT 

WIND  TUNNEL 

PRESSURE,  p'0  (atm) 

21.13 

21.50 

TEMPERATURE,  Te'  (°K) 

8632 

8730 

DENSITY,^'  (ama) 

0.1811 

0.1701 

SPECIES  (MOLES/gm  OF  MIXTURE) 

n2 

1.289x10"2 

1.209x10*2 

02 

9.153x10-® 

8.355x10-3 

N 

2.767x10-2 

2.928x10-2 

0 

l.lOlxlO"2 

l.lOlxlO"2 

NO 

1.059x  10"4 

3. 801x 10"4 

e" 

7.88x10"5 

8. 710x1  O'5 

H2+ 

1.939x10“6 

5.580x10-3 

°2+ 

8.981. 10’8 

9.255xl0-8 

N+ 

2.122xlO'S 

3.001xl0-5 

0+ 

9.088x)0-6 

I.055xl0-5 

N0+ 

1.051x IQ"8 

1. 115x10-3 
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Table  1-3 

COMPARISON  OF  STAGNATION-POINT  CONDITIONS 
Care  m,  =  35,000  fps 


THERMODYNAMIC  QUANTITY 

FLIGHT 

WIND  TUNNEL 

PRESSURE,  p0  (atm) 

0.6202 

0.6199 

TEMPERATURE,  T0  (°K) 

11,190 

1 1 , 230 

DENSITY,  y3/( ama) 

0.0703 

0.0698 

SPECIES  (MOLES/gm  OF  MIXTURE) 

h2 

2.642x!0-8 

2. 500x1 O"5 

°2 

3. 1 36x10“® 

3.047x1 0"8 

N 

•T.  901  x  1 0-2 

“1.885x1  O'2 

0 

l .  365x10" 2 

1. 363x1 0"2 

NO 

1.365xl0"6 

1.313xl0"6 

e“ 

5. 688x10" 3 

5.883x1 0~3 

N2+ 

2.044xl0"6 

2.024x1 0"6 

o 

K> 

4- 

1.704xl0"8 

1 .703x1 0“8 

N+ 

4.854xl0'3 

5.020x1 O"3 

0+ 

S.293xl0"4 

8. 573x1 0"4 

N0+ 

3.369x10-6 

3.297xl0"6 
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In  comparing  the  flight  and  wind-tunnel  conditions  it  is  seen  that  all  of  the 
thermodynamic  state  variables  [fa' ,  T0' t  and  p»  '  differ  by  less  than  4 
percent  for  the  three  cases.  This  is  a  direct  result  of  the  flow  being  Mach 
number  independent.  On  the  other  hand,  mismatches  of  up  to  40  percent 
occur  for  the  species  concentrations  but  only  in  the  case  of  species  having 
very  small  concentrations.  These  large  mismatches  result  from  being  in  a 
temperature -density  region  where  certain  species  concentrations  are  very 
sensitive  to  temperature  (e.g.  ,  see  Ref.  45). 

2.  Body  Streamtube,  Case  I  (  (J^  =  15,000  fps).  Typical  results  of 
the  body  streamtube  calculations  for  Case  I  are  shown  in  Figs.  8-11.  In 
each  figure,  the  four  curves  represent  the  frozen,  equilibrium  and  finite - 
rate  (or  nonequilibrium)  solutions  for  the  flight  case.  The  symbols  in 
these  figures  represent  the  nonequilibrium  solutions  for  the  wind-tunnel 
case.  In  these  figures  the  results  are  shown  as  a  function  of  distance  along 
the  cone  surface  in  units  of  nose  radii  (S/R^).  When  nondimensionalized 
in  this  manner,  note  that  the  frozen  and  equilibrium  solutions  are  independent 
of  body  nose  radius  while  the  nonequilibrium  solutions  are  not. 

The  temperature  (T/T^  )  is  shown  as  a  function  of  S/R^  in  Fig.  8. 
Presenting  the  temperature  as  T/T^  masks  the  effects  of  the  small  differ¬ 
ences  in  Tq  between  flight  and  wind  tunnel,  and  emphasizes  the  differences 
resulting  from  the  mismatch  in  M^  .  Comparing  the  flight  and  wind-tunnel 
nonequilibrium  solutions,  it  is  noted  that  the  temperature  ratios  T/T^ 
differ  by  less  than  10  per  cent  for  either  nose  bluntness  even  though  the  free- 
stream  Mach  numbers  differ  by  over  40  percent  and  the  pressures  along  the 
cone  differ  by  up  to  25  percent  (Fig.  7).  The  importance  of  flow  chemistry 
and  its  dependence  on  nose  bluntness  are  clearly  evident  in  Fig.  8  when  the 
flight  solutions  for  frozen,  finite -rate  and  equilibrium  flows  are  compared. 
Note  that  along  the  cone  the  equilibrium  temperature  is  about  60  percent 
higher  than  the  frozen-flow  temperature.  Although  equilibrium  and  frozen 
solutions  are  not  pn  sented  for  the  wind-tunnel  case,  they  lie  about  10  per¬ 
cent  above  the  corresponding  flight  solutions. 

The  density  distributions  (  p / p^  )  are  shown  in  Fig.  9.  In  general, 
the  density  follows  the  same  pattern  as  the  temperature  except  that  the 
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relative  levels  of  the  frozen  and  equilibrium  solutions  are  inverted.  Again, 
the  wind-tunnel  densities  for  nonequilibrium  flow  are  within  iO  percent  of 
the  corresponding  flight  value?;. 

Along  the  spherical  nose  (S/RN  -  1.4),  where  the  flow  is  Macn 
number  indc  per  dent,  the  wind-tunnel  temperatures  and  densities  are  essen¬ 
tially  identical  to  the  flight  values. 

Although  the  velocity  distributions  are  not  shown,  they  increase 
rapidly  along  the  spherical  nose  and  lie  between  60  percent  and  70  percent 
of  U0 0  along  the  conical  afterbody  for  all  cases.  Typically,  the  wind- 
tunnel  values  lie  about  1  percent  to  2  percent  above  the  flight  case  for  non¬ 
equilibrium  flow.  The  frozen-flow  velocities  are  about  7  percent  lcv/er 
than  the  equilibrium  velocities. 

In  examining  the  distributions  of  species  concentrations,  it  is  inter¬ 
esting  to  note  that  nonequilibrium  distributions  are  not  necessarily  bounded 
by  the  frozen  and  equilibrium  solutions  as  are  the  velocities  and  thermo¬ 
dynamic  state  variable  (T  andfl).  This  is  the  case,  for  example,  for  nitric 
oxide  (NO)  which  is  shown  in  Fig.  10.  The  nonequilibrium  concentrations 
fall  below  the  frozen  level  while  the  equilibrium  concentrations  lie  above 
the  frozen-flow  value.  A  comparison  of  the  flight  and  wind-tunnel  nonequili¬ 
brium  concentrations  indicates  that  the  initial  differences  are  approximately 
maintained  up  to  the  point  where  the  concentrations  start  to  approach  equili¬ 
brium. 

The  atomic  oxygen  (O)  concentrations  shown  in  Fig.  11  are  an  example 
where  the  nonequilibrium  case  is  bounded  by  the  frozen  and  equilibrium 
distributions.  The  effect  of  nose  bluntness  is  quite  striking  here.  For 
=  0.2  in.  ,  the  oxygen  atom  concentrations  essentially  remain  frozen 
at  their  initial  value  well  downstream  of  the  nose,  while  for  =1.5  in. 
the  concentrations  show  a  steady  decrease  toward  the  equilibrium  case. 

At  Un  =  15,000  fps,  molecular  nitrogen  is  relatively  undissociated. 
The  total  variation  from  the  stagnation-point  concentrations  is  only  4  percent 
for  all  cases.  The  other  species  (N,  O-,,  NO+  and  e  )  do  not  exhibit  any 
unique  behavior  and  are  not  presented  either. 
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3.  Coaical-Sb&ck  Calcalatioss,  Case  I  {  U»-  15,090  fps).-  The  relatively 

small  variation  in  pressure,  density  and  temperature  between  the  conical 
shock  wave  and  ccne  surface  has  been  cited  previously  as  the  justification  tor 
comparing  only  the  conditions  just  behind  the  conical  portion  of  the  shock  wave 
for  the  wind  tunnel  and  flight  cases.  This  comparison  is  shown  in  the  following 
table.  For  convenience,  a  comparison  of  seme  of  the  free-stream  conditions 
are  tabulated  also.  The  subscript  "s"  refers  to  equilibrium  conditions  just 
behind  the  shock  wave.  For  Case  1,  the  effect  of  nozule-flow  nonequlilbriam 
on  the  thermodynamic  state  variables  is  less  than  1  perce*%  and  has  not  been 
considered  here. 

Table  T-? 

CCHPiSlSSS  OF  F2EE-ST5EJ&  IS 3 
COEICil-SdSCI  CGS31T1CSS 
Case  I  (£»„  =  15. COD  fps) 


?JW“J*ET£5 

iRLISHif 

'  ¥«*s  iwxa  . 

«- 

15-3« 

9.0S2 

fm  (***) 

2.7S~zV3~2 

S.UtolO-2 

T.{*«) 

220-9 

£50.0 

-  ] 

I*?*) 

1S.765 

1«.7*9  1 

ft 

0- 1265 

0.QSS5 

fs 

0.2519 

0.3005 

t5  (*«) 

5W.0 

959.0  i 

The  above  table  indicates  aooreciable  differences  exist  between  the 
flight  and  wind  tunnel  flews  downstream  of  the  conical  shock.  Clearly,  this 
flow  s«  not  Mach  number  independent  as  was  the  entrooy- layer  flow  typified 
by  the  body  streamtube  calculation.  In  the  present  case,  if  Atf^sin  6c 
rather  than  {/sin  S  were  matched  between  flight  and  wind  tunnel,  the 

GO  £ 

shock  ratios  such  as  p®  •  Pspo  »  ar-°  %  ‘  Teo  would  be  closely 

matched  also.  However,  then  the  temperature  behind  the  shock  T_  in  the 
wind  tunnel  would  be  about  1500  CK  and  real-gas  effects  would  begin  to  be 
felt. 


34 


AEDC-TR-49-26 


4.  Body  Streamtube,  Case  II  (U  g>  =  25,  000  fps).  Typical  results  of  the 
body  streamtube  calculations  for  Case  II  are  shown  in  Figs.  12  -  1 6.  The 
figures  follow  the  same  format  as  the  set  for  Case  I;  namely,  the  four  curves 
are  for  the  flight  cases  and  the  two  sets  of  symbols  represent  the  nonequiii- 
brium  wind-tunnel  cases  (R^.  =  0.  2  in.  and  1. 5  in. ).  Again  the  abscissa  in 
each  figure  is  distance  along  the  cone  in  units  of  nose  radii  (S/R^.). 

The  temperature  distributions  (T/T^  )  are  shown  in  Fig.  12.  The 
general  features  are  the  same  as  Case  I  (Fig.  8);  the  flight  solutions  orient 
themselves  in  the  order  equilibrium  flow’,  nonecui  librium  flow  for  the  blunter 
nose  (R^.  =  1.5  in. ),  nonecuiiibrium  flow  for  the  sharper  nose  (R^  =  0.  2  in.  ), 
and  finally  frozen  flow.  In  comparing  Case  I  and  Case  II,  it  is  noted  that  for 
Case  II  the  nonequilibrium  solutions  lie  closer  to  the  equilibrium  solution 
thar.  in  Case  I  and  actually  attain  the  equilibrium  levels,  also  the  frozen  and 
equilibrium  solutions  lie  farther  apart  than  in  Case  L  Ko  simple  explanation 
can  be  given  for  the  closer  agreement  of  the  equilibrium  and  nonequilibrium 
solutions  for  Case  EL  It  results  from  complex  differences  in  the  air  chemistry 
for  the  two  cases.  Some  of  these  differences  will  be  noted  in  the  discussion 
of  the  species  concentrations.  The  larger  differences  between  the  frozen 
and  equilibrium  solutions  in  Case  II  compared  to  Case  I  occur  because  the 


higher  stagnation -point  temperature  in  Case  «i  produces  a  greater  degree 
of  dissociation  which  is  maintained  in  the  frozen  flow  calculation. 


Comparing  the  flight  and  wind-tunnel  nonecuiiibrium  solutions  in 
Fig.  12,  it  is  observed  that  the  wind-tunnel  temperatures  are  at  most  about 
S  percent  higher  than  the  flight  values.  Referring  to  Fig.  7,  it  is  noted  that 
the  pressures  in  the  wind-tunnel  case  are  23  percent  higher  than  the  flight 
case  and  the  free -stream  Mach  numbers  differ  by  a  factor  of  iwc. 

The  density  distributions  ( p  p')  for  Case  II  are  shown  in  Fig.  13. 
The  behavior  is  analogous  to  that  of  the  temperature.  The  nonecuiiibrium 
solutions  do  not  depart  from  equilibrium  as  much  as  in  Case  I  while  the 
frozen  and  equilibrium  solutions  are  appreciably  farther  apart.  Along  the 


cone,  the  density  levels  in  the  wind-tunnel  case  are  about  1C  percent  to  20 
percent  higher  than  in  the  flight  case. 
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Typical  species  concentrations  are  shown  in  Figs.  14  -  16.  Molecular 
nitrogen  (N^)  is  shown  in  Fig.  14.  It  i3  seen  that  in  Case  II,  the  nonequilibrium 
concentrations  lie  between  the  frozen  and  equilibrium  levels.  Also,  the  wind 
tunnel  concentrations  uniformly  lie  below  the  corresponding  flight  values  by  an 
amount  essentially  equal  to  the  initial  differences.  Finally,  the  variation  of 
N?  concentrations  along  the  cone  is  much  greater  than  in  Case  I  (60  percent 
versus  4  percent). 

The  nitric  oxide  (NO)  concentrations,  shewn  in  Fig.  15,  exhibit  a 
somewhat  different  behavior  from  Case  I  (Fig.  10).  The  first  notable  differ¬ 
ence  is  in  the  equilibrium  case.  In  Case  I,  equilibrium  chemistry  resulted  in 
a  net  production  cf  NO  above  the  initial  concentration  along  the  cone,  whereas 
in  Case  II  NO  is  consumed  and  falls  below  the  initial  concentration.  The  non- 
equiiitrium  solutions  show  an  even  greater  reduction  in  NO  concentration 
followed  by  a  return  to  the  equilibrium  levels.  For  the  blunter  nose,  the 
equilibrium  chemistry  produces  fluctuations  in  the  NO  concentrations  just 
downstream  of  the  nose.  These  are  observed  in  both  the  flight  and  wind- 
tunnel  solutions.  The  agreement  between  the  flight  and  wind-tunnel  solutions 
is  quite  good  up  to  about  5/ =  2.  After  that,  larger  differences  are  noted 
until  the  final  approach  to  equilibrium. 

For  Case  II,  oxygen  is  essentially  fully  dissociated  and  as  a  result 

the  concentrations  vary  by  less  than  4  oercent  for  all  cases  throughout  the 

_  _?  _? 
flow  lield.  Tne  concentrations  lie  m  the  range  1.40  x  10  to  1.45  x  10 

moles-  gm  of  mixture  and  are  not  shown. 

The  atomic  nitrogen  (N)  concentrations  are  shown  in  Fig.  16.  Like 
molecular  nitrogen,  the  nonequilibrium  concentrations  lie  between  the  frozen 
anc  equilibrium  cases,  and  approach  the  equilibrium  levels  at  large  values 
of  S'R^..  Also,  the  initial  differences  in  concentration  between  the  wind- 
tunnel  and  flight  cases  tend  to  persist  throughout  the  flow.  For  brevity,  the 
distributions  of  the  other  species  are  omitted. 

5.  Conical-Shock  Calculations,  Case  G  (Uco  =  25,000  fps).  For  Case  n, 
the  differences  in  7^,,  and  tAco  between  wind  tunnel  and  flight  are  approx¬ 

imately  twice  as  great  as  in  Case  I.  in  addition,  the  effects  of  nozzle-flow 
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nonequilibrium,  while  small  compared  to  the  mismatches  in  free- stream 
conditions,  are  not  negligible  as  they  we:e  for  Case  I.  In  taking  nozzle-flow 
nonequilibrium  effects  into  account,  two  limiting  cases  exist:  either  the  flow 
equilibrates  immediately  behind  the  conical  shock  or  the  composition  remains 
frozen  at  the  free-stream  values.  In  both  cases,  the  reduction  in  free-stream 
pressure  and  temperature  due  to  nozzle  nonequilibrium  is  more  significant 
than  the  increase  in  M  m  .  As  a  result,  the  effect  of  nozzle-flow  nonequili¬ 
brium  is  to  improve  the  agreement  between  wind  tunnel  and  flight.  The 
following  table  compares  free-stream  and  conical-shock  conditions  for  the 
case  of  equilibrium  nozzle  flow  and,  hence,  represents  the  maximum  mis¬ 
matches  that  could  occur. 


Table  1-5 


COMPARISON  OF  FREE-STREAH  AND 
CONICAL-SHOCK  CONDITIONS 
Case  II  [Um=  25.000  fps) 


j  PACKET'S 

FIS6HT 

¥1X0  TUSSEL  j 

i 

•  Mm  !  25.57 

(ais)  j  2.765x10'^ 

T,  (’*}  |  220.9 

;  Us  {*?*)  75.625 

!  pi  (ana)  0.1737 

j>s  (at n)  0.6536 

i  Ts  (•*)  |  >032 

n.25 

l.513x!0-1  ! 

,209 

25.53C 

0.H39 

0.7907 

,957 

As  ir.  Case  I,  the  differences  between  flight  and  wind-tunnel  conditions 
are  appreciable,  although  less  than  the  free-stream  differences.  Again,  if 

sin  8c  rather  than  sin  8c  were  duplicated  in  the  wind  tunnel,  the 

shock  ratios  (psJ  p^  etc.)  would  be  much  closer.  However,  the  temperature 
T  in  the  wind  tunnel  case  would  then  be  about  4000  °K  and  real-gas  effects 
would  be  important. 
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6.  Body  Streamtube,  Case  III  (Ua  =  35,  000  fps).  Before  presenting  the 
calculated  results  for  this  case,  a  strong  note  of  caution  in  accepting  the 
absolute  values  as  accurate  must  be  made.  At  180,  000  ft  altitude,  the  relax¬ 
ation  zone  behind  a  normal  shock  wave  at  Uoo  =  35,  000  fps  is  on  the  order 
of  1/2  inch  thick*  which  is  larger  than  the  shock  stand-off  distance  for  both 
nose  radii  of  interest.  Furthermore,  Fig.  3  indicates  that  low-density  viscous 
effects  will  be  present  in  the  flow  about  the  blunt  nose.  Also  viscous  interaction 
effects  will  exist  in  the  flow  along  the  conical  afterbody.  Thus,  an  accurate 
prediction  of  such  a  complex  flow  must  include  the  combined  effects  of  thermo¬ 
chemical  nonequilibrium  and  viscous  interactions  both  in  the  region  of  the 
blunt  nose  and  along  the  conical  afterbody.  Undertaking  such  a  solution  was 
not  felt  justifiable  in  view  of  the  relatively  small  differences  between  the 
flight  and  wind-tunnel  free- stream  flows  at  these  velocity  and  altitude  con¬ 
ditions  (Table  II-3).  In  fact,  when  the  wind-tunnel  free-stream  conditions 
are  corrected  for  the  estimated  nozzle  nonequilibrium  effects  (Table  II-4), 
the  flight  and  wind-tunnel  conditions  are  nearly  identical.  This  is  shown  in 
the  following  table  which  compares  flight  conditions  with  equilibrium  and 
estimated  nonequiiibrium  wind-tunnel  conditions. 


Table  1-6 

CGHPARIS0X  OF  FLIGHT  AMD  HIKD-TUKXEL  COXDITIONS 
Case  IIX  f/j  =  35,000  fps  180,000  ft  Ait. 


PARAMETER 

FLI6KT 

EQUILIBRIUM 
KIKD  TUKREL 

K0XEQU1LI3KIUH 
VINO  TUXHEL 

VELOCITY.  i/0  (fps) 
DEXSiTY,  (tua) 

PRESSURE,  jp  Ut») 

rco 

TEMPERATURE,  Ta  (*K) 
HACH  KUK5ER,  M 

CO 

35,000 

5.S13X10"’ 

?. 292x10"" 

265.9 

32.63 

35,000 
?.  413x10"^ 
6.135x10"^ 

380.0 

27.3? 

3?, COO 

?.<WxlG"? 

SJ.295X10"’ 

255 

33.1 

Thus  for  Case  III  the  primary  difference  between  flight  conditions  and 
the  wind  tunnel,  when  nonequilibrium  effects  are  taken  into  account,  is  in  the 
composition  of  the  test  gas  (Table  II-4). 
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In  order  to  show  the  maximum  differences  that  could  be  expected  for 
this  case,  the  equilibrium  wind-tunnel  conditions  have  been  used  rather  than 
the  estimated  nonequilibrium  conditions.  The  body  streamtube  calculations 
were  made  in  a  manner  completely  analogous  to  the  two  preceding  cases.  The 
equilibrium  stagnation -point  conditions  have  been  compared  in  Table  V-3  and 
indicate  very  little  difference  between  flight  and  wind  tunnel.  Thus,  while  the 
calculated  results  may  not  be  quantitatively  accurate  the  relative  differences 
are  probably  representative  of  the  actual  flow  situation.  Typical  results  for 
Case  III  are  presented  in  Figs.  17  to  20.  As  before  frozen,  nonequilibrium 
and  equilibrium  solutions  are  shown  for  the  flight  case  (solid  curves)  and 
nonequilibrium  solutions  only  for  the  wind-tunnel  (symbols).  The  abscissa 
in  each  figure  again  is  distance  along  the  cone  in  units  of  nose  radii. 

The  temperature  distributions  (T/Tj  )  shown  in  Fig.  17  are  very  much 
like  those  for  Case  I  except  that  the  difference  between  the  frozen  and  equi¬ 
librium  solutions  is  more  pronounced.  However,  the  wind-tunnel  solutions 
are  essentially  identical  to  the  flight  solutions.  This  is  to  be  expected  in  view 
o:  the  small  differences  in  tree-stream  conditions  between  wind  tunnel  and 
flight.  Also,  we  note  that  the  effect  of  nose  radius  on  the  nonequilibrium  flow 
is  much  smaller  for  Case  III  than  for  Case  II  (Fig.  12). 

The  density  distributions  (  p/p/  )  for  Case  HI  are  shown  in  Fig.  18. 
The  behavior  is  very  similar  to  Case  H  (Fig.  13)  except  that  the  wind-tunnel 
and  flight  solutions  are  in  much  better  agreement  than  in  Case  H.  Again,  the 
very  small  effect  of  nose  radius  on  the  nonequilibrium  flow  is  noted. 

In  the  assumed  equilibrium  stagnation  region  for  Case  HI,  nitrogen 

and  oxygen  are  essentially  completely  dissociated.  Thus,  atomic  nitrogen  is 

the  dominant  species  and  atomic  oxygen  is  the  second  most  abundant.  Little 

recombination  occurs  in  the  expansion  along  the  body  so  that  the  concentrations 

of  these  species  do  not  vary  greatly.  The  atomic  nitrogen  (K)  concentrations 

are  shown  in  Fig.  19,  where  it  is  noted  that  the  total  variation  is  less  than 

10%  of  the  stagnation-point  concentration.  For  equilibrium  and  nonequilibrium 

flow  the  variation  of  the  concentrations  with  S/R„  is  auite  different  in  Case  HI 

i\' 

than  in  Case  H  (Fig.  16).  In  Case  II,  both  equilibrium  and  nonequilibrium 
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solutions  showed  a  decrease  in  concentration  along  the  body.  In  Case  III,  the 
nonequilibrium  solutions  show  a  general  increase  in  atomic  nitrogen  con¬ 
centration  with  increasing  S/R^  ,  while  the  equilibrium  solution  yields  a 
very  nonmonotonic  behavior.  Because  the  atomic  oxygen  concentrations  only 
vary  +5%  from  the  stagnation-point  concentration,  they  are  not  shown. 

The  atomic  nitrogen  is  about  10%  ionized  at  the  conditions  for  Case  III 
and  is  the  largest  source  of  electrons.  The  atomic  nitrogen  ion  (N  )  concen¬ 
trations  are  shown  in  Fig.  20.  Within  a  distance  of  100  nose  radii  the  non¬ 
equilibrium  concentrations  decrease  by  a  factor  of  3  to  6  (depending  on  R^) 
below  the  initial  values,  whereas  the  equilibrium  solution  shows  a  decrease 
of  well  over  an  order  of  magnitude.  The  electron  concentrations  show  a 
similar  behavior. 

It  can  be  seen  in  Table  V-3  that  the  other  species  are  much  less 
abundant;-  hence,  their  concentrations  are  not  shown. 

7.  Conical-Shock  Calculations,  Case  III  (U^  =  35,000  fps).  Using  the 
equilibrium  wind-tunnel  conditions,  which  provide  the  largest  differences 
from  the  flight  conditions,  the  equilibrium  conditions  behind  a  conical  shock 
wave  nave  been  calculated  and  are  shown  in  the  following  table.  The  differ¬ 
ences  are  seen  to  be  very  small. 


Table  Y-7 

COMPARISON  OF  COX ICAL-SHGCK  CONDITIONS 
Case  LU  (  £4=  35,000  fps) 


PARAMETER 

FLiSfiT 

|  X I M0  TUMMEL 

us 

35,500 

!  35,500 

* 

Ps  (“a) 

G. 0025s 

!  0.00251 

j>s  (at*) 

0.0161 

\  0.0167 

Ts  (*K) 

1710 

j  1650 
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D.  DRAG,  SKIN- FRICTION  AND  HEAT  TRANSFER  COMPARISONS 

The  results  of  the  preceding  calculations  can  be  used  to  estimate  the 
drag,  skin  friction  and  heat  transfer  to  a  blunted  9°  cone.  The  total  drag 
consists  of  pressure  drag  and  skin  friction  which  are  in  turn  influenced  by 
contributions  due  to  viscous  interactions.  Other  influencing  factors  are  boun¬ 
dary  layer  transition  and  the  question  of  when  the  entropy  layer  near  the  body 
is  engulfed  by  the  boundary  layer.  Zakkay  and  Krause  have  given  a  general 
expression  for  predicting  this  "swallowing"  distance  XSw  ■  Using  their  result 
the  swallowing  distance  has  been  estimated  for  the  conditions  of  present 
interest  and  is  given  in  the  following  table. 


Table  V-8 

ENTROPY  LAYER  SWALLOWING  DISTANCE  (Xsw,  FT) 


rh 

=  0.2  IN 

«N 

=  1.5  IN 

X 

SW/rN 

CASE 

FLIGHT 

WIND  TUNNEL 

FLIGHT 

WIND  TUNNEL 

FLIGHT 

WIND  TUNNEL 

I 

2.54 

3.74 

35 

52 

278 

409 

n 

2.50 

6.1 

35 

63 

273 

494 

m 

0.15 

0. 15 

2.0 

2.3 

15.9 

18.0 

We  note  in  all  cases  the  swallowing  distance  for  the  wind  tunnel  case 
is  equal  to  or  greater  than  the  distance  for  the  corresponding  flight  case.  Also, 
in  terms  of  nose  radii  the  swallowing  distance  is  greater  than  the  range  of  the 
results  presented  in  Figs.  8-16  for  Cases  I  and  II.  Boundary  layer  calcula¬ 
tions  based  on  the  body  streamtube  results  can  be  considered  valid  for  dis¬ 
tances  less  than  the  swallowing  distance. 

Clearly,  accurate  predictions  of  the  total  skin  friction  drag  or  total 
heat  transfer  to  a  blunted  cone  must  take  account  of  the  transition  of  the  boun¬ 
dary  layer  from  laminar  to  turbulent  flow.  In  the  present  study,  however, 
our  primary  objective  is  a  comparison  of  flight  and  wind  tunnel  flow  fields. 

For  this  purpose  we  restrict  our  attention  to  the  laminar  boundary  layer  region 
only. 
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i.  Pressure  Drag.  The  effect  of  nose  bluntness  on  the  pressure  drag  of  a 

cone  is  not  dependent  on  per  se  but  rather  on  R^/R^  and  on  the  cone 

angle.  Specifying  R^/R^  for  a  given  cone  angle  8c  implies  that  the  length 

is  known.  If  we  take  a  length  of  88  inches,  then  the  nose  radii  of  0.2  in.  and 

1.  5  in.  correspond  to  R^/R^  =  0.  015  an''  0.  10.  For  a  9*  cone  and  such 

small  values  of  R../R„,  the  pressure  drag  coefficient  is  essentially  equal  to 

N  B 

that  of  a  sharp  cone^'  and  the  length  is  immaterial.  We  can  then  compare 
the  sharp  cone  pressure  drag  coefficients  for  the  three  flight  and  wind  tunnel 
cases  considered  here.  As  we  have  seen  in  Fig.  7  they  will  not  be  identical. 
The  pressure  drag  coefficients  obtained  from  the  correlations  reported  in 
Ref.  44  are  tabulated  below. 


Table  Y-9 

PRESSURE  DRAG  COEFFICIENTS 


CASE 

FLIGHT 

XI  HD  TUKHEL 

I 

0.0528 

0.0550 

n 

0.0518 

0.0540 

nr 

0.0516 

0.0517 

The  differences  between  flight  and  wind  tunnel  are  only  4%  even  when  the 
free-stream  Mach  numbers  differ  by  a  factor  of  2  (Case  II). 

2.  Viscous  Drag.  At  the  low  altitude  condition  (Cases  I  and  II)  boundary- 
layer  transition  must  be  considered.  Although  the  local  Reynolds  numbers 
and  Mach  numbers  obtained  from  the  body  streamtube  calculations  are  such 
that  one  would  not  predict  transition,  once  the  entropy  layer  is  swallowed  by 
the  boundary  layer,  transition  would  be  predicted.  The  present  calculations 
do  not  permit  determining  these  details.  One  can  make  comparisons  though 
for  the  portion  of  the  cone  where  the  boundary  layer  remains  laminar  assuming 
this  region  to  be  governed  by  the  entropy  layer. 

Taub  has  correlated  a  large  amount  of  wind  tunnel  date  for  blunted 
48 

cones  .  He  finds  that  the  drag  correlates  as  versus  TT  ,  where 
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and 


CD  "  CP/(I  +  0.315/3*)  s i**  9e 


u 


(/  +  0.-4-  TJTq  )c06  &c 
(/  +  0-37 5  f?)  sin.1  ec 


RujR* 

S  in  Qc 

Here  Cx  ”  /*-*  where  %  -  ~~  (/  +  3  T^Ta)  .  Based  on 

Taub's  correlations,  we  observe  that  for  the  present  conditions  the  drag  coef¬ 
ficients  obtained  in  the  hypothetical  wind  tunnel  would  be  the  same  as  in  flight 
if  the  wall  temperature  Tw  were  duplicated. 

3.  Skin  Friction  and  Heat  Transfer.  For  the  laminar  boundary  layer 
region  on  a  blunted  slender  cone,  the  heat-transfer  distribution  can  be  cal¬ 
culated  by  the  local  similarity  method  of  Lees^.  For  hypersonic  flow  over 
a  blunt-nosed  cone,  Lees  gives  the  ratio  of  the  local  surface  heat-transfer 
rate  to  the  stagnation-point  heat  transfer  as 


1) 

where  jp  -  the  local  static  pressure 

/  =  stagnation-point  pressure 
Ut  -  local  inviscid  flow  velocity 
Ua,  =  free-stream  velocity 

=  cross-sectional  radius  of  the  body 
S  =  distance  along  body  surface 
F$h  =  nose  radius 

d  6/t  / d.  (s/r*)0  =  velocity  gradient  at  the  stagnation  point 


The  body  streamtube  calculations  can  be  used  to  compute 
Eq.  V-l. 


from 
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The  result  for  Case  II,  which  represents  the  largest  mismatch  in 
free -stream  conditions  is  shown  in  Fig.  21a.  We  note  that  the  distributions 
for  equilibrium  and  frozen-flow  differ  by  only  15  to  20  per  cent  along  the  cone. 
The  insensitivity  to  chemistry  is  further  borne  out  by  the  closeness  of  the 
nonequilibrium  flow  curves  for  the  two  nose  radii.  Along  the  cone  wind  tunnel 
values  of  <^wj  <£.  are  about  10%  higher  than  the  flight  case. 

Use  of  Eq.  V-l  implies  a  fully  catalytic  wall.  In  this  case  there  is 
little  influence  of  chemical  reactions  in  the  boundary  layer  on  surface  heat 
transfer*.  Inger  has  shown  that  blunt  bodies  having  a  metallic-type  surface 
behave  as  very  nearly  perfectly  catalytic  at  all  altitudes  below  300,  000  ft  for 
U—  =  26,  000  fps.  Therefore,  the  present  results  are  deemed  reasonable. 


Skin-friction  distributions  can  be  obtained  from  the  heat-transfer 
distributions  through  the  modified  Reynolds  analogy  relation 


C-h 


2 


%  Pr 


-Z/3 


(V-2) 


i 

!■ 

=• 


where 


£  h  ~  Quj  f  ft  (  PtuiiM  ~  Hw) 

Cf  =  X  /  7*  p, 

Pr  =  Prandtl  number 

pc  =  inviscid-flow  density  at  edge  of  boundary  layer 
P(A.uj-  Pz  *  Pr  *  Ue.  ,  the  adiabatic-wall  enthalpy 

Hu  =  static  enthalpy  at  the  wall 
-  shear  stress  at  the  wall 


Computations  of  C ^  based  on  Eq.  V-2  and  the  results  contained  in  Fig.  21a 
are  shown  in  Fig.  21b  for  Case  II.  The  frozen  and  equilibrium  flow  distribu¬ 
tions  are  shown  only  for  the  flight  case  with  =  0.  2  in.  The  corresponding 
distributions  for  R  =  1.  5  in.  can  be  obtained  from  these  curves  by  multi¬ 
plying  by  the  ratio  of  the  stagnation-point  heat-transfer  rates.  The  difference 
between  the  equilibrium  and  frozen-flow  curves  varies  between  15%  and  25% 
along  the  body.  The  nonequilibrium  wind-tunnel  results  are  in  good  agreement 
with  the  corresponding  flight  case  for  both  nose  radii.  In  these  calculations 
the  wall  temperatures  have  been  taken  as  1000  °K  for  the  flight  cases  and 
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300  °K  for  the  wind  tunnel  cases.  Thus,  the  sensitivity  to  wall  temperature 
is  weak. 

The  stagnation-point  heating  rates  calculated  for  Case  II  are  given  in 
the  following  table. 


Table  1-10 

STAGNATION-POINT  HEAT  TRANSFER 
Case  n  -  25,000  fps  80,000  ft  Alt. 


Rh  (IN.) 

%  ( 

btu)/ft2-sec 

FLIGHT 

WIND  TUNHEL 

0.2 

19,700 

21,900 

1.5 

7340 

8160 

The  higner  heat  transfer  for  the  wind  tunnel  is  a  result  of  the  higher  total 
enthalpy  required  to  provide  velocity  and  density  altitude  duplication.  It  was 
shown  in  Ref.  1  that  duplicating  total  enthalpy  and  density  altitude  leads  to 
smaller  differences  in  stagnation-point  heat  transfer.  Recalling  that 
for  the  wind  tunnel  case  was  about  10%  higher  than  the  flight  case  (Fig.  21a), 
the  higher  values  of  in  the  wind  tunnel  will  result  in  actual  surface  heating 

rates  that  are  about  22%  higher  in  the  wind  tunnel  than  in  flight. 

E.  SUMMARY 

In  this  section,  a  comparison  has  been  made  of  wind  tunnel  and  flight 
flow  fields  about  a  blunted  slender  cone  at  three  flight  conditions.  The  effects 
of  flow-field  nonequilibrium  and  nose  bluntness  have  been  investigated  as  well 
as  the  effects  cf  appreciable  differences  in  free-stream  pressure,  temperature 
and  Mach  number.  A  hypothetical  wind  tunnel  performance  (Fig.  2)  was  used 
to  define  test-section  conditions.  The  comparisons  have  been  restricted  to 
the  body  streamtube,  to  conditions  behind  the  conical  portion  of  the  bow  shock 
wave,  and  to  other  results  derivable  from  these  flow  conditions.  It  has  been 
assumed  that  the  flight  free-stream  density  and  velocity  W'.-  duplicated  in  the 
hypothetical  wind  tunnel  and  that  the  bodies  were  geometrically  identical  -  both 
requirements  of  the  Mach  number  independence  principle  when  dissapative 
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effects  are  presect. 

Because  p.  and  were  duplicated,  and  since  the  free -stream  Mach 
numbers  were  fairly  large,  the  Mach  number  independence  principle  indicates 
that  the  flow  in  the  region  of  the  blunt  nose  should  be  the  same  for  both  the 
flight  and  wind  tunnel  cases.  A  comparison  of  the  stagnation-point  conditions 
and  the  distributions  of  the  various  quantities  along  the  spherical  nose  indicates 
that  this  is  indeed  the  case. 

In  reviewing  Cheng's  extension  of  hypersonic  small-disturbance  simil¬ 
itude  to  real  gases  (Sec.  IV),  we  have  pointed  out  that,  when  the  free-stream 
thermodynamic  state  is  not  duplicated,  it  is  preferable  to  replace  Mach  number 
with  free-stream  velocity  in  the  similitude  parameters.  Thus,  in  the  present 
study,  U*,  sin  9C  and  d-v/x  are  duplicated  rather  than  sin  &c 

and  Ku  “  "JcDJ  d-ufa  as  given  by  Cheng.  This  has  resulted  in  relatively 
small  differences  in  the  flight  and  wind-tunnel  flows  along  the  conical  afterbody 
when  the  quantities  are  referenced  to  stagnation-point  conditions  rather  than 
free-stream  conditions.  The  overall  effects  of  the  large  mismatches  in  -pm  , 
and  are  to  produce  25  percent  to  35  percent  differences  in  the 

pressure  ratio  jo  J jo',  ,  up  to  8  percent  differences  in  temperature  ratio 
T/Tb  and  up  to  20  percent  differences  in  density  ratio  p /$'o  at  a  given 
distance  along  the  conical  afterbody  (Figs.  7-9,  12,  13,  17,  18). 

It  may  be  concluded  that  the  entropy-layer  flow  is  relatively  insensi¬ 
tive  to  mismatches  in  jom  ,  7 ^  and  M ^  .  On  the  other  hand,  it  has  been 
shown  that  conditions  behind  the  conical  shock  wave  are  much  more  sensitive 
to  mismatches  in  the  above  parameters  for  the  cases  considered  here. 

Finally,  the  drag,  skin  friction  and  heat  transfer  for  the  case  of  a 
/ 

laminar  boundary  layer  growing  within  the  entropy  layer  have  been  investigated. 
Again  the  flight  and  wind  tunnel  solutions  show  relatively  small  differences; 
however,  the  turbulent  boundary-layer  case,  which  most  likely  is  governed 
by  the  shock  layer  flow  outside  the  entropy  layer,  has  not  been  investigated. 
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SECTIOH  VI 

ANALYSIS  OF  THE  COMBINED  EFFECTS  OF  NOSE  BLUNTNESS 
AND  NONEQUILIBRIUM  ON  SLENDER-BODY  FLOWS 


A.  INTRODUCTORY  REMARKS 

The  purpose  of  the  present  section  is  to  develop  a  theoretical  model 
for  the  study  of  the  combined  effects  of  nose  bluntness  and  nonequilibrium 
on  slender-bcdy  flows  The  studies  in  this  section  are  concerned  with  flight 
conditions  at  higher  altitudes  where  the  assumed  performance  of  the  hypo- 
che;ical  wind  tunnel  under  consideration  is  such  that  it  could  fully  duplicate 
the  ambient  altitude  thermodynamic  state  if  there  were  no  nozzle  nonequi¬ 
librium  effects.  The  major  part  of  the  present  analysis  is  confined  to  the 
treatment  of  inviscid  flows;  however,  viscous  effects  are  subsequently 
accounted  for  through  a  boundary- layer  approach. 

One  of  the  primary  objectives  of  the  present  analysis  is  to  check  the 

influence  of  nonequilibrium  effects  on  the  pressure  distributions  on  blunt 

nosed  slender  bodies.  It  will  be  recalled  that  in  the  phase  I  part  of  the 

present  study  (which  was  concerned  with  flight  conditions  below  the  assumed 

full -duplication  boundary),  it  was  assumed  that  the  pressure  distribution 

was  insensitive  to  real-gas  effects.  However,  in  the  flight  regime  of  present 

51 

interest  this  assumption  is  questionable.  For  example,  Whalen  studied 
the  case  of  a  wedge -nosed  flat  plate  and  found  that  freezing  through  the 
Prandtl-Meyer  expansion  would  produce  pressures  on  the  afterbody  that 
were  significantly  lower  than  the  pressures  for  the  equilibrium  or  perfect- 
gas  cases. 

B.  REVIEW  OF  BLUNT  NOSED  SLENDER  BODY  FLOWS 

We  will  now  review  the  existing  theoretical  analyses  for  predicting 

the  flow  fields  around  blunt  nosed  slender  bodies.  However,  no  exhaustive 

review  will  be  attempted  here  since  several  comprehensive  reviews  of  the 

52 

various  aspects  of  the  problem  are  available  (Mirels  ,  Guiraud,  Vallee 
and  Zolver  ,  Hayes  and  Probstein"^,  Cox  and  Crabtree^,  Cheng^, 
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Chernyi^,  Sedov^,  Stanyukovich  Probstein^  and  Lees^).  On  the 
other  hand  we  will  revir-w  the  analyses  that  are  directly  relevant  tc  the 
present  study  in  some  detail. 


Most  analytical  studies  of  blunt  nosed  slender  bodies  are  based  on 

either  an  explicit  or  implicit  use  of  the  hypersonic  equivalence  principle  of 
24 

Hayes  .  In  general,  use  is  also  made  of  the  hypersonic  small  disturbance 

61  23  62 

approximations  *  and  the  Newtonian  slender  body  approximations 

analogy  between  the  flows  past  blunted  slender  bodies  and  the  unsteady  flows 


The 


due  to  strong  explosions  was  pointed  out  independently  by  Cheng  and  Pailcne 

54 

and  Lees  and  Kubota  .  These  authors  obtained  the  solutions  for  steady 
hypersonic  flows  past  blunted  flat  plates  and  cylinders  by  applying  the  equiv¬ 
alence  principle  to  the  corresponding  unsteady  flow  problems  of  explosions 

with  cylindrical  and  spherical  symmetries.  The  solutions  to  unsteady  prob- 
,  65  66  67  68 

lems  had  been  obtained  earlier  by  L,in  ,  Taylor  ,  Sedov  and  Sakurai 

Other  authors^  have  studied  the  unsteady  'similar1  flows  corresponding 
to  shock  waves  propagating  according  to  various  power  laws  and  the  corre¬ 
sponding  steady  flows  past  equivalent  power-law  bodies.  An  approximate 
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method  of  solving  for  flows  past  power-law  bodies  has  been  given  by  Mirels 

73 

and  Mirels  and  Thornton  have  investigated  bodies  associated  with  shock 
shapes  which  deviate  slightly  from  a  power  law. 


72 


The  method  of  generating  solutions  for  flows  past  power-law  bodies 
from  equivalent  unsteady,  similar  flows  is  not  applicable  to  the  flows  past 
blunted  wedges  and  cones  since  these  flows  are  basically  nonsimilar.  That 
is,  the  time  evolution  of  the  flows  in  the  unsteady  cases,  and  the  space  evolu¬ 
tion  of  the  flows  in  the  steady  cases,  do  not  obey  similarity  requirements. 

7^4  76  3  1 

Chernyi  *  and  Cheng  *  have  developed  approximate  methods  to  study 

the  flows  past  blunted  wedges  and  cones. 


The  method  used  by  Chernyi  is  an  integral  one  in  which  the  conservation 
relations  of  mass,  momentum  and  energy  are  applied  in  an  integral  form  across 
the  disturbed  region  between  the  b  .'dy  a:  d  the  shock  wave.  The  effect  of  the 
blunt  nose  on  the  flow  downstream  of  it  is  assumed  to  be  describable  in  terms 
of  an  energy  in  the  energy  conservation  equation  and  an  impulse  in  the  norhnal 
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momentum  conservation  equation.  He  further  assumes  that  the  bulk  of  the 

disturbed  gas  is  concentrated  in  a  thin  region  next  to  the  shock  wave,  and 

obtains  approximate  solutions  for  the  flows  past  blunted  wedges  and  cones. 

77  78 

Chernyi  ’  has  also  described  general  integral  formulations  applicable 
to  hypersonic  flows  containing  strong  shock  waves. 

The  approach  used  by  Cheng^*  is  somewhat  different  in  that  he  accounts 
for  the  so-called  :entropy  layer’  in  a  specific  way.  The  entropy  layer  consists 
of  those  streamlines  which  have  crossed  the  nearly  normal  portions  of  the 
shock  wave  and  are  associated  with  much  larger  values  of  entropy  than  the 
streamlines  which  have  crossed  the  more  oblique  (and  therefore  weaker) 
parts  of  the  shock  wave.  By  making  estimates  of  the  relative  orders  of  mag¬ 
nitudes  of  the  various  flow  quantities  within  and  outside  the  entropy  layer, 
Cheng  was  able  to  reduce  the  partial  differential  equations  governing  the 
problem  to  a  single  ordinary  differential  equation  relating  the  pressure  dis¬ 
tribution  on  the  body  to  the  cross-sectional  area  of  the  entropy  layer.  This 
differential  equation  can  be  solved  for  the  shape  of  the  outer  edge  of  the 
entropy  layer  (and  the  body  pressure  distribution)  by  assuming  that  the  pres¬ 
sure  distribution  is  related  to  the  shape  of  the  entropy-layer  edge  through  the 
Newton-Busemann  pressure  distribution  law.  Thus  Cheng's  approach  clearly 
brings  out  the  fact  '.nat  the  influence  of  a  blunt  nose  on  the  flow  downstream 
of  it  appears  essentially  as  a  displacement  effect. 

The  analyses  of  Chernyi  and  Cheng  are  both  based  on  thin-shock-layer 
concepts,  that  is,  both  these  methods  utilize  the  fact  that  a  major  portion  of 
the  disturbed  gas  is  concentrated  in  a  thin  region  close  to  the  shock  wave. 

The  thin-shock-iayer  approximation  becomes  progressively  more  accurate 
as  the  specific  heat  ratio,  1  ,  of  the  gas  approaches  unity.  In  fact,  as  far 
as  the  leading  approximation  for  small  e  =  (/  -  l)  j  (V  +  1)  and  (M0)  ^  is  con¬ 
cerned,  the  two  analyses  are  identical;  they, differ  only  when  higher-order 
approximations  are  considered.  There  seems  to  be  some  disagreement  in 

the  literature  regarding  the  relative  accuracies  of  the  two  methods.  Thus 
3 1 

Cheng  states  that  Chernyi's  method  is  unlikely  to  lead  to  consistent  higher- 

54 

order  approximations,  while  Hayes  and  Probstein  argue  that  Cheng's  method 
will  be  less  accurate  (than  Chernyi's)  since  he  uses  the  Newton-Busemann 
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pressure  formula  in  terms  of  the  entropy-layer  edge  rather  than  the  shock 
shape.  The  approach  we'^hall  adopt  for  the  present  analysis  will  be  a  hybrid 
one  incorporating  features  from  both  methods. 


79 

A  method  similar  to  that  of  Chernyi  has  also  been  used  by  Luniev 
80-83 

Several  authors  have  also  investigated  the  validity  of  the  'blast-wave 

analogy'.  The  requirements  for  the  validity  of  the  equivalence  principle  fail 

not  only  in  the  vicinity  of  the  blunt  nose^but  also  in  the  high-entropy  layer 

next  to  the  body.  Thus,  when  the  effect  of  a  blunt  nose  on  the  flow  downstream 

of  it  is  replaced  by  that  of  a  concentrated  energy  release  at  the  nose,  the 

density  and  temperature  distributions  near  the  body  are  predicted  incorrectly. 

31  r 

However,  within  the  accuracy  of  the  assumptions  made  by  Cheng  ,  the 
velocity  and  pressure  fields  are  predicted  satisfactorily  throughout  the  whole 
flow  field. 


All  the  analyses  described  above  are  applicable  only  for  the  perfect- 

gas  case,  that  is,  for  the  case  in  which  the  specific  heat  ratio  is  constant 

throughout  the  whole  flow  field.  There  are  very  few  analyses  in  the  literature 

which  are  applicable  to  the  real-gas  flows  past  blunt-nosed  slender  bodies. 

25 

Cheng  has  described  a  similitude  for  equilibrium  real-air  flows  past  slender 

bodies,  and  this  similitude  has  been  extended  to  nonequilibrium  flows  by 
26 

Inger  .  Analyses  which  are  valid  for  special  geometries  and  simplified 

84  51 

flow  chemistries  have  been  reported  by  Luniev  ,  Whalen  and  Vagalio- 
,  85 

Launn  and  Bloom 


Analyses  of  the  combined  effects  of  nose  bluntness  and  boundary  layer 

displacement  effects  have  been  reported  by  Cheng^*,  Oguchi^  and  others. 

The  similitude  requirements  for  viscous  slender  body  flows  with  and  without 

30 

nose  bluntness  effects  have  been  considered  by  Hayes  and  Probstein  , 

87  88  3 1 

Luniev  '  and  Cheng  .  Numerical  solutions  for  the  blunt-leading-edge 

8  9 

problem  have  been  reported  by  Van  Hise  '  for  the  perfect-gas  case,  by 

90  40 

Feldman  for  the  equilibrium  flow  case  and  by  Curtis  and  Strom  for  the 

nonequilibrium  flow  case. 
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C.  F ORM U LA T ION  OF  A  THEORETICAL  MODEL 

In  the  present  section  we  wilt  formulate  a  theoretical  r  odel  to  study 
the  'cbmbined  effects  of  nose  bluntness  and  flow  nonequilibrium  slender 
body  flow  fields.  The  model  described  is  considered  to  be  the  o..e  most- 
likely  to  further  qualitative  understanding  of  the  problem  and  at  the  same 
time  give  reasonable  quantitative  answer.  For  the  inviscid  case,  following 
the  practice  in  the  perfect-gas  analyses,  the  flow  field  downstream  of  the 
blunt  nose  is  divided  into  an  inner  low-density  entropy  layer  near  the  body 
and  an  outer  shock  layer.  The  entropy  layer  consists  of  streamlines  cross¬ 
ing  the  nearly  normal  portion  of  the  shock  wave  in  the  vicinity  of  the  blunt 
nose,  while  the  shock  layer  consists  of  streamlines  crossing  the  more 
oblique  and  weaker  parts  of  the  shock  wave.  It  is  clear  that  becatse  of 
the  relatively  larger  values  of  the  entropy  and  lower  values  of  the  densi- / 
in  the  entropy  layer,  nonequilibrium  effects  are  more  likely  to  be  important 
there.  In  the  proposed  simplified  model  the  nonequilibrium  effects  will  be 
accounted  for  only  in  the  entropy  layer,  and  the  outer  shock  layer  will  be 
assumed  to  be  described  by  the  classical  Newtonian  shock  layer  theory. 

The  accuracy  of  the  assumption  that  real-gas  effects  are  more 
important  in  the  entropy  layer  than  in  the  shock  layer  can  be  verified  as 
follows:  For  the  sake  of  simplicity  let  us  consider  the  case  when  the  flow 
is  assumed  to  be  in  equilibrium  everywhere.  Then,  an  estimate  of  the 
relative  importance  of  real-gas  effects  in  the  entropy  and  shock  layers  can 
be  obtained  by  comparing  the  dissociation  fractions  at  the  body  and  immed¬ 
iately  behind  the  shock  wave  for  various  axial  locations.  For  the  purpose 
of  the  present  study,  the  dissociation  fraction,  oc  ,  for  air  can  be  assumed 
to  be  equal  to  Z-l  where  Z  is  the  compressibility  factor. 

Consider  the  blunt-nosed  slender  body  shown  in  Fig.  22.  The  notation 
that  we  shall  adopt  is  also  shown  on  this  figure.  All  the  flow  quantities  in 
the  nose  region  will  be  denoted  by  the  subscript,  ’n1  while  all  the  quantities 
on  the  body  and  immediately  behind  the  shock  wave  will  be  denoted  by  the 
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subscripts  ‘b1  and  's*,  respectively.  The  symbol  ‘e*  will  be  used  to  denote 
the  edge  of  the  entropy  layer  while  the  subscript  '&>'  will  be  used  to  indicate 
free-stream  values. 

Now,  because  of  the  equilibrium  as  sumption, the  flow  will  be  particle 
isentropic,  and  the  entropy  on  the  body  will  be  that  corresponding  to  con¬ 
ditions  behind  a  normal  shock  wave  at  the  same  free-stream  conditions. 

Also,  for  order  of  magnitude  estimates,  we  may  assume  that  the  pressure 
on  the  body  is  of  the  same  order  as  that  behind  the  shock  wave  at  the  same 
axial  location,  i.  e.  -f>k  *  o(j bj).  Hence,  for  a  given  set  of  free-stream  con¬ 
ditions,  an  estimate  of  the  dissociation  fraction  on  the  body  downstream  of 
the  blunt  nose  can  be  obtained  by  reading  off  a  Mollier  diagram  for  air  the 
value  of  oc b  corresponding  to  the  pressure  -f»s  and  the  entropy  . 
Figure  23  shows  a  comparison  of  the  dissociation  fraction  on  the  body  and 
behind  the  oblique  shock  wave  for  various  axial  locations  (corresponding  to 
different  values  of  the  local  shock  angle)  and  for  various  free-stream  con- 
ditions.  It  can  be  seen  from  the  figure  that  the  dissociation  fraction  on  the 
body  is  always  much  larger  than  the  dissociation  fraction  behind  the  oblique 
shock  wave,  and  that  the  dissociation  levels  behind  the  shock  wave  reach 
significant  levels  only  at  very  high  velocities  and  altitudes.  Even  at  the 
higher  velocities  and  altitudes  it  can  be  seen  that  the  dissociation  levels  in 
the  shock  wave  are  small  for  local  shock  inclinations  of  20°  or  less.  Thus, 
for  slender  bodies  the  real-gas  effects  in  the  shock  la^er  will  be  small  a 
few  nose  radii  downstream  of  a  blunt  nose.  On  the  other  hand,  real-gas 
effects  will  always  be  relatively  more  important  in  the  entropy  layer  (for 


#  We  will  not  use  any  explicit  specification  of  the  edge  of  the  entropy  layer, 
since  this  is  not  necessary  within  the  context  of  the  present  analysis. 

Cox  and  CrabtreeSS  define  the  edge  of  the  entropy  layer  as  being  composed 
of  that  streamline  whose  local  inclination  after  crossing  the  shock  wave  is 
approximately  20°.  Other  definitions  have  also  been  used  in  the  literature. 

**The  equilibrium  normal  shock  wave  parameters  reported  in  Ref.  91  were 
used  in  these  calculations.  Ref.  91  uses  the  1959  ARDC  model  atmosphere. 
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the  ranges  of  free-stream  conditions  under  consideration)  even  for  small 
values  of  the  local  shock- wave  inclination. 

It  should  be  pointed  out  that,  even  though  the  above  estimations  were 
made  only  for  the  case  of  fully  equilibrium  flow,  the  conclusions  will  be 
valid  even  when  the  flow  is  not  in  equilibrium.  This  is  because  in  nonequi¬ 
librium  flow  the  dissociation  fractions  on  the  body  will  be  higher  than  the 
corresponding  equilibrium  values  since  the  relaxation  process  here  is  essen¬ 
tially  one  of  recombination.  Therefore,  at  any  axial  location  downstream 
of  the  blunt  nose,  the  dissociation  fractions  in  regions  close  to  the  body  will 
always  be  higher  than  those  in  regions  close  to  the  shock  wave.  Henceforth 
we  will  assume  that  nonequilibrium  effects  are  important  only  in  the  entropy 
layer  aqd^we  will  approximate  the  region  outside  the  entropy  layer  by  a 
constant  specific  heat  ratio,  ideal-gas  region. 

Before  we  can  proceed  with  estimations  of  the  relative  magnitudes 
of  the  other  flow  variables  within  the  entropy  and  shock  layers,  it  is  neces¬ 
sary  to  formulate  a  simplified  model  for  the  flow  chemistry.  Because  of 
the  extremely  complicated  nature  of  physical-chemical  processes  involved 
in  hypersonic  flows  past  even  such  simple  bodies  as  spherically  blunted 
slender  cones,  some  simplifying  assumptions  have  to  be  resorted  to.  Several 
su  ;h  simplifying  assumptions  have  been  used  in  the  literature.  One  of  the 
more  commonly  used  approximations  is  to  assume  that  the  flow  remains 
effectively  in  equilibrium  along  each  streamline  up  to  some  suitably  chosen 
boundary  and  that  the  gas  composition  'freezes'-  (and  remains  frozen)  sub¬ 
sequently  along  each  streamline  at  the  value  coi responding  to  that  at  the 
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boundary.  For  example,  Whalen  has  studied  the  hypersonic  flow  past  a 
wedge -nosed  flat  plate  by  assuming  that  the  gas  composition  is  constant 
along  each  streamline  and  is  given  by  the  equilibrium  value  immediately 
behind  the  shock  envelope.  In  Refs.  1  and  2  the  influence  of  real-gas  effects 
on  the  sensitivity  of  various  afterbody  flow  quantities  of  a  sharp-nosed 
slender  body  to  mismatches  in  free -stream  conditions  was  studied  by 
assuming  that  the  flow  was  in  equilibrium  in  the  nose  region  and  subsequently 
freezes  at  the  junction  of  the  nose  and  the  afterbody.  In  Ref.  1  this  model 
was  termed  che  equilibrium-frozen  model. 
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When  there  exists  a  sharp  discontinuity  in  the  body  slope  in  the 
shoulder  region  between  the  nose  and  the  afterbody,  then  it  is  logical  to 
assume  that  the  boundary  between  the  equilibrium  and  frozen  flow  regions 
passes  through  this  sharp  corner  since  the  flow  experiences  a  rapid  expan¬ 
sion  in  this  region.  Thus,  for  the  flows  past  wedge-nosed  flat  plates  or 
conical-nosed  cylinders,  one  can  assume  that  the  flows  remain  in  thermo¬ 
dynamic  equilibrium  upstream  of  the  expansion  fan  emanating  from  the 
sharp  shoulder  junction  and  that  they  freeze  downstream  of  it.  Such  a 
logical  choice  does  not  exist  when  there  is  no  sharply  defined  boundary 
between  the  nose  and  the  afterbody  as  is  the  case,  for  example,  with  spher¬ 
ically  blunted  cones. 

85 

Vaglio-Laurin  and  Bloom  have  studied  the  influence  of  real-gas 
effects  on  the  flowfields  past  several  typical  body  shapes  by  using  a  sim¬ 
plified  model  for  the  flow  chemistry.  In  this  model,  the  authors  assume 
that  the  flow  chemistry  remains  frozen  downstream  of  the  forward  facing 
characteristic  that  intersects  the  body  at  the  sonic  point;  upstream  of  this 
characteristic  the  flow  is  assumed  to  remain  in  chemical  equilibrium. 
Vaglio-Laurin  and  Bloom  also  assume  that  the  gas  composition  remains 
constant  (at  the  values  corresponding  to  the  equilibrium  values  behind  the 
shock  wave)  along  the  streamlines  that  cross  the  shock  wave  downstream 
of  the  characteristic  mentioned  above.  With  the  above  assumptions,  Vaglio- 
Laurin  and  Bloom  have  integrated  the  simplified  differential  equations  gov¬ 
erning  the  problem  by  using  the  method  of  characteristics,  and  they  have 
obtained  the  distributions  of  various  flow  quantities  on  the  body  as  well  as 
the  shock  layer. 

In  support  of  their  sonic-point-freezing  assumption,  Vaglio-Laurin 

and  Bloom  point  out  that  for  typical  reentry  configurations  flying  at  high 

velocities  and  at  altitudes  in  excess  of  125,  000  feet  the  flows  indeed  show 

a  strong  tendency  to  freeze  at  or  near  the  sonic  point  on  the  body.  Moreover 

90 

exact  numerical  calculations,  such  as  the  ones  performed  by  Feldman  for 
an  equilibrium  flow  over  a  blunt  nosed  slender  body,  indicate  that  the  gas 
close  to  the  bedy  at  the  sonic  point  expands  drastically  as  it  flows  around  the 
body.  In  fact,  Feldman  gives  an  example  in  which  a  streamline  which  is 
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only  a  distance  of  0.  012  nose-radii  away  from  the  body  upstream  of  the  sonic 
line  eventually  ends  up  2  nose-radii  away  from  the  body  when  far  downstream 
of  the  nose.  This  drastic  expansion  around  and  beyond  the  sonic  region  will, 
of  coux-aq^tend  to  freeze  the  flow  chemistry  along  streamlines  close  to  the 
body  beyond  the  sonic  point. 

In  the  present  analysis  we  will  assume  that  the  flow  chemistry  within 

the  entropy  layer  freezes  (and  remains  frozen)  downstream  of  the  sonic  point 

station.  It  should  be  recalled  here  that  in  the  present  model  the  flow  outside 

the  entropy  layer  is  assumed  to  behave  like  an  ideal  gas  (with  a  specific  heat 

ratio  equal  to  that  of  the  free  stream).  On  the  other  hand,  in  the  models  used 
51  85 

by  Whalen  and  Vaglio-Laurin  and  Bloom  ,  the  authors  assume  that  in 
addition  to  the  flow  being  frozen  downstream  of  a  specified  boundary,  the 
gas  composition  remains  constant  (at  the  values  corresponding  to  equilibrium 
values  immediately  behind  the  shock)  along  all  streamlines  that  cross  even 
the  relatively  weaker  portions  of  the  shock  wave.  This  difference  between 
the  assumptions  of  the  above  authors  and  those  of  the  present  model  is 
important  since  it  will  be  crucial  to  the  analysis  that  follows. 

Having  formulated  a  simplified  model  for  the  flow  chemisti /,  we 
can  now  proceed  to  make  estimations  of  the  various  flow  quantities  within 
the  entropy  layer.  One  of  the  crucial  assumptions  involved  in  the  so-called 
thin-shock-layer  theories  is  that  the  density  within  the  entropy  layer  in 
general,  and  on  the  body  in  particular,  is  much  smaller  than  the  free- stream 
density.  This  assumption  can  be  checked  for  the  present  model  in  fairly 
general  terms  without  any  specific  considerations  of  the  nose  or  body  shapes. 

Now,  for  a  given  set  of  free -stream  conditions,  the  entropy  behind 

the  shock  wave, as  well  as  the  pressure  at  the  stagnation  point, can  be  deter- 
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mined  from  equilibrium-normal- shock  tables  .  Also,  since  the  flow  in  the 
nose  region  is  assumed  to  be  in  equilibrium  and  since  the  pressure  distri¬ 
bution  in  the  nose  region  is  fairly  insensitive  to  real-gas  effects,  one  can 
calculate  the  pressure  at  the  sonic  point  by  using  the  isentropic-expansion 
relations  with  the  specific  heat  ratio,  /  equal  to  1.4.  That  is,  one  can  use 
the  relation  =  0.  528  Pst,  where  the  subscripts  'p 1  and  ' st 1  refer  to  the 
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sonic  point  and  stagnation  point,  respectively.  Sin_a  the  entropy  and  the 
pressure  at  the  sonic  point  are  known, all  other  state  variables  can  be  easily 
read  off  a  Mollier  diagram;  in  particular,  the  density  and  the  effective  dis¬ 
sociation  fraction  at  the  sonic  point  c^a  be  determined. 

In  the  present  model  the  flow  within  the  entropy  layer  is  assumed  to 
freeze  downstream  of  a  line  passing  through  the  sonic  point  on  the  body.  Thus 
the  expansion  along  the  body  surface  will  take  place  with  an  (constant)  effective 
specific  heat  ratio  determined  by  the  conditions  at  the  sonic  point.  Following 
Whalen  ,  this  effective  specific  heat  ratio  will  be  taken  to  be 
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The  density  distribution  along  the  body  corresponding  to  a  given  pressure 
distribution  can  then  be  determined  by  using  the  isentropic  relation 


h  J 


(VL  2) 


Fig.  24  shows  a  plot  of  the  variation  of  density  with  the  free-stream 
velocity  (at  an  altitude  of  200,  000  ft)  for  the  station  on  the  body  at  which 
-f>b  .  It  can  be  seen  from  the  figure  that  the  density  at  this  station 

is  always  considerably  smaller  than  the  free-stream  density.  The  small¬ 
ness  of  the  density  on  the  body  is,  of  course,  a  direct  consequence  of  the 
large  (compared  to  free  stream)  value  of  the  entropy  along  the  body  stream¬ 
line.  For  the  same  reason,  the  density  within  the  entire  entropy  layer  will 
also  be  much  smaller  than  the  free-stream  value. 

It  follows  immediately  that  the  low-density  entropy  layer  will  be 
incapable  of  supporting  any  transverse  pres,  are  gradients,  and  that  the 
pressure  at  any  axial  station  within  the  entropy  layer  can  be  taken  as  equal 
to  the  value  at  the  outer  edge  of  the  entropy  layer.  It  also  follows  from 
simple  mass-flow  arguments  that  the  cross-sectional  area  of  the  shock  layer 
will  be  much  smaller  than  that  of  the  entropy  layer.  These  are  the  criteria 
that  are  necessary  for  the  successful  application  of  thin- shock-layer  methods 
for  the  present  problem. 
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Before  the  proposed  model  can  be  utilized  to  study  the  problem,  it  is 

also  necessary  to  reexamine  the  validity  of  the  equivalence  principle  for  the 

present  case,  since  the  assumption  of  its  validity  is  implicit  in  both  the  methods 

of  Cheng  and  Chernyi  .  Inger  has  pointed  out  that  the  equivalence 

principle  has  very  general  validity  within  the  framework  of  the  assumptions 

of  small-disturbance  theory,  and  that  it  is  valid  for  nonequilibrium  flows 

irrespective  of  the  form  of  the  equation  of  state  and  the  reaction  rates.  The 

small-disturbance  approximations  are,  of  course,  based  on  the  assumption 

that  u  x  U- &  and  ^  <r  ~  cr  However,  as  has  been  pointed 
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out  by  Whalen  ,  the  validity  of  the  former  of  the  two  assumptions  is  depen¬ 
dent  on  the  value  of  the  dissociation  fraction,  since  the  energy  conservation 
relation  can  be  written  in  the  form 


(VI.  3) 


Thus  the  assumption  that  u  %  will  be  valid  only  in  those  regions 

of  the  flow  in  which  the  energy  contained  in  the  dissociation  mode  is  small 
compared  to  the  free -stream  kinetic  energy. 

Since  in  the  present  model  real-gas  effects  are  not  considered  in  the 

shock  layer,  the  assumption  u  ~  u  ^  will  be  valid  within  the  shock  layer 

to  the  same  degree  of  approximation  that  it  is  valid  in  the  corresponding 
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ideal-gas  case.  Cheng  has  shown  that  for  the  ideal-gas  case  the  assump- 

/  ■? 

tion  u  x  u  will  be  valid  if  -L-  cr  ' - - J 


Is  negligible  compared  to 


unity,  where  <r  can  be  a  typical  shock  or  flow  angle.  In  the  present  model, 
the  above  requirement  should  also  ensure  the  validity  of  the  equivalence 
principle  within  the  shock  layer.  However,  it  's  clear  that  in  regions  close 
to  the  wall  the  approximation  u  ~  um  cannot  be  valid,  since  the  dissoci¬ 
ation  fractions  in  these  regions  need  not  be  small.  However,  this  per  se 
does  not  introduce  additional  complications,  since  the  equivalence  principle 
fails  in  regions  close  to  the  body  even  in  the  ideal-gas  case.  However,  as 
is  well  known,  the  pressure  field  on  the  body  and  the  shock  shape  are  still 
predicted  correctly  in  spite  of  this  local  failure  of  the  equivalence  principle. 

To  summarize  briefly,  the  proposed  model  is  such  that  real-gas  effects 
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are  included  when  calculating  the  displacement  effect  of  the  nose  bluntness 
and  the  resulting  effective  body  shape,  but  not  in  calculating  the  evolution  of 
the  shock  wave  due  to  this  effective  body  shape.  In  other  words,  the  outer 
edge  of  the  entrcpy  layer  is  related  to  the  body  shape  including  real-gas 
effects,  but  the  shock-wave  shape  corresponding  to  the  outer  edge  of  the 
entropy  layer  is  calculated  using  the  classi  .i  Newton-Busemann  relations. 

The  various  approximations  and  restrictions  that  apply  to  the  present 
analysis  are  listed  below. 

a)  The  flow  outside  the  entropy  layer  is  assumed  to  behave  like  that  of 
a  perfect  gas  with  a  specific  heat  ratio  equal  to  the  free- stream  value.  The 
flow  inside  the  entropy  layer  is  particle  isentropic,  but  the  specific  heat 
ratio  is  different  on  different  streamlines. 

b)  The  flow  is  assumed  to  remain  in  chemical  equilibrium  in  the  nose 
region  and  subsequently  freeze  along  a  line  passing  through  the  sonic  point 
on  the  body. 

c)  The  pressure  distribution  on  the  blunt  nose  is  assumed  to  be  given  by 
the  modified  Newtonian  pressure  formula. 

d)  The  free-stream  is  assumed  to  be  homogeneous  in  flow  properties 
and  undissociated. 

e)  The  local  streamline  inclination  to  the  free-stream  direction  is 
assumed  to  be  small  everywhere  except  in  a  small  region  near  the  nose. 

f)  The  vibrational  degree  of  freedom  is  assumed  to  remain  in  equilib¬ 
rium  with  the  translational  temperature  within  the  entropy  layer. 

g)  The  ratio  of  specific  heats  is  assumed  to  be  sufficiently  close  to  unity 

y  j 

so  that  terms  of  order  — - —  can  be  neglected. 
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A  RESTRICTED  SIMILITUDE  APPLICABLE  TO  REAL-GAS  FLOWS 
PAST  BLUNTED  SLENDER  BODIES 


As  pointed  out  earlier,  a  similitude  applicable  to  equilibrium  flows 

25 

past  blunt-nosed  slender  bodies  was  given  by  Cheng  ,  and  this  was  extended 

26 

to  nonequilibrium  flows  by  Inger  .  The  similitude  developed  by  Cheng 
requires  that  for  proper  similitude  between  two  flows  past  two  geometrically 
similar  bodies,  the  hypersonic  similarity  parameter,  u w  T  ,  a  bluntness 
parameter, as  well  as  the  free-stream  thermodynamic  strte  have  to  be  the 
same.  The  last  of  these  conditions  is  quite  stringent,  and  the  free-stream 
thermodynamic  state  corresponding  to  low-altitude,  hypervelocity  flight  can¬ 
not  be  reproduced  in  existing  test  facilities,  except  perhaps  in  such  advanced 
(and  as  yet  in  a  developmental  stage)  facilities  as  the  Isentropic  Compression 
Tube  ’  .  Therefore,  it  would  be  extremely  useful  if  even  a  restricted 

similitude  can  be  developed  in  which  the  requirement  of  the  duplication  of  the 
free-stream  thermodynamic  state  can  be  relaxed.  Such  a  similitude  will  be 
derived  in  the  present  section. 

Now  Cheng's  similitude  (for  the  pressure  field)  for  a  blunt-nosed 
slender  body  in  an  equilibrium  flow  can  be  written  as 


r  l  T 


M  i 


(VI.  4) 


In  the  above  equation  L  is  a  characteristic  body  length,  r  is  the  thickness 
ratio,  Ji  is  the  nose  drag  coefficient,  and  the  free-stream  has  been  assumed 
to  be  undissociated  for  the  sake  of  simplicity.  The  quantity  j)  is  to  be 
taken  as  zero  for  two-dimensional  flows  and  unity  for  axisymmetric  flows. 

In  the  above  similitude  the  requirement  of  duplication  of  the  free- 
stream  temperature  becomes  unnecessary,  by  virtue  of  the  Mach-nurnber- 


*  In  principle, free -stream  dissociation  can  be  accounted  for  without  additional 
difficulties.  For  example,  it  might  be  possible  to  use  the  'Subtraction  Rule' 
concept  developed  by  Gibson v  in  conjunction  with  the  present  formulation. 
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independence  principle,  in  the  limit  when  7  — 1 -  (or  more  appro¬ 
priately,  when  M  T  >>  /  ).  In  this  case,  the  similitude  reduces  to 


Alternately,  when  Mco  t  -  ocn  ,  the  real-gas  effects  become 
negligible  throughout  the  flow  field,  and  the  similitude  becomes 


JL  _  £  [±  JL  .  M  r  U  L  A. 
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(VI.  5) 


(VI.  6) 


In  the  intermediate  region  where  /»f  t'  is  neither  so  large  that 
the  Mach-number  independence  principle  becomes  operable,  nor  so  small 
that  real-gas  effects  can  be  neglected,  the  full  similitude  given  by  Eq.  VI.  4 
has  to  be  used.  However,  it  was  shown  in  Section  III  that  for  a  significant  range 
of  free-stream  conditions,  the  real-gas  effects  are  important  only  in  part  of 
the  flow  field,  and  that  the  remainder  of  the  flow  field  behaves  like  an  ideal 
gas  one.  Thus  in  the  flow  past  a  blunt-nosed  slender  body,  the  displacement 
effect  due  to  the  blunt  nose  is  influenced  by  real-gas  effects  while  the  evolu¬ 
tion  of  the  shock  wave  due  to  the  resulting  effective  slender-body  shape  is  not. 
Heuristically  it  would  appear  that  the  above  observation  should  lead  to  a  sim¬ 
plification  of  the  similitude  given  in  Eq.  VI.  4.  We  will  now  investigate  this 
possibility. 

It  should  be  pointed  out  here  that  many  of  the  restrictions  that  were 
prescribed  for  the  model  formulated  in  the  previous  section  are  not  necessary 
for  the  derivation  of  the  similitude  given  below  (these  restrictions  are  necessary 
only  for  the  development  of  the  zeroth-order  theory  given  in  the  next  section). 
For  example, the  ratio  of  the  specific  heats  need  not  be  close  to  unity,  and  the 
nose  pressure  distribution  need  not  be  of  any  specified  form.  Also,  the  deri¬ 
vation  given  below  is  not  restricted  to  any  specific  model  for  the  flow  chemistry 
within  the  entropy  layer.  In  particular,  the  analysis  will  be  valid  both  for  the 
case  when  the  flow  is  assumed  to  be  in  equilibrium  throughout  the  entropy 
layer,  and  for  the  case  when  the  flow  is  in  equilibrium  in  the  nose  region  and 
freezes  downstream  of  some  point  near  the  nose-afterbody  junction.  The 
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above  two  models  have  been  termed  equilibrium-equilibrium,  and  equilibrium- 
frozen  models  in  Refs.  1  and  2.  The  primary  requirement  for  the  validity  of 
the  similitude  described  below  is  that  real-gas  effects  be  important  only  in  the 
high-entropy  layer  close  to  the  body  and  that  they  be  negligible  outside  this 
layer. 

Before  we  proceed,  one  other  point  needs  to  be  made  here.  Cheng's 
25 

similitude  for  blunt-nosed  slender  bodies  is  based  on  the  postulate  that  the 
downstream  influence  of  a  blunt  nose  is  equivalent  to  that  of  a  concentrated 
force  of  appropriate  magnitude  applied  at  the  nose  in  a  direction  transverse 
to  the  free-stream  direction.  The  magnitude  of  the  force  is  assumed  to  be 
determined  uniquely  by  the  nose  drag  coefficient  and  is  independent  of  the 
detailed  nose  shape.  Therefore,  like  the  blast-wave  analogy  itself,  the  sim¬ 
ilitude  is  necessarily  asymptotic  in  nature  and  becomes  applicable  only  at  a 
considerable  distance  downstream  of  the  blunt  nose.  However,  in  many 
practical  problems  similitude  is  required  not  only  in  the  region  considerably 
downstream  of  the  nose,  but  in  the  nose  region  itself.  In  such  cases  Eq.  VI.  4 
is,  of  course,  inadequate. 

Similitude  in  the  nose  region  can  be  assured  by  invoking  the  Mach- 
number-independence  principle.  Thus,  for  the  nose  region  one  can  write 
(assuming  geometrically  similar  boundaries), 


The  simultaneous  requirement  of  the  conditions  in  Eqs.  VI.  4  and  VI.  7  implies 
that  essentially  no  similitude  exists;  that  is,  the  quantities  /o  ,  ~f  ,  u 
and  7  have  all  to  be  duplicated  simultaneously.  The  only  freedom  of  choice 
that  exists  is  in  connection  with  the  length  scale  L  ,  and  even  this  freedom 
is  lost  if  nonequilibrium  effects  or  viscous  effects  are  important.  For  such 
cases,  a  similitude  which  relaxes  the  requirement  of  the  duplication  of  the 
free-stream  temperature  would  be  of  great  interest. 

Consider  two  blunt-nosed  slender  bodies  (shown  below)  which  have  geo¬ 
metrically  similar  nose  and  afterbody  shapes,  but  are  characterized 
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by  different  nose  radii,  thickness  parameters,  characteristic  lengths  and  are 
immersed  in  different  hypersonic  flows.  Now,  as  mentioned  earlier,  the 
downstream  effect  of  a  blunt  nose  on  the  afterbody  flow  field  appears  essentially 
as  a  displacement  effect.  This  displacement  is,  of  course,  dependent  on  the 
initial  conditions  at  the  nose  as  well  as  on  the  afterbody  shape.  The  conditions 
at  the  nose  can  be  made  identical  for  the  two  cases  by  requiring  that  the  free- 
stream  density  and  velocity  be  the  same  for  both.  Then,  by  virtue  of  the  Mach- 
number- independence  principle,  all  quantities  in  the  nose  region  of  the  two 
bodies  will  be  not  merely  similar  but  identical.  In  particular,  the  nose  drag 
(and  therefore  the  energy  and  impulse  imparted  to  the  transverse  flow  field) 
will  scale  in  a  simple  geometric  way. 

Note  that  the  above  statements  hold  true  whether  or  not  the  flow  in  the 
region  of  the  nose  is  in  chemical  equilibrium,  and  whether  or  not  the  flow  is 
assumed  to  freeze  at  some  point  in  the  neighborhood  of  the  sonic  point.  Irre¬ 
spective  of  the  exact  nature  of  the  flow  chemistry,  the  displacement  effect 
(which  will  be  assumed  to  be  given  by  Ta  )  due  to  the  blunt  nose  must  be 
governed  solely  by  the  initial  conditions  at  the  nose  and  the  afterbody  thick¬ 
ness  parameter  T  .  Thus  one  can  write 

r.  s  O  •  (VI.  8) 


Moreover,  under  the  assumptions  of  hypersonic  small  disturbance  theory 
Eq,  VI.  8  can  be  written  as 


(VI.  9) 


62 


AEDC-TR-69-36 


The  flow  field  due  to  the  effective  body  shape  governed  by  the  paramete 
Te  can  be  expressed  in  terms  of  the  free-stream  <  onditons  using  the  ideal- 
gas  relations  alone,  since  by  assumption  real-gas  eftects  are  negligible  out¬ 
side  the  entropy  layer  in  the  present  analysis.  In  particular,  the  pressure 

25 

distribution  can  be  expressed  in  the  form 
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Eqs.  VI.  9  and  VI.  10  can  be  combined  to  yield 
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A  practically  more  useful  form  of  Eq.  VI.  1 1  is  given  by 


(VI.  10) 
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(VI.  12) 


Similar  relations  can  be  written  down  for  the  shock  shape  as  well  as  the  other 
flow  quantities. 


Eq.  VI.  12  can  have  signficant  applications  in  wind-tunnel  testing, 
since  it  affords  a  means  of  taking  into  account  mismatches  in  free-stream 
temperature  in  the  flows  over  blunt-nosed  slender  bodies.  It  was  shown  in 
Ref.  1  that  for  certain  flight  conditions  the  free-stream  temperature  in  a  wind 
tunnel  could  be  as  high  as  a  factor  of  five  larger  than  the  corresponding  flight 
value,  when  the  density-altitude  and  velocity  are  duplicated  in  the  tunnel.  For 
such  cases,  Eq.  VI.  12  provides  a  means  of  obtaining  proper  similitude 
between  the  flight  and  tunnel  flows  by  choosing  the  thickness  ratio  of  the 
afterbody  appropriately.  As  mentioned  earlier,  the  range  of  free-stream 
conditions  over  which  the  assumption  that  real-gas  effects  are  important  only 
in  the  entropy  layer  (and  hence  the  range  of  validity  of  Eq.  VI.  12),  is  sig¬ 
nificant  and  certainly  nontrivial. 


The  similitude  expressed  by  Eq.  VI.  12  can  also  be  derived  in  a  more 
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rigorous  manner  from  the  basic  differential  equations  governing  the  flow.  For 
the  sake  of  completeness  we  will  indicate  this  derivation  below. 


Under  the  small-disturbance  approximations  that 


and 
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(VI.  13) 


(where  u.  and  or  are  the  perturbation  velocity  components  and  0  is  the  total 
axial  velocity)  the  continuity,  momentum  and  entropy  equations  can  be  written 
in  the  form 
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(VI.  14) 
(VI.  15)  . 
(VI.  16) 

(VI.  1 7) 


Note  that  the  form  of  the  first  three  of  the  above  equations  is  the  same  whether 
or  not  any  real-gas  effects  are  present. 


The  boundary  condition  to  be  satisfied  by  the.  above  equations  at  the 
body  is  simply  that 


(VI.  18) 


Since  in  the  present  model  real-gas  effects  are  assumed  to  be  negligible 
within  the  shock  layer,  the  boundary  conditions  on  the  quantities  or  ,  and 
P  at  the  oblique  shock  wave  ^  -  Ys  (*)  can  be  written  as 
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and 
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(VI.  19) 

(VI.  20) 

(VI.  21) 
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Cheng  has  shown  that  the  downstream  influence  of  the  blunt  nose  can  be 
accounted  for  within  the  framework  of  small-disturbance  theory  by  using  the 
equation 
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as  the  initial  condition  in  the  nose  region,  that  is, as  4. — >0 
equation.  represents  the  total  inviscid  drag  of  the  nose. 


(VI.  22) 


In  the  above 


The  above  equations  are  the  same  as  those  given  by  Cheng  except 
for  the  simplified  form  of  the  boundary  condition  equations  at  the  shock  wave 
in  the  prt-,ent  case  (Eqs.  VI.  19  to  VI.  21).  The  system  of  equations  given 
by  Eqs.  VI.  14  to  VI.  17  and  the  appropriate  boundary  and  initial  conditions 
will  constitute  a  closed  system  if  the  functional  dependence  of  the  entropy  on 
the  other  state  variables  is  prescribed. 


Within  the  shock  layer,  the  entropy  can  be  expressed  in  terms  of  the 
local  values  of  pressure  and  density  using  the  perfect  gas,  isentropic  relation 
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(VI.  23) 


where  the  value  of  the  constant  is  unimportant  since  only  changes  in  entropy 
from  some  reference  state  are  of  interest  here.  Thus  within  the  shock  layer, 
the  entropy  Eq.  VI.  17  can  be  written  as 


(VI.  24) 
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Note  that  Eq.  VI.  24  is  valid  only  for  4.  y  4  Yg  • 

For  the  range  given  by  Yfc  <  y  <  ,  a  more  general  form  of 

the  entropy  equation  is  necessary.  Now,  the  entropy  at  any  point  in  the  entropy 
layer  is  primarily  dependent  on  the  conditions  in  the  nose  region.  In  particular, 
because  of  the  particle  isentropic  condition,  the  entropy  at  any  point  is  equal 
to  the  value  corresponding  to  that  on  the  streamline  passing  through  that  point, 
and  this  latter  quantity  is  given  by  its  value  immediately  behind  the  shock  wave. 
Thus  the  entropy  at  any  point  can  be  determined  if  the  quantity  1  y#'  corre¬ 
sponding  to  that  point  (see  sketch  below)  can  be  expressed  in  terms  of  other 


local  variables  and  the  free-stream  conditions.  The  desired  type  of  relation 
is  provided  by  the  mass  conservation  relation,  which  is 
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(VI.  25) 


where  U  is  the  total  axial  velocity  component. 

The  entropy  within  the  entropy  layer  can  now  be  expressed  in  the  form 
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(VI.  26) 


Eq.  VI.  26  is,  of  course,  a  direct  consequence  of  the  Mach-number-indepen- 
dence  principle. 


For  the  purposes  of  studying  the  similitude  expressed  by  the  above 
equations,  the  transformation 
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is  made. 

Under  this 

transformation  Eqs. 

VI.  14  to  VI.  16  become 

\ 

Bfi 

3* 

-  -h 

4^ 

J£_. 

~ 

-  /° 

9  -y. 

*  j  rl 
j  J 

(VI.  28) 

iw.x  f 

3u 

4- 

/lA 

9£ 

1 

9 

(VI.  29) 

i 

3* 

T 

__ 

P 

'  dZ 

_  _i- 

*✓* 

d^r 

/ 

(VI.  30) 

7 

df\ 

£ 

i 

The  boundary 

conditions 

are 

4T 

on 

/•  / 

=  rfc  (z> 

(VI.  31) 

t 

cCf 

and 


'V 


”(''?)  OLt 

’  *  0  ‘  T')  (t!)*  <r-')  <"„rf 


p 


(VI.  32) 

(VI.  33) 

(VI.  34) 


on  y  «  y  (?) 


The  initial  condition  at  the  nose  given  by  Eq.  VI.  22  becomes 
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where  the  specific  internal  energy  ' e '  has  been  nondimensionalized  with 

/  ’fi' 

respect  to  fhe  free-stream  value,  given  by 

For  the  shock  layer  (that  is, for  Yg  4  -y  4  Ys  )  the  entropy 
equation  becomes 
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For  the  entropy  layer  (that  it,  Y  4  <  Y  )  the  entropy  equation  is 
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(VL  38) 


la  Zq.  VL  3S.  the  total  axial  velocity  has  been  noncimensionalized  by  the 
tree-stream  velocity  .  so  that  U  =  .  (Note  that  under  the  small- 

disturbance  assumptions  U  ~  and  u  ~  i ) - 


In  principle,  Eos.  VL  28  to  VL  38  are  sufficient  to  solve  for  the  entire 
flow  field.  The  parameters  appearing  in  these  equations  are 
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After  some  manipulation  these  parameters  can  be  reduced  to 
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where  the  nose  drag  coefficient  *  4  *  is  defined  as 
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i>  .ce  the  present  analysis  is  confined  to  bodies  with  geometrically  similar 
nose  shapes,  the  nose  crag  coefficient,  4  .  need  not  be  included  in  the  si~iil- 
itude  criteria. 


The  similitude  relation  for  the  pressure  distribution  can  be  written  in 
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with  similar  relations  being  applicable  for  the  other  flow  variables. 

The  similitude  derived  above  for  a  general  blunted  slender  body  can 
be  specialized  for  blunted  wedges  and  cones.  For  these  geometries,  the 
equation  of  the  surface  may  be  written  as 


(VI.  43) 


where  3.  is  the  ha  If -angle  of  the  wedge  or  cone  under  consideration  and  plays 
the  role  of  7  m  the  more  general  case.  Since  for  the  wedge  or  cone  case 
a  characteristic  body  length.  L  ,  does  not  exist,  we  can  eliminate  this  ouan- 
tity  from  the  similitude  by  choosing  it  to  be  equal  to  d  .  men 

Eq.  VL  42  becomes 
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Similarly,  the  shock  shape  can  be  expressed  in  the  form 
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Note  that  in  Eqs.  VI.  42,  VI.  44  and  VI.  45,  the  various  combinations 
of  the  parameters  of  the  problem  that  have  been  used  are  not  nondimensional. 
These  equations  have  been  written  in  the  form  given  mainly  to  illustrate  in  a 
clear  manner  the  effects  of  a  mismatch  in  free-stream  temperature.  As 
mentioned  earlier,  the  effects  of  mismatches  in  free-stream  temperature 
are  of  primary  concern  in  hypervelocity  testing.  However,  if  desired,the 
above  equations  can  be  easily  rewritten  in  a  more  conventional,  dimensionally 
consistent  form  by  making  use  of  the  fact  that  products  and  ratios  of  similitude 
parameters  are  also  similitude  parameters.  Thus,  Eq.  VI.  44  can  also  be 
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In  order  to  verify  the  similitude  given  above,  some  numerical  solutions 

have  been  carried  out  using  the  computer  program  developed  by  Curtis  and 
40 

Strom  .  The  calculations  were  performed  for  a  spherically  blunted  cone 
immersed  in  an  equilibrium  hypersonic  flow  of  25,  000  ft/sec  velocity  and  a 
free-stream  density  corresponding  to  250,  000  ft  altitude.  Two  cases  were 
considered-  In  the  first  case,  tl--e  free-stream  temperature  chosen  was  that 
corresponding  to  the  atmospheric  value  at  an  altitude  of  250,  000  ft,  and  the 
cone  half-angle  for  this  case  was  9*  -  In  the  second  case,  the  free-stream 
temperature  was  chosen  to  be  a  factor  of  four  higher  than  the  value  corre¬ 
sponding  to  an  altitude  of  250,  000  ft  {as  mentioned  earlier,  for  the  flight 
conditions  that  are  being  considered  here  mismatches  in  the  free-stream 
temperature  of  this  order  are  likely  to  occur  in  most  hypersonic  test  facil¬ 
ities).  The  similitude  derived  above  indicates  that  for  proper  similitude 
between  the  two  cases,  the  cone  half-angle  for  the  second  case  has  to  be 
chosen  as  IS*. 

The  calculated  results  are  shown  in  Figs.  25  and  26  correlated  in 

$ 

terms  of  the  similitude  derived  above.  Figure  25  illustrates  the  shock  shape, 
while  Fig.  26  illustrates  the  calcula  .ed  surface  pressure  distribution.  It  can  be 
seen  that  in  spite  of  the  fairly  large  mismatches  in  the  free-stream  temper¬ 
ature  the  correlations  obtained  ar  *  fairly  good. 

E.  ANALYSIS  OF  BLUNT-NOSED  SLENDER  BODY 

In  the  present  section  we  will  analyze  the  hypersonic  flow  past  a  blunt- 
nosed  slender  body  using  the  approximate  physical -chemical  model  developed 
in  Section  C.  The  essential  features  of  this  model  are  that  the  hypersonic 
small-disturbance  approximations  are  made  throughout  the  flow  field  except  in 
a  region  close  to  the  blunt  nose,  and  that  the  flow  field  is  divided  roughly  into 
two  regions  (a  high  entropy  region  close  to  the  body  where  real-gas  effects 


*  The  authors  express  their  thanks  to  Dr.  J.  Curtis  of  the  Aerodynamic 
Research  Department  for  performing  these  calculations. 
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are  included  and  a  relatively  low  entropy  region  close  to  the  shock  wave  where 
the  gas  behavior  is  assumed  to  be  ideal). 


1.  The  Basic  Eauations.  The  method  of  analysis  used  in  the  present 

- : -  ^  74,75 

section  (which  closely  follows  the  methods  of  Cheng  and  Chernyi  *  )  is 

essentially  integral  in  nature;  that  is,  the  partial  differential  equations  gov¬ 
erning  the  problem  are  reduced  to  a  single  ordinary  differential  equation  by 
making  various  assumptions  regarding  the  transverse  distributions  of  various 
quantities  within  the  disturbed  region  of  the  flow.  Here,  we  shall  start  with 
a  consideration  of  the  energy  conservation  integral  for  the  transverse  flow 

field  including  the  contributions  due  to  the  blunt  leading  edge.  The  energy 

54 

integral  for  the  present  problem  may  be  written  as 
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where  e  is  the  specific  internal  energy  of  the  gas.  Eq.  VI.  47  expresses 
the  conservation  of  energy  for  an  observer  who  is  fixed  in  the  undisturbed 
fluid. 


Ir:  Eq.  VI.  47  the  specific  internal  energy  e  is  given  by 
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for  ye  <  y  <  Ys 


,  and 


(VI.  49) 


for  <  y-  <  Ye  ,  where  is  the  energy  contained  in  the  dissoci¬ 

ation  mode  (within  the  entropy  layer)  and  j'.  is  the  frozen-flow  specific  heat 
ratio  corresponding  to  the  local  value  of  the  dissociation  fraction. 

If  in  Eq.  VI.  47  the  small-disturbance  approximation  U  ~  uv  is 
made,  one  obtains 
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(ff  Yj*’  d*  =  J  *  f(e+j*'a)  {jT'f)  dy 


K 

r-i 


Y 


frY*)* 

l+>> 


(VI.  50) 


The  second  term  on  the  left-hand  side  of  Eq.  VI.  50  represents  the  contribution 
of  the  afterbody  drag  to  the  overall  momentum  balance.  In  the  large  Mach 
number  regime,  the  second  term  in  the  right-hand  side  of  Eq.  VI.  50  may  be 
neglected.  Then  Eq.  VI.  50  can  be  written  in  the  form 


Zrr’’ 


+ 


Yk  d* 


(VI.  51) 


This  relatively  simple  expression  will  form  the  basis  for  much  of  the  following 
analysis. 


Before  we  proceed,  a  point  of  clarification  is  in  order  here.  As  men¬ 
tioned  earlier,  the  analyses  of  Cheng  and  Chemyi  *’  have  considerable 
similarities,  and  both  analyses  use  an  equation  of  the  form  given  by  Eq.  VI.  51 
However,  in  his  analysis,  Cheng  neglected  the  afterbody  and  kinetic 
energy  integrals  since  these  are  of  higher  order  (in  the  parameter 

+  than  the  remaining  terms.  On  the  other  hand,  Chernyi  included 

these  terms  in  his  analysis.  The  pressure  distributions  calculated  by  Cheng 

for  blunted  wedges  and  cones  show  an  oscillatory  behavior  far  downstream 

7a  75  31 

of  the  nose,  whereas  the  ones  calculated  by  Chernyi  *’  do  not.  Cheng 

notes  that  while  this  oscillatory  behavior  is  not  unbelievable,  it  is  of  doubtful 

54 

physical  reality.  Hayes  and  Probstein  *  point  out  that  the  oscillatory  behavior 
in  Cheng's  results  may  be  caused  by  the  neglect  of  the  kinetic -energy  terms. 


31 


In  the  present  analysis,  following  Cheng,  we  will  neglect  the  afterbody 
and  kinetic -energy  integrals  since  this  is  consistent  with  small  disturbance 
approximations  and  with  the  other  assumptions  involved  in  the  present  model. 
However,  later  we  will  briefly  indicate  the  application  of  Cnernyi's  method  to 
the  present  problem. 


When  the  afterbody  and  kinetic -energy  integrals  are  neglected  in  Eq.  VI.  51, 
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and  when  the  resulting  expression  is  combined  with  Eqs.  VI.  48  and  VI.  49, 
one  obtains 
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(VI.  52) 


Since  the  pressure  in  the  entropy  layer  is  of  the  same  order  of  magnitude  as 
that  in  the  shock  layer  and  since  the  lateral  extent  of  the  shock  layer  is  small 
compared  to  that  of  the  entropy  layer  {that  is.  Yc  -  Y  «  Yc  -  Yfa  ).  the 
second  term  in  the  right-hand  side  of  Eq.  VI.  52  can  be  neglected  compared 
to  the  first.  Moreover,  as  shown  earlier,  the  pressure  inside  the  entropy 
layer  is  equal  (at  each  axial  location)  to  that  at  its  outer  edge.  Therefore 
Eq.  VI.  52  reduces  to 


2ir' 


■t ' 


(VI.  53) 


The  second  term  in  the  left-hand  side  of  Eq.  VI.  53  merely  represents 
the  total  energy  that  is  frozen  out  of  the  flow  within  the  entropy  layer,  and  by 
the  assumptions  of  the  present  model  (that  is,  ’freezing’  of  the  flow  along  a 
line  passing  through  the  sonic  point)  this  quantity  must  be  independent  of  the 
distance  downstream  of  the  nose.  Thus,  as  noted  by  Whalen3  ,  one  of  the 
principle  influences  of  real -gas  effects  on  blunt-nosed  slender-body  flows  is 
to  cause  an  effective  reduction  in  the  nose  bluntness  and  hence  also  its  down¬ 
stream  influence.  The  net  energy  introduced  into  the  transverse  flow  field 
by  the  nose  bluntness  is  determined  by  the  effective  drag  ,  defined  by 
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(VI.  54) 


Note  that  this  quantity  is  always  positive  since  the  right-hand  side  of  Eq.  VI.  53 
is  always  positive. 


A  net  drag  coefficient  &  can  he  defined  by  the  relation 


Jk  =  J). 
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(VI.  55) 
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However,  note  that  tinlike  in  the  equivalent  ideal-gas  case,  is  not  inde¬ 

pendent  cf  the  free- stream  state  and  velocity  since  the  amount  of  energy  that 
is  frozen  out  at  the  nose  as  chemical  energy  is  a  function  of  the  free-stream 
conditions. 

The  pressure,  ,  at  the  outer  edge  of  the  entropy  layer  can  be 

related  to  the  free-stream  conditions  and  the  shape  of  the  outer  edge  of  the 
entropy  layer  through  the  Newtonian-plus -centrifugal  pressure  formula 

■p  (*)  -  U  2  (y'Z  +  “%-)  (VI.  56) 

e  '  a>  CO  V  e  /*'*'/ 

where  the  primes  denote  differentiation  with  respect  to  the  axial  coordinate, 

2. 


Finally.  Eqs.  VI.  54-56  can  be  combined  with  Eq.  VI.  53  to  give 


Eq.  VI.  57  is  the  fundamental  equation  which  determines  the  rate  of  growth 
of  the  entropy  layer  and  thus,  also, the  displacement  effect  of  the  blunt  nose  on 
the  afterbody  flow  field  including  the  effects  of  flow  nonequilibrium. 

As  mentioned  earlier  (and  as  can  be  seen  from  Eq.  VI.  57), one  of  the 
principle  effects  of  flow  nonequilibrium  on  the  afterbody  flow  field  appears 
through  a  net  reduction  in  the  energy  release  at  the  nose  to  the  transverse 
flow  field.  Moreover,  the  effective  drag  coefficient  Ji!  (unlike  the  equivalent 
ideal-gas  drag  coefficient  Jk, )  is  a  function  of  the  free-stream  conditions, 
and  in  general  it  decreases  with  increasing  free-stream  velocity.  That  is, 
at  high  velocities  the  amount  ox  energy  that  is  frozen  out  as  chemical  energy 
by  the  rapid  expansion  of  the  flow  in  the  shoulder  region  is  a  significant 
fraction  of  the  energy  that  is  introduced  into  the  transverse  flow  field  by  the 
blunting.  Thus  when  real-gas  effects  are  included, the  downstream  influence 
of  a  blunt  nose  can  be  significantly  smaller  than  that  calculated  using  ideal- 
gas  theories. 
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Another  important  consequence  of  real-gas  effects  on  the  downstream 
influence  of  a  blunt  nose  appears  in  Eq.  VI.  57  through  the  integral  on  the 
lef'.-hand  side.  It  can  be  seen  from  this  integral  that  the  local  rate  of  growth 
of  the  entropy  layer  is  influenced  by  the  distribution  of  the  frozen-flow  specific 
heat  ratio,  7^  ,  within  the  entropy  layer.  The  local  value  of  the  frozen-flow 
specific  heat  ratio  (which  is  related  to  the  local  value  of  the  dissociation 
fraction,  oc  )  is,  of  course,  a  function  of  the  details  of  the  flow  field  in  the 
nose  region.  Thus  when  real-gas  effects  are  present  the  rate  of  evolution 
of  the  entropy  layer  and  the  afterbody  pressure  distribution  cannot  be  deter¬ 
mined  independently  of  the  details  of  the  nose  shape. 

In  ideal-gas  theories  of  blunted  slender  body  flows,  it  is  assumed  that 

the  effects  of  the  blunt  nose  on  the  afterbody  flow  field  can  be  represented  by 

an  appropriate  concentrated  force  applied  at  the  nose;  the  details  of  the  nose 

shape  are  assumed  to  be  of  little  importance.  Even  though  this  so-called 

'blast-wave  analogy1  leads  to  erroneous  results  for  the  density  and  entropy 

distributions  near  the  afterbody;  it  is  known  that  the  afterbody  pressure  field 

and  the  shock  wave  shape  are  predicted  correctly.  However,  j*  can  be  seen 

from  Eq.  VI.  57  that  when  real-gas  effects  are  present  even  the  afterbody 

pressure  field  and  shock  shape  cannot  be  solved  for  independently  of  the  details 

of  the  flow  in  the  nose  region.  This  point  has  been  mentioned  earlier  by 
31 

Cheng 

Eq.  VI.  57  is  also  consistent  with  the  similitude  derive'4  in  the  previous 
section  according  to  which  similitude  between  two  blunted  slender  body  flow 
fields  can  be  obtained  only  if  the  flows  in  the  nose  region  of  the  two  cases  are 
also  similar.  Thus  even  if  it  were  possible  for  two  slender  bodies  with  dif¬ 
ferent  nose  shapes  to  have  the  same  values  of  the  effective  drag  coefficient, 

Jfd  ,  over  at  least  a  small  range  of  free-stream  conditions,  it  is  highly  unlikely 
that  the  two  bodies  will  also  have  similar  distributions  of  the  frozen-flow 
specific  heat  ratios  within  their  entropy  layers.  Therefore, 'blast-wave 
analogy'in  the  usual  sense  (that  is,  in  the  sense  that  the  downstream  influence 
of  a  blunt  nose  can  be  considered  to  be  independent  of  the  detailed  nose  shape) 
is  not  applicable  to  the  present  case. 
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2.  Evaluation  of  the  Integral  in  Eq.  VI.  57.  Before  quantitative  answers 
can  be  obtained  from  Eq.  VI.  57,  the  integral  appearing  in  its  left-hand  side 
has  to  be  evaluated  explicitly.  However,  as  mentioned  above,  this  task  can¬ 
not  be  accomplished  without  knowing  the  details  of  the  flow  in  the  nose  region. 
On  the  other  hand,  consistent  with  the  integral  nature  of  the  analysis  one  can 
assume  various  profile  shapes  for  the  distribution  of  the  effective  frozen-flow 
specific  heat  ratio  within  the  entropy  layer.  Any  arbitrary  constants  used  in 
the  profile-shape  equations  can  be  evaluated  by  using  the  appropriate  boundary 
conditions  and  (or)  various  integral  conservation  equations. 

The  simplest  assumption  one  can  make  to  evaluate  the  integral  in 
Eq.  VI.  57  is  to  assume  that  the  entropy  layer  can  be  characterized  by  an 
average  frozen-flow  specific  heat  ratio,  .  Then,  the  integral  can  be  eval¬ 
uated  immediately,  and  Eq.  VI.  57  reduces  to 


(VI.  58) 


where  T  _  rt  ~ 1  _.  ~ 1  (VI.  59) 

t  -  - 5 —  £2  — = — — 

2  rf  +1 
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Equation  VI.  58  is  identical  to  the  equation  given  by  Cheng  for  the  xdeal-gas 

case  except  that  in  the  present  case  both  6  and  A!  are  functions  of  the  free- 

stream  conditions.  Moreover,  note  that  in  the  present  case  6  may  be  a 

function  of  .  In  the  present  approximation  each  axial  station  in  the  entropy 

layer  is  assumed  to  be  characterized  by  an  average  frozen-flow  specific  heat 

ratio;  however,  this  average  specific  heat  ratio  may  vary  from  station  to 
* 

station. 


Since  the  value  of  the  frozen  flow  specific  heat  ratio  depends  only  on 
the  local  value  of  the  dissociation  fraction,  c<  ,  the  value  of  the  average 
specific  heat  ratio  at  any  station  can  be  determined  by  matching  for  the  two 
cases  the  value  of  the  total  energy  contained  in  the  dissociation  mode  at  that 


*  The  shock  layer  is,  of  course,  assumed  to  be  characterized  by  a  constant 
(ideal  gas)  value. of  the  specific  heat  ratio. 


76 


AEDC-TP.-69.36 


station.  That  is,  the  value  of  the  average  specific  heat  ratio  is  assumed 

to  be  such  that  the  dissociation  energy  it  represents  is  equal  to  that  corre¬ 
sponding  to  the  actual  distributions  of  the  flow  quantities  at  the  station  under 
consideration.  Thus,  the  condition  determining  f ^  is  represented  by 


P<*- 


(VI.  60) 


The  left-hand  side  of  Eq.  VI.  60  is,  of  course,  independent  of  the  axial  position 
and  is  equal  to  the  value  at  the  nose. 

Now,  within  the  accuracy  of  the  small  disturbance  approximation 
U  ~  ,  the  integral  conservation  of  mass  is  governed  by  the  equation 

/>„  v  (* if  =  V  p  fcrf  dr  ■  <VL  61) 

Thus  Eqs.  VI.  60  and  VI.  61  can  be  combined  together  to  yield 


\  5  dt '  y.  (p'yJ  *  •  (vi.  62) 

To  be  consistent  with  our  previous  approximations,  the  approximation  Ys  ~ 
has  been  made  in  Eq.  VI.  62.  The  average  specific  heat  ratio,  ,  is  related 
to  the  average  dissociation  fraction,  oC  ,  through  the  relation 


3Z  i-  7 
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(VI.  63) 


Equations  VI.  58,  59,  62  and  63  determine  the  shape  of  the  entropy 
layer  (for  a  given  body  shape,  )  in  terms  of  the  integral  in  the  left-hand 
side  of  Eq.  VI.  62.  Note  that  the  above  integral  also  determines  the  value  of 
the  effective  drag  coefficient,  Jt!  (see  Eq.  VI.  54).  Thus  in  principle,  given 
sufficient  details  near  the  nose  for  the  evaluation  of  the  above  integral,  one 
can  solve  for  the  shock  shape  and  the  body-surface  pressure  distribution.  The 
evaluation  of  this  integral  will  be  discussed  later. 
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It  is  appropriate  here  to  consider  other  possible  approximations  that 
can  be  used  in  conjunction  with  Eq.  VI.  57.  The  simplest  of  these  is  to  assume 
that  the  whole  entropy  layer  can  be  characterized  by  a  single  average  value 
(independent  of  itf)  of  the  specific  heat  ratio  which  is  given  by  the  arithmetic 
mean  of  ;he  values  at  the  inner  and  outer  edges  of  the  entropy  layer.  Thus, 

_  +  * 

r(  =  - £ -  (VI.  64a) 


or  alternately, 


(VI.  64b) 


Of  course,  apart  frpm  its  very  attractive  simplicity,  the  assumption  embodied 
in  Eqs.  VI.  64a  and  b  does  not  really  have  a  sound  rational  basis. 


Another  possible  assumption  is  to  take  the  transverse  distribution  of 
the  effective  froifen-flow  specific  heat  vauo  within  the  entropy  layer  to  be  a 
linear  function  of  the  cross-sectional  area  of  the  flow  field.  Thus,  one  can 


write 


»)♦•  >>+ 1 
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(VI.  65) 


The  choice  of  the  cross-sectional  area  of  the  flow  field  as  the  appropriate 

parameter  is  prompted  by  the  fact  that  in  all  blunted  slender-body  flow  analyses 

the  lateral  coordinate  always  enters  in  the  form  -2^  or  its  derivative  (this 

quantity  is  of  course  proportional  to  the  cross-sectional  area  of  the  flow  field). 

..94 

The  above  property  has  been  exploited  by  Lady  zhenskii  who  derived  an 
area  rule  for  blunted,  slender  three-dimensional  bodies  which  have  the  same 
axial  distribution  of  body  cross-sectional  area. 

Note  that  Eq.  VI.  65  satisfies  the  appropriate  boundary  conditions  at 
the  inner  and  outer  edge3  of  the  entropy  layer.  Equation  VI.  65  can  be  regarded 
as  being  composed  of  the  zeroth  and  first-order  terms  in  the  small  parameter 
(  — 7 ..  )  and,  in  principle,  more  terms  can  be  included  if  desired.  However, 

'7  .  y  _  y  . 

lor  air, the  parameter  (  — Lp —  j  has  a  possible  maximum  value  of  only  about 
0.  19.  Thus,  the  errors  involved  in  the  terms  neglected  in  Eq.  VI.  65  would 

—t— — J  ,  or  of  the  order  of  four  percent. 
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Therefore  inclusion  of  additional  terms  in  Eq.  VI.  65  does  not  appear  to  be 
worthwhile. 

When  Eq.  VI.  65  is  introduced  in  Eq.  VI.  57,  one  obtains 


vf)  - 


;vi.  66) 


Equation  VI.  66  represents  the  'pressure -area1  relation  governing  th  ;  entropy 
layer  of  any  arbitrary  plane  or  axisymmetric  blunted  slender  body. 


Equation  VI.  66  reduces  to  the  corresponding  ideal-gas  relation  given 
3? 

by  Cheng  '  in  the  limit  ih — »/  (or  more  appropriately  ■/  -  y  <£<.  y-  /) 
since 
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(VI.  67) 


Several  approximate  methods  have  been  discussed  above  to  simplify 
the  basic  equation  (Eq.  VI.  57)  governing  the  entropy  layer.  All  these  methods 
were  based  on  simplifying  assumptions  on  the  nature  of  the  tra.nsverse  distri¬ 
bution  of  the  various  flow  quantities  within  the  entropy  layer.  Of  the  three 
methods  described,  those  represented  by  Eqs.  VI.  62  and  VI.  66  are  likely  to 
be  more  accurate  than  that  represented  by  Eq.  VI.  64  since  the  former  are 
based  on  rational  simplifications  and  satisfy  appropriate  boundary  and  con¬ 
servation  conditions.  Other  approximations  are  possible,  of  course,  but  we 
shall  confine  our  attention  to  the  three  given  above  and  discuss  them  in  more 
detail  for  the  specific  case  of  a  spherically-blunted  slender  cone.  In  particular, 
we  will  discuss  Eq.  VI.  66  in  some  detail  since  it  clearly  illustrates  the  various 
effects  of  real-gas  phenomena. 

Eq.  Vi.  66  can  be  rewritten  in  the  more  convenient  form  with  i)  -  1 
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(VI.  68) 
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me  above  equation  is  identical  to  the  corresponding  ideal-gas  equation  except 
for  the  terms  occurring  cn  the  right-hand  side  within  the  parenthesis.  L«ei 


(VL  69) 


Obviously,  CZ  represents  a  correction  factor  which  accounts  for  the  real-gas 
effects.  It  can  be  seen  from  Eq.  VL  68  that  influence  of  real-gas  effects  on 
the  flow  field  can  be  visualized  as  being  that  due  to  a  change  in  the  nose  radius 
from  d  to  an  effective  value  of  c  cL. 


The  value  of  the  correction  factor  is  governed  by  two  different  effects. 
The  first  effect,  which  is  represented  bv  the  factor  X  /X  ,  arises  because 
only  part  of  the  energy  introduced  into  the  transverse  flow  field  by  the  blunt 
nose  is  felt  by  the  afterbody  flow  field,  the  rest  being  frozen  our  of  the  flow 
as  dissociation  energy.  Since  the  quantity  -k'/X  is  always  less  than  unity,  the 
influence  of  this  effect  is  to  cause  an  apparent  reduction  of  the  nose  diameter 
as  fa.  as  the  afterbody  flew  field  is  concerned.  The  second  effect,  which  is 
represented  by  the  parenthesized  term  in  Eq.  VI.  69,  arises  from  the  fact 
that  the  effective  frozen  flow  specific  heat  ratio  within  the  entropy  layer  is 
higher  than  the  corresponding  ideal-gas  value.  This  factor  can  be  rewritten 
in  the  form 


c 

v777 ■> 


t  - 


Jizl. 

r- ! 


I VI.  70) 


For  air,  the  quantity  t  has  a  maximum  possible  value  of  0.667.  Therefore, 
the  second  factor  in  Eq.  VI.  69  is  always  greater  than  unity,  its  maximum 
value  for  air  being  around  1.  3. 

Thus  it  appears  that  the  effect  of  the  increased  specific  heat  ratio  within 
the  entropy  layer  is  to  cause  an  increase  in  the  downstream  influence  of  the 
blunt  nose,  the  effect  being  manifested  as  an  apparent  increase  in  the  nose 
diameter  of  up  to  14%.  Hence  the  two  effects  appearing  m  Eq.  VI.  69  display 
opposing  trends  in  their  influence  on  the  slcr.de i -oody  flow  field.  However, 
the  first  effect  is  usually  the  dominant  one,  and  the  influence  of  real-gas 
effects  on  slender-bodv  flows  is  in  general  to  reduce  the  downstream  influence 
of  the  blunt  nose. 
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3.  Solution  for  a  Spherically  Blunted  Cone.  For  a  cone,  the  body  surface 
is  represented  by  the  equation  Y  =  9  -z  where  &c  i3  the  cone  half-angle. 
For  this  case,  Eq.  VI.  63  can  be  transformed  to  the  form 
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{VL  72} 


(VL  73) 
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The  solution  of  the  ordinary  differentional  Eq.  VI.  71  has  been  given  by  Cheng' 

It  can  be  seen  clearly  from  Eos.  VI.  71-73  that  the  influence  of  real-gas  effects 
(contained  in  the  parameter  c  }  is  to  reduce  the  length  scale  for  the  downstream 
influence,  this  length  scale  being  represented  by 
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(VL  74) 


In  order  to  calculate  the  real-gas  correction  factor,  c  ,  the  frozen- 
dissociation-energy  integral  mentioned  earlier  (see  Eqs.  VI.  54  and  VI.  62) 
has  to  be  evaluated.  The  value  of  this  integral,  of  course,  depends  on  the 
details  of  the  flow  field  in  the  nose  region.  Now,  if  a  numerical  solution  were 
available  for  the  nose  region  (that  is  for  the  flow  past  a  sphere)  for  the  free- 
stream  conditions  of  interest,  then  c  can  be  computed  from  this  and  the  solution 
can  be  continued  over  the  afterbody  with  the  help  of  Eq.  VI.  71.  However,  if  a 
numerical  solution  is  not  available,  then  additional  assumptions  are  necessary 
to  evaluate  C  . 

Firstly,  in  the  nose  region  the  shock  wave  will  be  assumed  to  be  parallel 
to  the  body  surface.  Numerical  calculations  show  that  this  assumption  is 
fairly  accurate  for  spherically  shaped  bodies.  This  assumption  is  sufficiently 
accurate  for  the  present  purpose  since  it  is  consistent  with  the  Newtonian- 
flow  assumption,  4  44/  ,  used  in  other  simplifications.  Secondly,  numerical 
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calculations  also  show  that  the  pressure  distribution  on  a  blunt  body  is 
influenced  very  little  by  real-gas  effects.  Therefore,  the  drag  coefficient 
of  the  spherical  nose  will  be  taken  to  be  equal  to  0.  92. 

Let  us  consider  the  spherical  nose  shown  below.  Let  the  shock 
stand-off  distance  be  denoted  by  the  symbol  A  .  Then  since  we  have 


SHOCK  WAVc 


assumed  that  the  shock  wave  and  the  body  surface  are  parallel  to  each  other 
in  the  region  of  interest,  the  thickness  of  the  shock  layer  will  be  equal  to  A 
everywhere  in  the  region  of  the  nose. 


For  a  spher.'cal  nose,  the  dissociation-energy  integral  can  be  written 
in  the  convenient  form 


x„  = 


?  f*,A 

pJbp  (2  fTsl)  d-A.  =  I  pApfe 


Tf/t. 


)  oLa. 


>4. 


(VI.  75) 


where  P  is  the  sonic  point  and  Q  is  the  point  on  the  shock  wave  radially 
opposite  to  P  .  By  assuming  that  the  quantity  pX p  has  a  suitably  defined 
average  over  the  shock  layer  between  P  and  Q  ,  Eq.  VI.  75  can  be  reduced 
to 


TL^-TT  pXD  A  (d+A)  ■ 


Eqs.  VI.  54  and  VI.  76  can  be  combined  together  to  yield 


/°*D 


(VI.  76) 


(VI.  77) 
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Now,  for  the  conditions  of  present  interest  the  dissociation  fraction 

2 


at  Q  is  usually  quite  small;  in  other  words,  Xs  /um  ^  ^ 
accuracy  of  the  present  analysis  we  can  neglect  (  J,s  /  u ^ 
to  unity  and  write  Eq.  VI.  77  as 


Ji  =  J>- 


P  ul 


d/2 


(z*£) 


Within  the 
)  compared 


(VL  78} 


where  -<Cr>p  represents  the  dissociation  energy  in  the  body  streamtube 

and  is  a  constant  downstream  of  the  sonic  point,  P  . 


The  method  of  calculating  ot p  (the  dissociation  fraction  at  the  sonic 
point)  has  been  described  earlier.  Essentially  one  assumes  that  the  pressure 
in  the  nose  region  is  insensitive  to  real-gas  effects,  so  that  the  pressure  at 
the  sonic  point  may  be  written  as  ■p'p  =  0.526  p*  t  •  Since  the  entropy 
along  the  body  streamline  is  also  known,  all  other  state  variables  at  the  sonic 
point  can  be  obtained  immediately  from  a  Mollier  diagram.  However,  the 
value  of  cannot  be  obtained  in  a  simple  manner  from  the  value  of 

since  air  cannot  be  represented  by  a  single,  constant  dissociation  energy. 

That  is,  if  one  writes 


-  d  D  , 


(VI.  79) 


then  D  is  not  a  characteristic  constant  for  air  (as  it  is  for  a  simple  ideal 
dissociating  gas,  for  example). 

The  value  of  can  be  obtained  from  the  energy  equation  along  the 

body  streamline.  The  energy  equation  can  be  written  as 


*  One  of  the  common  definitions  of  the  outer  edge  of  the  entropy  layer  is 
that  it  is  the  streamline  for  which  (d  Ys  /  d  /  at  the  shock. 

Also, for  a  sphere  the  sonic  point  lies  between  40*  and  50°  from  the  axis. 
Therefore,  it  can  be  seen  that  this  approximation  is  consistent  with  our 
earlier  assumption  that  real-gas  effects  are  negligible  outside  the  entropy 
layer. 
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(VL  80) 


T*uC  iV 

and 


J-w  is  the  total  enihalov  of  the  flow.  At  the  sonic  point,  p  =  0.528 p- 
*  .  '  ?  rJf 

u.2  -  iff  Pp/PpJ  •  Thus  Eo-  VL  80  can  be  written  as 


-2 

■1  I 


2.5  M  1+0.2 -0.075Z 
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m. 

rr‘ 
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From  this  the  value  of  the  dissociation  energy  frozen  in  the  body  streamtube 
can  be  calculated  for  any  given  free -stream  conditions. 

The  variations  of  Ji,v  Jua  and  aLb/u2  with  the  free- 
stream  velocity  are  shown  in  Fig.  27  for  an  altitude  of  250,  000  ft.  It  can  be 
seen  from  this  figure  that  a  significant  fraction  of  the  free- stream  kinetic 
energy  remains  frozen  along  the  body  stream  tube.  It  can  also  be  seen  that 
the  value  of  D  varies  significantly  over  the  range  of  free-stream  conditions 
of  interest.  The  dip  that  is  observed  in  both  the  curves  around  a  f~ee~stream 
velocity  of  16,  000  ft/ sec  is  due  to  more  or  less  complete  dissoci_cion  of 
oxygen,  while  at  this  point  nitrogen  dissociation  is  still  relatively  unimportant. 

It  can  be  seen  from  Eq.  VI.  76  that  the  value  of  the  shock  stand-off 
distance  has  also  to  be  known  in  order  to  evaluate  the  frozen-dissociation- 
energy  integral.  In  ideal-gas  flows  various  approximate  methods  (such  as 
the  constant-density  and  thin-shock-layer  approximations)  exist  for  deter¬ 
mining  the  shock  stand-off  distances  on  blunt  bodies.  However,  simple 
formulae  for  the  shock  stand-off  distance  do  not  exist  when  real-gas  effects 
are  present.  In  the  present  analysis  we  will  use  the  expression  for  the  shock 


*  The  subscript  b  has  been  used  earlier  to  denote  values  on  the  afterbody. 
Note  that  since  oc  p  and  J^Dp  are  both  constant  downstream  of  the  sonic 
point,  the  subscripts  P  and  b  can  be  used  interchangeably  for  these 
variables. 
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stand-off  distance  giver,  by  the  constant-density  approximation  (for  the  ideai- 

gas  case),  with  the  density  ratio  being  interpreted  as  that  corresponding  to 

an  equilibrium  shock  wave.  The  consistency  between  the  constant  density 

.  .  -  54 

and  Newtonian  approximations  has  been  pointed  out  by  Hayes  ana  Probstem 
Thus,  since  the  Newtonian  approximation  has  been  used  in  Section  C  as  the 


basis  for  the  ordering  of  the  va-ious  quantities,  it  is  consister.I  to  use  the 
constant  density  approximation  «.o  evaluate  the  shock  stand-off  distance  in 
Eq.  VI.  78.  Therefore,  we  will  assume  that 


4 

d/2 


JL  L.  (8  ] 

's  L  v3  pJ  J 


(VI.  82) 


With  the  heln  of  Eos.  VI.  81  and  82,  the  quantity  J>‘ /Ji  can  now  be 
calculated  from  Eq.  VI.  78.  The  variations  of  the  quantities 
c,(=  4'fjt)  ,  c2  (=  lzl-/£eyc  IkL and  c  (=  I~ITI~2)  with  the 

free-stream  velocity  are  shown  in  Fig.  28  for  a  spherical  r.cse  for  an  altitude 
of  250,  000  ft.  It  can  be  seen  that  the  opposing  behaviors  of  and  are 
such  that  the  value  of  c  remains  close  to  unity  over  the  range  of  free-stream 
velocities  considered.  Clearly  the  length  scale  for  the  downstream  influence 
of  the  blunt  nose  (as  defined  by  Eq.  VI.  74)  is  reduced  by  real-gas  effects  by 
less  than  ten  percent  over  the  range  of  free-stream  conditions  under  consider¬ 
ation. 


The  pressure  distribution  and  shock  shape  (as  given  by  Eqs.  VI.  71  and  56) 
corresponding  to  the  conditions  Up  =  25,  000  f.  p.  s.  ,  h  =  250,  000  ft  and  = 
are  shown  in  Figs.  29  and  30.  The  value  of  c  for  these  conditions  was  very 
close  to  unity,  and  the  real-gas  distributions  are  indistinguishable  from  the 
ideal-gas  ones.  For  the  sake  of  comparison  the  results  obtained  by  the  com¬ 
puter  program  of  Ref.  40  are  also  shown  on  Figs.  29  and  30.  Even  in  the 
numerical  calculations  the  ideal  and  nonequilibrium  flow  solutions  are  found 
to  be  almost  identical. 


Before  proceeding  further  let  us  consider  the  two  other  approximations 
that  were  suggested  for  simplifying  Eq.  VI.  57.  The  only  influence  on  the 
above  results  of  using  the  approximation  suggested  in  Eq.  VI.  64  in  place  of 
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the  one  given  in  Eq.  VI.  65  is  to  slightly  modify  the  value  of  the  parameter  . 
It  is  found  that  the  difference  in  the  value  of  cz  predicted  by  the  two  approx¬ 
imations  is  in  general  less  than  one  percent.  Thus  even  the  crude  approxima¬ 
tion  that  is  represented  by  Eq.  VI.  64  gives  fairly  accurate  results. 


The  other  approximation  that  was  suggested  is  represented  by 
Eqs.  VI.  59  and  62.  Equation  VI.  62  can  be  simplified  in  a  manner  similar 
to  that  used  in  simplifying  Eq.  VI.  75,  with  the  result  one  car.  write 


(VI.  83) 


Again,  the  only  effect  of  using  the  above  approximation  (instead  of  that  given 
in  Eq.  VI.  65)  will  be  to  modify  the  value  of  Cz  .  However,  in  the  present 
case  c.  will  be  a  function  of  and  Eqs.  VI.  83  and  VI.  71  have  to  be  solved 
simultaneously.  On  the  other  hand,  it  does  not  appear  that  this  procedure 
will  drastically  modify  the  qualitative  conclusion  (namely,  that  real-gas  effects 
do  not  markedly  influence  the  pressure  distribution  and  shock  shape)  reached 
above.  Therefore,  the  rather  tedious  procedure  involved  in  carrying  out  the 
simultaneous  solution  of  Eqs.  VI.  83  and  VI.  71  has  not  been  attempted  here. 

As  mentioned  earlier,  the  procedure  used  above  to  account  for  real- 
gas  effects  can  also  be  used  in  conjunction  with  Chernyi's  method.  In 
Chernyi's  method  no  explicit  account  is  taken  of  the  entropy  layer,  the  major 
assumption  being  that  a  major  bulk  of  the  gas  that  crosses  the  shock  wave  is 
contained  in  a  thin  layer  next  to  the  shock  envelope.  The  energy  and  momen¬ 
tum-integral  equations  are  written  in  the  form 


Equations  VI.  84  and  85  determine  the  shock  shape  and  the  afterbody 
pressure  distribution.  The  integral  on  the  right-hand  side  of  Eq.  VI.  84  can 
be  determined  using  Eq.  VI.  65.  Thus  even  within  the  framework  of  Chernyi's 


86 


analysis  real-gas  effects  can  be  accounted  for  in  a  fairly  simple  manner.  The 
details  of  the  analysis  are  straightforward  and  will  not  be  shown  here. 

Now,  it  should  be  pointed  out  that  the  accuracy  of  integral  methods 
depends  quite  critically  on  the  appropriate  choice  of  the  profiles  chosen  for 
representing  the  distribution  of  the  flow  quantities.  The  success  of  integral 
methods  in  boundary  layer  theory  is,  to  a  large  measure,  due  to  the  appro¬ 
priate  choice  of  profiles.  For  example,  the  demonstrated  success  of  the 
95 

Lees-R.eeves  method  in  treating  separated  flows  is  partly  due  to  the  use 

96 

of  the  Cohen-Reshotko  similar  family  of  profiles.  Therefore,  even  in  the 
present  case  an  integral  method  based  on  the  similar  solution  profiles  for 
power  law  piston  motions^’ ^  would  seem  to  hold  promise. 


The  possibility  of  using  similar  solution  profiles  with  integral  methods 

78  52 

has  been  suggested  by  Chernyi  and  Mirels"*  for  the  ideal-gas  problem. 
However,  these  suggestions  have  not  been  exploited  so  far.  A  possible 
approach  to  the  study  of  the  combined  effects  of  nonequilibrium  and  bluntness 
on  the  flows  over  slender  bodies  using  integral  methods  based  on  similar 
solution  profiles  will  be  described  below.  Like  in  the  integral  methods  of 
boundary  layer  theory,  one  can  assume  for  the  present  case  that  at  any  given 
axial  station  the  distributions  of  all  the  flow  variables  normal  to  the  body  are 
given  by  the  similar  solution  profiles  ’  /  corresponding  to  some  unspec¬ 

ified  values  of  the  power-law  index  and  the  specific -heat  ratio.  The  power- 
law  index  and  specific-heat  ratio  characterizing  the  profiles  are  of  course 
assumed  to  vary  from  one  axial  station  to  another,  this  variation  itself  being 
a  part  of  the  solution  sought.  This  approach  is  equivalent  to  describing  the 
flow  at  each  axial  station  with  an  average  value  of  the  specific -heat  ratio, 
with  the  average  specific -heat  ratio  being  different  for  different  axial  positions. 
The  variation  of  the  power-law  index  and  the  average  specific -heat  ratio  with 
axial  position  can  be  determined  from  the  various  conservation  equations  and 
the  shock  relations. 


Of  course,  like  in  all  other  integral  methods  the  above  approach  would 
involve  considerable  numerical  calculations.  However,  incorporation  of  some 
of  the  assumptions  of  thin  shock  layer  theory  into  the  above  scheme  might  lead 
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to  useful  and  tractable  solutions.  As  mentioned  earlier,  the  use  of  the  similar 
solution  profiles  to  describe  the  distribution  of  at  least  some  of  the  flow  quan¬ 
tities  will  lead  to  considerably  better  results  than  using  any  other  simple 
profiles.  An  investigation  of  this  area  would  be  worthwhile. 

F.  VISCOUS  EFFECTS 


We  will  now  consider  the  combined  effects  of  nose  bluntness,  boundary 

layer  displacement  and  flow  nonequilibrium  on  slender-body  flow  fields. 

Viscous  effects  can  be  incorporated  into  the  inviscid  analyses  given  above 

31 

by  the  model  proposed  by  Cheng  .  In  this  model  a  thin  and  distinct  boundary 

layer  is  assumed  to  exist  submerged  within  the  entropy  layer.  In  particular, 

the  eventual  'swallowing'  of  the  entropy  layer  by  the  boundary  layer  (which  is 

known  to  occur  at  high  altitudes)  is  not  accounted  for.  An  accounting  of  this 

phenomenon  does  not  seem  possible  wituin  the  scope  of  integral  methods.  In 

31 

addition  to  the  assumptions  made  by  Cheng  ,  in  the  present  case  it  has  been 
assumed  that  the  chemical  composition  inside  the  boundary  layer  remains 
frozen. 


We  will  first  consider  the  similitude  for  blunt-nosed  slender  bodies 

when  both  viscous  and  real-gas  effects  are  present.  The  viscous  hypersonic 

similitude  for  slender  bodies  without  nose  biuntness  has  been  given  by  Hayes 
30  .  87 

and  Probstein  and  Luniev  ,  while  the  similitude  including  the  effects  of 

,,  31  88 

nose  bluntness  has  been  given  by  Cheng  and  Luniev  .  The  above  similitudes 

show  that  for  hypersonic  viscous  flows  past  slender  bodies  which  have  similar 

shapes  and  surface -temperature  distributions  of  the  form 


"  |  -  >(B 


*  In  the  present  section,  the  subscript  denotes  quantities  ,  the  wall, 

while  the  subscript  b  denotes  the  inner  edge  of  the  entropy  I  /er. 
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The  similitude  parameters  {in  addition  to  those  that  are  relevant  to  the 
corresponding  inviscid-flow  case)  are  the  prandtl  number,  Pr  ,  and  the 
quantity  V  Re^  .where 

P  U  L 

=  -----  •  (VI.  87) 

'  XV 


The  quantity  has  to  be  evaluated  at  some  suitable  reference  point  on 

the  body. 

The  requirement  of  the  equality  of  the  parameter  T  Re  ^  for 
similarity  between  two  flows  is  merely  a  statement  of  the  condition  that  the 
distribution  of  the  displacement  thicknesses  on  the  two  (affinely  related) 
bodies  must  follow  the  same  affine  law  as  their  thickness  distributions. 
Thus,  the  basic  requirement  for  similarity  can  be  stated  as 


Jl 

rL 


(VI.  88) 


For  a  linear  viscosity-temperature  relation  of  the  form 

/'v  _  £  Tar 

To  ' 


(VI.  89) 


the  requirement  expressed  by  Eq.  VI.  88  reduces  to  the  condition  that  the 
parameter 


-  _ 

71 

'  a o 


(VI.  90) 


tO 

must  be  an  invarient,  while  for  a  general  viscosity  law  of  the  form  M  rj~T  , 

ZrO  /  - 

the  condition  is  that  the  quantity  M  /  yRe  must  be  invarient. 


With  the  above  remarks  in  mind  we  now  turn  our  attention  to  the 
extension  of  the  similitudes  given  in  Section  VI.  D  to  include  viscous  ef.tects. 

Now,  if  the  flow  in  the  nose  region  of  a  blunted  slender  body  is  required  to 
be  simulated,  then  similarity  can  be  obtained  by  means  of  the  Mach-number- 
independence  principle.  As  pointed  out  earlier,  the  Mach-number- independence 
principle  is  applicable  even  when  viscous  effects  are  included;  however,  now 
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the  appropriate  length  scale  of  th*»  nose  has  to  be  duplicated  in  addition  to  the 
free-stream  density  and  velocity.  Thus,  for  the  nose  region  one  can  write 


When  the  parameters  on  the  right-hand  side  of  Eq.  VI.  91  are  the  same  for 
two  flows,  then  .hey  are  not  merely  similar  but  essentially  identical. 

Now,  the  solution  of  a  hypersonic  boundary  layer  is  completely  deter¬ 
mined  if  the  appropriate  initial  conditions,  the  wall  conditions  and  the  quantities 
,  H[,  (-  h'a, )  ,  and  U  b  at  the  outer  edge  of  the  boundary  layer  are 

duplicated?  The  appropriate  initial  conditions  are,  of  course,  given  by 
Eq.  VI.  91  while  the  requirements  at  the  wall  are  the  no- slip  condition  and 

the  wall-temperature  specification.  Moreover,  under  the  hypersonic  small 

/  z 

disturbance  approximation,  ufe  Uco  and  H<o  ~  T  ■ 

Thus  in  a  manner  completely  analogous  to  Eq.  VI.  8,  the  evolution  of  the 
boundary  layer  displacement  thickness  can  be  written  in  the  form 
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(VL  92) 


v  For  a  blunt-nosed  slender  body  the  subscript  b  refers  to  tne  outer  edge  of 
the  boundary  layer  which  is  the  same  as  the  inner  edge  of  the  entropy  layer. 
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where  Tq  is  the  stagnation  temperature  of  the  flow. 

Under  the  hypersonic  small  disturbance  approximation,  Eq. 
reduces  to  the  form 


VI.  92 


(VI.  93) 


In  the  above  derivation,  consistent  with  the  classical  boundary  layer  approach, 
we  have  assumed  that  the  pressure  distribution  can  be  specified  arbitrarily 
and  that  it  is  part  of  the  conditions  that  are  specified.  However,  in  the  present 
case  there  is  an  interaction  between  the  boundary-layer  growth  and  the  pressure 
field,  and  it  is  precisely  this  interaction  that  we  are  seeking  to  simulate. 


Clearly,  the  pressure  distribution  at  the  outer  edge  of  the  boundary 
layer  of  a  blunt-nosed  slender  body  will  be  determined  by  an  equation  of  the 
form  given  by  Eq.  VI.  11,  but  with  an  effective  body  shape  parameter 
in  place  of  f  .  Here  8^  is  the  value  of  the  displacement  thickness  at  some 
conveniently  chosen  reference  point  (the  mid-chord  of  a  wing,  for  example). 
Thus,  one  can  write 
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For  bodies  with  similar  noses,  the  nose  drag  coefficient  Ji  can  be 
dropped  from  the  similitude  and  Eqs.  VI.  93  and  VI.  94  can  be  combined 
together  to  give 
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Thus,  it  can  be  seen  from  Eq.  VI.  95  that  for  the  hypersonic  viscous  flows 
past  blunt-nosed  slender  bodies  which  obey  Eq.  VI.  86  the  only  additional 
(that  is,  not  included  in  the  corresponding  inviscid  case)  requirement  for 
similarity  is  the  duplication  of  the  nose  length  scale. 

Finally,  the  similitude  relations  for  the  surface  pressure,  shock  shape 
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and  the  heat-transfer  rate  can  be  written  in  the  form 


p  «  *  r 


«.nA 


M  £. 


/  V 

(VL  96} 


In  Eq.  VI.  96  the  boundary  layer  assumption  ^  -5,  has  been  made, 

and  the  parameters  /  ,  sj  and  Pr  (which  are  us-wlly  included  in  general 
viscous  similitudes)  are  not  included  since  we  are  assuming  here  that  only 
flows  of  the  same  gas  are  being  studied;  these  parameters  can  be  included 
if  desired. 

Unfortunately  no  solutions  have  been  available  to  us  for  testing  the 
accuracy  and  the  range  of  validity  of  the  above  similitude.  However,  since 
the  basic  premise  on  which  the  similitude  is  based  (namely,  that  real-gas 
effects  are  important  only  in  that  part  of  the  flow  field  which  is  very  near  the 
body)  has  been  ver  iied  for  the  invisid  case,  one  can  expect  that  the  viscous 
similitude  will  be  valid  over  approximately  the  same  range  of  conditions  as 
the  corresponding  inviscid  similitude. 

Cheng  has  also  developed  a  theory  for  the  combined  influence  of 
boundary-layer  displacement  and  leading  edge  effects  on  slender  body  flow 
fields.  Cheng's  development  i3  based  on  die  observation  that  the  local  flat- 
plate  similarity  approximation  is  consistent  with  thin-shock-layer  concepts. 
Starting  with  the  above  observation,  Cheng  included  boundary  layer  displace¬ 
ment  effects  into  Eq.  VI.  57  (for  the  plane  flow  case,  -  0  )  by  replacing 
Ytr  by  Yw  where 


Y«.  =  Yfc  +  S* 


(VI.  97) 


where  the  displacement  thickness  '  S*  is  related  to  the  wall  temperature  and 
the  free-stream  conditions  through  the  Newton-Busemann  pressure  formula 
given  in  Eq.  VI.  56.  Cheng's  analysis  is  valia  only  for  plane  flows  since 
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transverse  curvature  effects  (which  are  often  dominant  in  axisymmetric  flows) 
are  not  incl»ded.  Since  we  are  mainly  concerned  with  axisymmetric  flows, 
we  have  not  pursued  the  generalization  of  Cheng's  viscous  flow  case  to  include 
real-gas  effects. 

G.  SUMMARY  AND  DISCUSSION  OF  THE  RESULTS 

The  combined  effects  of  nose  bluntness  and  chemical  nonequilibrium 
phenomena  on  slender  body  flow  fields  was  studied  by  using  a  simplified  model 
for  flow  chemistry.  The  model  was  based  on  the  observation  that  over  a 
significant  range  of  free-stream  conditions  real-gas  effects  are  more  impor¬ 
tant  within  the  high  entropy  layer  near  the  body  than  outside  it.  Thus  in  the 
model  that  was  used,  real-gas  effects  were  neglected  in  the  thin  layer  close 
to  the  shock  wave  (and  outside  the  entropy  layer),  and  the  gas  behavior  in 
this  region  was  assumed  to  be  ideal  with  a  constant  specific  heat  ratio  equal 
to  that  of  the  free  stream. 

It  was  also  assumed  that  the  flow  in  the  nose  region  was  effectively 
in  equilibrium  and  that  the  flow  freezes  aloig  a  line  passing  through  the  sonic 
point  on  the  body  due  to  the  rapid  expansion  in  the  region  of  nose-afterbody 
junctions.  Under  this  approximation,  the  entropy  layer  is  characterized  by 
different  values  of  entropy  and  effective  (frozen)  specific -heat  ratio  along 
each  streamline. 

Based  on  the  above  model,  the  effects  of  real-gas  phenomena  on 

blunted  slender  body  flow  fields  was  studied  using  a  zeroth-order  thin-shock- 

layer  theory;  that  is,  quantities  of  the  order  of  and  (M  &)  were 

2f+l 

assumed  to  be  negligible.  The  analysis  was  essentially  an  extension  of  the 
theory  developed  by  Cheng ^  for  the  corresponding  ideal-gas  case.  It  was 
found  that  the  effect  of  real -gas  phenomena  on  the  downstream  influence  of 
a  blunt  nose  can  be  throught  of  as  causing  a  net  apparent  reduction  in  the  nose 
diameter. 

The  real-gas  correction  factor  for  the  nose  diameter  was  found  to  be 
influenced  by  two  different  effects,  one  due  to  the  smaller  amounts  of  momen¬ 
tum  and  energy  that  are  imparted  to  the  transverse  flow  field  in  the  real-gas 
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case  than  in  the  corresponding  ideal-gas  case,  and  the  other  due  to  an  increased 
value  of  the  average  specific  heat  ratio  within  the  entropy  layer.  The  first 
effect  tends  to  reduce  the  downstream  influence  of  the  blunt  nose  while  the 
second  effect  tends  to  increase  the  downstream  influence.  Therefore,  the 
value  of  the  correction  factor  remains  close  to  unity  over  the  range  of  free- 
stream  conditions  of  interest,  and  the  shock  shape  and  the  pressure  distribu¬ 
tion  on  the  body  were  found  to  be  influenced  very  little  by  the  presence  of 
real-gas  effects. 

To  check  the  above  conclusion,  numerical  computations  were  carried 

out  for  both  ideal  and  real-gas  flows  over  blunted  slender  cones  by  using  the 

40 

computer  program  developed  by  Curtis  and  Strom  .  The  numerical  results 
confirmed  that,  for  the  free-stream  conditions  under  consideration,  the  effects 
of  real-gas  phenomena  on  the  surface -pressure  distributions  and  shock  shapes 
on  blunted  slender  bodies  are  indeed  quite  small. 

51 

The  above  conclusion  is  in  contrast  to  the  result  obtained  by  Whalen 
who  found  that  the  pressure  distributions  on  wedge-nosed  flat  plates  were 
significantly  influenced  by  real-gas  effects.  In  fact,  Whalen  .concluded  that 
the  blast-wave  phenomena  may  be  essentially  eliminated  when  real-gas  effects 
are  present.  A  somewhat  analogous  conclusion  was  reached  in  Refs.  1  and  2 
in  connection  with  the  study  of  the  sensitivities  of  afterbody  flow  quantities  on 
wedge-nosed  flat  plates  to  mismatches  in  free-stream  conditions.  In  Refs.  1 
and  2  it  was  found  that  the  pressure  immediately  downstream  of  the  nose- 
afterbody  junction  (of  a  wedge-nosed  flat  plate)  was  considerably  smaller  than 
the  corresponding  ideal-gas  value,  so  that  it  would  appear  that  in  this  case 
the  downstream  influence  of  the  blunt  nose  would  be  significantly  reduced  by 
the  presence  of  real-gas  effects. 

Therefore,  we  can  conclude  that  the  effect  of  real-gas  phenomena  on 
the  downstream  influence  of  a  blunt  nose  will  depend  on  the  shape  of  the  nose. 
This  is  not  surprising  since  real-gas  effects  are  essentially  nonsimilar,  and 
the  flow  downstream  of  the  nose  cannot  be  expected  to  be  independent  of  the 
details  of  the  flow  in  the  nose  region  when  real-gas  effects  are  present.  In 
general,  it  appears  that  when  there  is  a  smooth  junction  between  the  nose  ar._ 
the  afterbody  (such  as  in  a  spherically  blunted  cone),  real-gas  effects  do  not 
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s ignificantly  influence  either  the  surface  pressure  distribution  or  the  shock 
shape.  On  the  other  hand,  when  the  junction  between  the  nose  and  the  after¬ 
body  is  a  sharp  corner  (such  as  in  a  wedge-nosed  flat  plate  or  a  flat-plate 
slender  body),  it  appears  that  real-gas  effects  will  significantly  reduce  the 
downstream  influence  of  the  blunt  nose,  and  in  some  cases  altogether  eliminate 
the  blast-wave  effect. 

A  restricted  similitude  applicable  to  real-gas  flows  past  blunt-nosed 

slender  bodies  was  also  derived  using  the  concept  that  over  a  significant 

range  of  free-stream  conditions  real-gas  effects  are  important  only  within 

the  entropy  layer  and  not  outside  it.  It  was  found  that  the  requirement  of  the 

duplication  of  the  free-stream  temperature  (which  is  part  of  the  similitudes 

25  26 

given  by  Cheng  and  Inger  )  can  be  relaxed  provided  that  the  nose-shape 

and  the  free-stream  velocity  are  duplicated.  This  similitude  is  especially 

useful  when  similitude  is  required  not  only  in  the  region  of  the  afterbody  but 

also  in  the  nose  region.  For  the  above  case  the  requirements  imposed  by  the 

25  26 

similitudes  given  by  Cheng  and  Inger  would  imply  that  there  is  essentially 
no  similitude;  that  is  all  conditions  have  to  be  duplicated.  The  similitude  is 
found  to  be  valid  for  the  viscous  case  also  provided  that  the  characteristic 
length  scale  of  the  nose  is  duplicated. 
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SECTION  VII 
CONCLUSIONS 


In  the  present  study,  we  have  examined  the  similitude  requirements 
for  hypervelocity,  slender-body  flows  when  nonequilibrium  and  nose-bluntness 
effects  are  present.  For  such  flows,  it  has  been  shown  by  Cheng  and  Inger 
that  proper  similitude  requires,  among  other  things,  complete  duplication  of 
the  free-stream  thermo-chemical  state.  However,  in  many  cases,  the  res¬ 
ervoir  conditions  required  to  duplicate  all  the  similitude  parameters  simul¬ 
taneously  in  a  wind  tunnel  are  so  severe  as  to  make  partial  simulation  the 
only  recourse.  Therefore  the  present  study  has  been  concerned  primarily 
with  the  effects  of  mismatches  in  certain  free-stream  conditions  on  slender- 
body  flow  fields. 

Two  related  studies  have  been  carried  out.  In  the  first  study,  direct 
comparisons  between  the  flight  and  wind-tunnel  flow  fields  were  made  for  a 
blunted  slender  cone  by  assuming  a  hypothetical  wind-tunnel  performance. 

The  pressure  distribution  on  the  cone  was  assumed  to  be  uninfluenced  by 
real-gas  effects,  and  the  influence  of  rsal-gas  effects  on  the  other  flow  quan¬ 
tities  was  calculated  by  using  a  streamtube  computer  program.  The  second 
study  consisted  of  an  analytical  investigation  of  real-gas  effects  on  blunted 
slender  body  flows  by  using  a  thin-shock-layer  approach.  This  study  was 
based  on  a  simplified  model  for  the  flow  chemistry  in  which  real- gas  effects 
were  included  in  the  high  entropy,  low  density  layer  near  the  body  but  was 
neglected  in  the  high  density,  lower  entropy  layer  near  the  shock  wave. 

In  the  first  study  it  was  found  that,  at  the  free-stream  conditions  under 
consideration,  the  flow  in  the  vicinity  of  the  blunt  nose  is  Mach  number  inde¬ 
pendent  even  when  real-gas  effects  are  present.  Thus  the  flow  within  the 
entropy  layer  (that  is,  the  layer  consisting  of  streamlines  that  have  crossed 
the  nearly  normal  portions  of  the  shock  wave)  also  remains  Mach  number 
independent.  On  the  other  hand,  the  flow  in  the  outer  portion  of  the  shock 
layer  (which  consists  of  streamlines  that  have  crossed  the  more  oblique  parts 
of  the  shock  wave)  is  more  sensitive  to  mismatches  in  free-stream  temperature, 
and  is  analogous  to  the  sharp-nosed  slender  body  case.  Thus  the  amount  of 
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nose  bluntness  has  an  important  bearing  on  the  similitude  characteristics 
since  it  determines  the  relative  extents  of  the  entropy  layer  and  the  outer 
regions  of  the  shock  layer.  In  the  spherically  blunted  slender  cone  cases 
considered  for  =  1.  5  in. ,  the  flow  along  the  cone  is  dominated  by  the 
blunt  nose  and,  hence,  is  insensitive  to  mismatches  in  Moo  •  However,  for 
R^-  =  0.  2  in,  ,  the  flow  over  most  of  the  conical  afterbody  will  behave  as  a 
sharp-nosed  body  and  mismatches  in  M „  can  be  important.  In  both  cases 
the  effects  of  nonequilibrium  phenomena  are  as  large  or  larger  than  the 
differences  between  the  flight  and  wind-tunnel  flow  fields.  Thus,  it  is  impor¬ 
tant  to  match  body  scale  as  well  as  (Jco  sc  that  L/Uto  is  duplicated. 

In  the  second  study,  which  was  based  on  thin-shock-layer  concepts, 
it  was  found  that  the  effect  of  real-gas  phenomena  on  the  downstream  influ¬ 
ence  of  a  blunt  nose  can  be  thought  of  as  causing  an  apparent  reduction  in  the 
diameter  of  the  nose.  However,  it  was  found  that  the  effect  was  fairly  small 
for  spherically  blunted  cones  over  the  range  of  free-stream  conditions  that 
were  consider  ’d.  However,  it  appears  that  the  effect  of  real-gas  phenomena 
on  blunted  slender  body  flows  will  be  a  function  of  the  details  of  the  nose  shape. 
Thus,  it  is  likely  that  real- gas  effects  would  be  quite  significant  for  bodies 
with  flat  noses,  or  where  there  is  a  sharp  corner  at  the  nose-afterbody 
junction. 

A  new  similitude  applicable  to  real- gas  flows  over  blunted  slender 
bodies  was  also  derived.  This  similitude  has  the  advantage  that  the  duplication 
of  the  free-stream  temperature  (which  is  required  by  existing  real-gas  sim¬ 
ilitudes)  has  been  relaxed.  The  extension  of  the  similitude  to  include  viscous 
effects  is  also  given. 

The  present  study  serves  to  illustrate  some  of  the  problems  associated 
with  partial  simulation  of  the  hypervelocity  flows.  An  analysis  of  the  type 
reported  here  should  be  made  for  each  hypervelocity  wind  tunnel  in  order  to 
delineate  the  flight  conditions  that  can  be  simulated  in  a  given  facility.  In 
this  regard  it  should  be  reiterated  that  the  hypothetical  wind-tunnel  perform¬ 
ance  assumed  here  far  exceeds  that  of  existing  facilities,  having  been  selected 
in  order  to  minimize  nozzle-flow  nonequilibrium  effects. 
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